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Two of the bands in the spectrum of H2S identified as 
v3 and » have been remeasured under high dispersion. 
The center of the band identified as v3 is found to be 
2685 cm~ and is a type of band showing no single Q branch 
at the center and is therefore ascribed to an oscillation of 
the electric movement along the axis of the intermediate 
moment of inertia. The fundamental », has been remeas- 
ured with two different gratings of different resolving 
power. It has been found to be of the type displaying a 


broad and partially resolved Q branch at the center and is 
henceforth identified with an oscillation of the electric 
moment along the axis of smallest moment of inertia. An 
estimate of the values of the moments of inertia has been 
made and these are found to agree well with those given 
by Cross, but the data here presented is interpreted as 
evidence for a model where the vertex angle is 85° rather 
than 92° as given by Cross. 





I. INTRODUCTION 


HE absorption in the infrared by the hydro- 

gen sulfide molecule, presumably a non- 
linear triatomic molecule similar to water vapor, 
has been studied by a number of investigators.! 
Absorption bands, enumerated in the order of 
their intensities, have been located at approxi- 
mately the following frequency positions: 1260 
em, 3795 cm=, 2685 cm-!, 9911 cm, 10,100 
cm, and 2615 cm. All of these except the last, 
which has been found only in the Raman spec- 
trum, have been measured under high dispersion 
and have been given the following identification, 
using Dennison’s notation :? ve, ve+v3, v3, 3v1 +73, 
3v3+v1 and »;, where the frequency »; is perpen- 
dicular to the bisector of the apex angle and 
where »; and v2 are along it. In general the 
frequencies y; and v3 take values relatively high 
1 A.H. Rollefson, Phys. Rev. 34, 604 (1929) ; W. Mischke, 
Zeits. f. Physik 67, 106 (1931); H. H. Nielsen and E. F. 
Barker, Phys. Rev. 37, 727 (1931); H. H. Nielsen and A, 
D. Sprague, Phys. Rev. 37, 1183 (1931); P. C. Cross, Phys. 


Rev. 47, 7 (1935). 
*-D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 


85 


compared with v2. From an analysis of the 
rotational structure of the high frequency combi- 
nation bands at 9911 and 10,100 cm Cross has 
arrived at certain values of the moments of 
inertia and has inferred that the vertex angle 
must be slightly in excess of 92°. This would mean 
that v3; would lie along the axis of least moment of 
inertia while vy; and v2 would lie along the axis of 
the intermediate moment of inertia. According 
to the charts prepared by Dennison on the 
appearance of infrared bands of asymmetrical 
type molecules, v; should then present a rota- 
tional structure containing a rather broad Q 
branch while v; and v2 should be expected to 
display an absence of such a Q branch. The band 
taken to be the fundamental v3 is one of the bands 
measured by Nielsen and Barker and does not 
appear to give evidence of any Q branch what- 
ever. It has been suggested that this band must 
have another half which might lie to higher or to 
lower frequencies. This belief was somewhat 
further strengthened by the measurements made 
by Rollefson on the band near 8.0 interpreted as 
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ve and which gave evidence of a minimum near 
the center. It has therefore seemed important to 
remeasure these two absorption regions and in 
the 3.74 region to extend the measurements well 
beyond the limits at both ends of the region as 
shown in the work of Nielsen and Barker. This 
has been the purpose of this investigation. 


II. EXPERIMENTAL WoRK 


The hydrogen sulfide gas used in making these 
measurements was obtained in tanks from the 
Ohio Chemical and Manufacturing Company in 
Cleveland. To insure that the gas was free from 
water vapor and carbon dioxide the gas was 
collected in a trap surrounded by liquid air and 
then fractionally distilled into the evacuated 
absorption cell. For both regions the same cell 
was used. It was made from a glass tube five 
centimeters in diameter and thirty centimeters 
long closed at the ends by windows of polished 
rocksalt. 

For making the measurements a spectrometer 
of the usual design was available. Three echellette 
gratings ruled by R. W. Wood were used, one 
with 3600 lines per inch for the 3.74 region and 
two others with 800 lines per inch and 2000 lines 
per inch for the region near 8.0u. It was found 
feasible in each case to use slit widths equivalent 
to a spectral breadth of about 0.5 cm and 
deflections were taken along the circle at these 
intervals. 


III. EXPERIMENTAL RESULTs. THE 3.74 REGION 


This is one of the two regions originally 
mapped by Nielsen and Barker. This region has 
here been carefully remeasured from 3.574 to 
3.83u and these measurements extend beyond the 


region measured by Nielsen and Barker at each 
end by about 50 cm~". In none of the sets of data 
taken could anything be found in addition to 
that shown in the original work by these authors. 
Bailey, Thompson and Hale, in a recent article* 
have reported an absorption maximum at a 
slightly shorter wave-length than this and have 
suggested that the —7 line reported in Nielsen 
and Barker’s work is to be interpreted to be a 
weak Q branch. This view we can definitely not 
ascribe to since then what we have observed 
must be only one branch. As we have already 
pointed out we can find no evidence of any other 
half of this band and what further absorption 
that may be found on either side of what we have 
measured is certainly far less intense and can 
hardly be identified as another branch. Moreover 
it would seem extremely difficult to account for a 
Q branch having a relative intensity as small as 
that of line —7 reported by Nielsen and Barker. 
Moreover no evidence of the frequency 2615 


TABLE I. The 8.0u region under low dispersion. Wave 
numbers and wave-lengths. 

















FREQUENCIES FREQUENCIES 

LINE No. IN cm™} LINE No. IN cmM™ 
1 1196.8 14 79.6 
2 1209.6 15 87.5 
3 21.1 16 93.8 
+ 26.3 17 1304.5 
5 35.9 18 16.8 
6 40.3 19 30.4 
7 44.2 20 41.2 
8 48.7 21 48.5 
9 54.1 22 56.4 
10 58.2 23 66.0 
11 62.3 24 74.6 
12 67.8 25 82.3 
13 71.1 26 94.5 

















3C. R. Bailey, J. W. Thompson and J. B. Hale. J. Chem. 


Phys. 4, 625 (1936). 
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cm~! was to be found. The band appears then to 
have its center at about 2685 cm instead of 
2630 cm~ and in our opinion very definitely does 
not show a Q branch. This band must therefore 
be ascribed to an oscillation of the electric 
moment along the axis of the intermediate 
moment of inertia. These measurements are not 
shown graphically since they repeat almost 
exactly those given in the work of Nielsen and 
Barker ; neither has any tabulation of frequency 
and wave-length positions of the lines been 
included, these also being in excellent agreement 
with the values given by them. 


THE 8.04 REGION 


This band, originally reported by Rollefson, 
has here been remeasured under two different 
dispersions, one about equivalent to that avail- 
able to Rollefson and one approximately three 
times as large. Under the lower dispersion the 
lines found in general agree well with the posi- 
tions of lines reported by Rollefson, although in 
many cases the intensities of these are quite 
different. In fact, the appearance of the band is 
an entirely different one, in the present case there 
being no evidence of any minimum near the 
center, but rather the opposite.‘ At what appears 
to be the center of the band there occurs a broad 
line of appreciably greater intensity than any of 
its neighbors. The average spacing between 
principal lines seems to be about 10 cm-, very 

4 In the work by Nielsen and Barker mention was made 
of the fact that in a set of preliminary grating measure- ¢ 
ments made by Professor Barker on this region an absorp- 
tion pattern was obtained, but in subsequent trials nothing 
was to be found. The above measurements have been 
compared with the aforementioned run and the two have 
been found to agree almost exactly with each other. The 
fact that nothing was to be found in the subsequent runs 
referred to in this earlier work can be interpreted only as 
meaning that the monochrometer preceding the spectrom- 


eter was out of proper adjustment while these measure- 
ments were made. 





nearly that found between lines in the band at 
2.74 by Nielsen and Barker and by Cross at 
10,100 cm-! between lines arising from his 
J;—J;+1 transitions. This spacing appears to 
be about ten percent larger than the spacing 
between lines at 3.74. The absorption pattern 
delineating those characteristics that repeated in 
the many runs made are shown in Fig. 1 and 
the frequencies of the principal lines measured 
under low dispersion are set down in Table I. It 
was found necessary to guard against falsification 
by water vapor in the atmosphere. This was 
accomplished by taking readings with the cell 
containing the gas in the beam, then with the 
evacuated cell in the beam, and comparing. Also 
as a further check, the deflections obtained with 
the evacuated cell in the beam were plotted 
underneath the percentage absorption curve to 
see that there was no consistent agreement 
between the absorption peaks apparently due to 
gas and the absorption peaks due to the atmos- 
phere. The impression which this band leaves 
with us is that it resembles the band measured by 
Nielsen and Barker at 2.74 in its rotational 
characteristics and also those reported by Cross 
much more than it resembles the band at 3.7u. 
As its center we have chosen 1290 cm, the 
position of the most intense of the lines in the 
band which we believe to be a Q branch. 

When this region is examined using the grating 
ruled with 2000 lines per inch a great many new 
details in the absorption pattern become dis- 
cernible. The absorption pattern, which has 
consistently repeated itself in the smallest details, 
is shown in Fig. 2. At a glance the two curves, 
Fig. i and Fig. 2, appear quite dissimilar but 
upon careful comparison of the one to the other, 
the general features of the two will be seen to be 
the same. Some uncertainty concerning the 
positions of some of the lines on the high fre- 
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quency side must be allowed, due again to 
falsification by atmospheric water vapor lines. 
This falsification is much more difficult to 
compensate for here than when the dispersion is 
lower. Table II lists the frequency and wave- 
length positions of the most important peaks in 
the band. 

Also here the most intense part of the band is 
at 1290 cm and more even than before does this 
line appear to be a broad and partly resolved Q 
branch. We shall therefore interpret this band as 
arising from a vibration along the axis of least 
moment of inertia. 


IV. DiscussIoN OF EXPERIMENTAL RESULTS 


The evidence obtained by these measurements 
seems to indicate that some modification of the 
model suggested by Cross, which has as the axis 
of least moment of inertia a line normal to the 
bisector of the apex angle, is necessary. There can 
be little doubt concerning the identification of 
the vibration frequency 1290 cm as v2 since this 
frequency must be a relatively small one. Since, 
however, this band shows a Q branch, v2 must be 
a vibration along the axis of least moment of 
inertia, a condition which can be fulfilled only 
provided the apex angle be less than 88° 16’. 
‘Moreover for the frequency 2685 cm='!, which in 
our opinion definitely shows no Q branch, to be 
v3 it must be required to be a vibration along the 
axis of the intermediate moment of inertia. This 
again indicates an apex angle of less than 88° 16’. 

It is now necessary to reexamine the identifi- 
cation of the other bands known to exist in the 
infrared spectrum of hydrogen sulfide. The 
frequency of any band will in general be given by 
the relation v=.1+n2v2-+n3v3 where, according 
to Dennison, when 7; is odd the electric moment 
will be induced in a direction normal to the 
bisector of the apex angle, but will be along the 
bisector when 7; is even. On the basis of such a 
model as we propose, namely one for which the 


vertex angle is less than 88° 16’, and taking the * 


Raman line 2615 cm~ as »;, 1290 cm~ as v2 and 
2685 cm—‘as vzit is evident that the band centered 
at 3790 cm and showing a Q branch can no 
longer be ve+v3, as suggested by Mecke,® and 
Cross, since with m3=1, an odd integer, the 


5 R. Mecke, Zeits. f. physik. Chemie 16, 431 (1932), 
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vibration would be along the axis of the inter- 
mediate moment of inertia. It would then show 
no Q branch. This band might, however, be set 
equal to vi+v2. Likewise the frequencies 9911 
cm~! and 10,100 cm~ are of the type in which Q 
branches occur and on the basis of our proposed 
model they must be redefined. We suggest that 
they may be set equal to 2v3+2v, and 2v2+3”, 
say, rather than to »3+3v, and 3v3+”1 as 
suggested by Cross. In fact these values agree a 
little better numerically, especially since v3 must 
be taken to be 2685 cm~ rather than 2630 cm™. 
In Table III is summarized the manner in which 
we have identified the observed bands. 

It would be of interest if the actual value of the 
apex angle of the molecule could be computed. 
This could of course be accomplished if the 


TABLE II. The 8.0u region under high dispersion. Wave 
numbers and wave-lengths. 























FREQUENCIES FREQUENCIES 
LINE No. IN cmM™} LINE No. IN cmM™! 
1 1203.9 39 1270.3 
2 04.7. —«, 40 71.9 
3 05.7 41 73.3 
4 08.1 42 76.1 
5 10.3 43 78.1 
6 12.3 44 79.9 
7 13.3 45 81.2 
8 16.1 46 84.3 
9 17.6 47 85.4 
10 19.9 48 86.6 
11 21.9 49 87.5 
12 23.0 50 90.0 
13 24.3 51 91.1 
14 25.3 §2 93.4 
15 27.6 S3 94.4 
16 29.4 54 96.7 
17 32.1 55 97.8 
18 34.4 56 99.8 
19 36.4 57 1302.1 
20 38.2 58 03.6 
21 38.9 59 05.0 
22 41.0 60 05.9 
23 43.2 61 08.2 
24 45.9 62 12.6 
25 47.1 63 15.5 
26 48.4 64 23.8* 
27 49.3 65 27.6 
28 52.5 66 30.2 
29 53.7 67 39.6 
30 54.7 68 41.4 
31 55.7 69 43.9 
32 57.0 70 52.6 
33 59.0 71 55.5 
34 60.9 72 64.6 
35 63.2 73 69.1 
36 64.4 74 74.2 
37 67.0 75 83.3 
38 68.5 











*From this point on, intense atmospheric water vapor absorption 
peaks occur from time to time and for this reason the positions of the 
absorption lines due to H2S as given may indeed be slightly in error. 
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TABLE III. Identification of observed bands. 

















FREQUENCY FREQUENCY 
COMBINATION IN cmM™! COMBINATION IN cM™! 
V1 2615 cm vitve 3970 cm7 
V2 1290 2v3+2r; 9911 
V3 2685 2vo+3r, 10100 











correct values of the moments of inertia could be 
determined. The values J4, Jp, and Je as de- 
termined by Cross from his analysis of the two 
high frequency combination bands seem very 
convincing indeed, also in the light of our data 
even though no exact rotational analysis has 
been made. This may be seen in the following 
manner. Two quite different types of bands have 
been observed, namely that at 3.74 and the other 
to which we believe those at 2.74 and 8.0u 
belong. The separation between the lines is, 
however, nearly the same in both types, about 
9.0 cm and 10 cm™, respectively. If we ap- 
proximate the hydrogen sulfide molecule by a 
rotator where the moments of inertia are 
A=B=C/2 we should expect to find all the 
bands in the spectrum to be of the perpendicular 
type with the spacings between their lines 
identical. From Dennison’s charts, however, it is 
seen that for a vibration along A there will 
develop a crowding together toward the center 
of the component lines in the band to form a Q 
branch as A is made less than B. Moreover, for an 
oscillation along B quite the opposite effect is to 
be observed as A is decreased with respect to B. 
Now clearly the value of Dennison’s parameter 
of asymmetry p defined as A/B cannot be much 
less than 0.9 since at least under not too high 
dispersion some semblance of regularity in the 
rotational structure is maintained. Approxi- 
mating the hydrogen sulfide molecule then by a 
symmetric rotator where A’= B’=1/[(A+B)/2] 
and C’=2A’ and remembering that prominent 
lines (i.e., the Q branches of the component bands 
of the perpendicular band) will arise from tran- 
sitions of the type AJ=0, AK=+1 we ob- 
tain as the spacing between any two lines 
dv=h/[4n?(A +B) ]. Taking as the value for Av 
the average spacing in the two types of ob- 
served bands and letting A=0.9B we obtain 
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A=2.7X10-" g cm’, B=3.110-" g cm?® and 
C=5.8X10~* g cm’, in good agreement with the 
values [4 = 2.667 X10-*° g cm?, Ip = 3.076 XK 10-*° 
g cm’, and J¢=5.845X10-* g cm? given by 
Cross. These values for the moments of inertia, 
when the axis of the smallest moment of inertia is 
taken along the bisector of the apex angle, lead to 
a value for the apex angle of about 85°.*: 7 

It is interesting to compare the fundamental 
bands here observed with the fundamentals 
observed in the spectrum of water vapor. In the 
water vapor molecule where the vertex angle 
appears to be about 104° the axis of smallest 
moment of inertia is normal to the bisector of the 
apex angle. The low frequency of vibration v2 
along the bisector of the vertex angle has here 
been identified as the observed frequency 1595 
cm~ and its position in the spectrum as well as 
its rotational character is consistent with such 
an identification. The frequency v3, a high fre- 
quency vibration normal to the bisector of the 
apex angle, has been taken as the observed 
frequency 3756 cm™, and also its position in the 
spectrum and its rotational character are in 
agreement with such an identification. This 
leaves the Raman frequency 3604 cm™ to be 
taken as v1. Hence the observed fundamentals 
found in the spectrum of hydrogen sulfide are 
entirely analogous to those found in the spec- 
trum of water vapor. 

6 Further evidence in support of such a model will be 
produced by A. H. Nielsen and H. H. Nielsen in an article 
soon to appear dealing with the spectrum of D.S and HDS. 

7 The value of the vertex angle of 110° given by Bailey, 
Thompson and Hale as determined from their work seems 
to us as hardly compatible with the evidence gained from a 
study of the rotational structure of these bands. Quite 
aside from the question of whether the molecule is an ob- 
tuse or an acute angle isosceles triangle, the ratio of the 
moments of inertia assumes an impossible value for such 
a model. Computing the value of the parameter p defined 
above from Bailey, Thompson and Hale’s 110° model one 
obtains a value p=0.45 and that this is a value much too 
small is readily seen by an examination of the diagrams 
prepared by Dennison. These show that for one of the 
band types, for a value of p=0.5 the Q branches are al- 
ready quite sharp and by far the most intense lines in the 
bands while in the other type of band a definite absence 
of lines at the center of the band is to be expected. This 
can hardly be said to be what is observed experimentally. 
Moreover, the regularity in the rotational structure ob- 
served in the H.S bands at low dispersion and especially 
the fact that the spacings between lines is almost the same 


in both types of bands seems incoherent with a value of 
p much less than 0.9, 
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The Analysis of the Raman Spectra of Si,Cl,(2) and of Si.H,(g)* 


FRED STITT AND Don M. Yost 
Caitfornia Institute of Technology, Pasadena, California 


(Received October 31, 1936) 


The Raman spectra of SizH¢(g) and of SisCis(/) have been photographed. Polarization meas- 
urements of the stronger lines of SisClz have been made. From a consideration of the relations 
among the modes of vibration of similar molecules having different symmetry, the analysis of 
the observed lines is facilitated by a knowledge of the assignments of the fundamental fre- 
quencies of SiH,, HSiCl; and BrSiCl;. A definite assignment is made for the three observed lines 
of Si,Hs. A tentative assignment of the observed lines of SizCls is made and some evidence is 
presented for the existence of internal rotation in this molecule at room temperature. 





N addition to determining and assigning some 
of the fundamental vibration frequencies of 
disilane and of disilicon hexachloride, we have 
attempted in this investigation to determine 
whether or not internal rotation is present in the 
SisClg molecule at room temperature. 


EXPERIMENTAL 


The disilane used was obtained from Professor 
Johnson of the University of Chicago, and we 
wish to express our gratitude to him for his kind- 
ness in making the material available to us. 
It was quite pure, showing, according to Profes- 
sor Johnson, a molecular weight of 62.6; the 
formula weight is 62.2. 

The disilicon hexachloride was prepared by 
treating ferrosilicon' with chlorine and fraction- 
ally distilling the product four or five times. It 
boiled at 145°+0.5°. 

The spectrum of gaseous disilane was photo- 
graphed with radiation from a mercury arc. A 
Pyrex Raman tube was used to avoid possible 
decomposition since Hg 2537A was found to 
decompose SiH,.? When the Si.Cls was photo- 
graphed in a Pyrex tube without a filter it 
showed some decomposition with the appearance 
of small brown flakes of silicon. An acid quinine 
chloride solution as filter absorbed almost com- 
pletely all lines of the mercury arc of higher fre- 
quency than the Hg 4358 group, and no decom- 
position occurred. As Raman lines excited by the 


* Contribution from the Gates and Crellin Laboratories 
of Chemistry of the California Institute of Technology, 
No. 574. , 

? Quig and Wilkinson, J. Am. Chem. Soc. 48, 902 (1926); 
Martin, J. Chem. Soc. 105, 2836 (1914). 

* Stitt and Yost, J. Chem. Phys. 4, 82 (1936). 
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Fic. 1. Experimental arrangement for polarization 
measurements. M, half-wave mica plate; N, Nicol prism; 
S, step weakener on slit of spectrograph; 7, Raman tube; 
L, condensing lens; B, parallel black screens. 


4047 group were observed in the absence of the 
filter, the decomposition was probably due to the 
Hg 3650 group. 

Polarization measurements were made on the 
stronger lines of SisCl, (1). The experimental ar- 
rangement used for these measurements is shown 
in Fig. 1. This arrangement suffers the disad- 
vantage that a separate exposure must be made 
for the parallel and perpendicular components,’ 
but possesses the two following important ad- 
vantages: (1) The use of a step weakener on the 
slit‘ automatically calibrates the photographic 
plate, the blackening along the length of each 
spectral line varying in steps of known relative 
intensity of illumination. Thus both the source 
and the exposure time for calibration are the same 
as for the image, an important factor where low 
intensities and long exposure times are involved. 


(2) The use of a half-wave mica plate® preceding 


’ An arrangement eliminating this disadvantage is ob- 
tained if, in Fig. 1, Mand N be removed and a large aper- 
ture Wollaston prism be introduced into the collimator tube 
of the spectrograph. Correction for polarization due to the 
spectrograph may then be eliminated by rotating the 
Wollaston 45° about the axis of the tube from the position 
in which the two beams emerging from it are in a vertical 
plane. The mercury arc is placed so that with the Raman 
tube it determines a plane at 45° to the vertical. This ar- 
rangement was not used because a suitable Wollaston 
prism was not available. 

‘ We are indebted to Professor Houston and Dr. C. M. 
Lewis for allowing us to use a step weakener made by them. 
See Phys. Rev. 44, 903 (1933). : 

5 We are indebted to Dr. John Strong for suggesting the 
use of a half-wave plate, 
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RAMAN SPECTRA OF Si2Cle 


the Nicol prism eliminates the correction for 
polarization due to the apparatus, and also al- 
lows each component to be photographed passing 
through the spectrograph with the electric vector 
in the horizontal plane, thus minimizing the loss 
due to reflection by the prism faces. 

For photographing one component (the ||- 
component for illumination in a vertical plane) 
the axis of the Nicol prism is horizontal and the 
axis of the mica plate parallel to that of the Nicol. 
The only alteration necessary for photographing 
the other component is to rotate the half-wave 
plate through 45°, thus rotating by 90° the 
electric vector of the light passing through it. 
Tests showed that a 0.048 mm mica sheet serves 
as a half-wave plate with practically no error 
from 4400A to 5000A. Only for high frequency 
shifts from Hg 4358A as an exciting line need a 
different thickness plate be used. 

The step weakener was calibrated at a number 
of different wave-lengths. For the calibration a 
low voltage Mazda lamp, such as used in street 
lights, served as constant source of illumination, 
and the small, spiral, vertical filament served as 
a good point source. The inverse square law of 
illumination was assumed throughout. Five posi- 
tions giving relative intensities of 1, 0.8, 0.6, 0.4, 
0.2 were used, each with exactly one hour expo- 
sure time. By microphotometering these con- 
tinuous spectra at any wave-length the calibra- 
tion for that wave-length could be made. Eight 
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TABLE I. Raman spectrum of Si.Hs (g). Pressure approx- 
imately 2 atmos. 








No. OF EXCITING 





Ay IN cm7! LINES DESCRIPTION 
2163 6 Strong; sharp. 
910 to 955 2 Weak; diffuse. 
434.5 3 Medium intensity; sharp; appeared 


also as an anti-Stokes line. 








steps were used having relative transmission 
values of 1.00, 0.85, 0.61, 0.45, 0.34, 0.26, 0.19, 
0.14 at 4384A. 

A Raman tube of rectangular cross section was 
employed in order to minimize reflection of inci- 
dent light. Parallel black screens for making the 
incident light approximately parallel were so 
spaced that only light within about 15° of the 
vertical could enter the Raman tube. For p=0.86 
the error in p is less than 0.02 for incident light at 
an angle of 20° with the vertical.* Eastman I —O 
plates were used. 

The mercury arc was run at constant voltage 
and showed little fluctuation in current over long 
periods of time. A 48-hour exposure of each com- 
ponent was taken, the illumination probably 
being very nearly the same for both components. 
However, it should be pointed out that, for this 
work, it is not necessary that the illumination be 
precisely the same for the two components for 
accurate results. For practically all compounds 

6G, Placzek, ‘‘Rayleigh-Streuung und Raman-Effekt,”’ 


Handbuch der Radiologie, Band I1, Teil 2 (Leipzig, 1934), 
p. 244. 


TABLE II. Raman spectrum of SizCle (1). (Anti-Stokes lines not included.) 











> - = RELATIVE POLARIZATION 
vy IN CM7! Av IN cm™! INTENSITY RATIO (p) DESCRIPTION 

22806 124 8 } 0.72+-0.02 Double line; diffuse. Not resolved in polarization measure- 

22814 132 5 ‘ : urements. 132 cm portion seems slightly less diffuse than 
124 cm™ portion. 

22759 179 0.4 Sharp. 

22726 212 4 0.86+0.02 Sharp. 

22685 354 <<i oes Hg 4047 as exciting line. 

22641 354 es | Hg 4339 as exciting line. 

22584 354 10 <0.05 Sharp. 

22561 377 } 1 Double; similar to 124, 132 cm doublet. 388 cm portion 

22551 387 seems slightly more intense and less diffuse than 377 cm™ 
portion. 

22517 421 <<<i Broad; too broad to be Av=425 cm of SiCly. 

22478 460 << <i Very weak; moderately diffuse. 

22451 487 <“—<s Very broad. Strongest of < < <1 lines. Too intense to be 
Av=588 cm from Hg 4339. 

22405 590 es | Hg 4347 as exciting line. 

22373 565 <<i Moderately diffuse. Too intense for Av =625 cm~ from Hg 
4347 as exciting line. 

22348 590 3 0.88 +0.03 Rather diffuse. 

22313 625 2 <0.5 Moderately diffuse. 

22232 706 <<<l Very weak; moderately diffuse. 
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TABLE III. Relations among the irreducible representations 
of certain symmetry groups. 
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some Raman lines have p=6/7 for ordinary 
incident light. If the exposure is somewhat differ- 
ent for the two components, one may find an ac- 
curate correction factor to use with the observed 
p values from the fact that the highest p is 
probably actually 6/7. In the present work, the 
observed values for two lines were 0.86 within 
narrow limits thus confirming our assumption 
that the illumination was the same for the two 
exposures. 

On plotting the microphotometer reading for 
each step as ordinate against the percentage 
transmission for that step as abscissa for each 
component, the depolarization p is the ratio of 
the abscissas for the same ordinate on the two 
curves. (Correction for background was unneces- 
sary as the background was negligible. Note the 
background radiation also is reduced in known 
ratios.) The values of p thus found choosiny 
different ordinates agreed quite well. 

The experimental results are summarized in 


Tables I and II. 


THEORETICAL 


Wigner’ has shown that the normal coordinates 
for vibrations of a system of point particles belong 
to (that is, have the linear transformation proper- 
ties of) definite irreducible representations (IR) 
of the point group representing the symmetry of 
the system. The methods for determining, for a 
given model, how many vibrational frequencies 
belong to each of the IR of the appropriate group 
and for determining the selection rules in the 
Raman and infrared spectrum are now well 
known.® 


7 Gottingen Nachrichten (1930), p. 133. 
8 See reference 6, 7; also E. B. Wilson, Jr., Phys. Rev. 


45, 706 (1934), J. Chem. Phys. 2, 432 (1934) and references 
given there. 
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It is instructive to follow (a) how the funda- 
mental frequencies for a model with a given 
number of atoms are distributed among the IR 
of the appropriate point group as the symmetry 
of the model is changed from that of one point 
group to that of another, and (b) and (c) how 
the new modes of vibration introduced when the 
number of atoms in a model is increased, with 
and without altering the point group symmetry, 
distribute themselves among the IR of the ap- 
propriate point group. The use of such considera- 
tions for correlation and assignment of frequen- 
cies is frequently very helpful. We will restrict 
ourselves here to a discussion of examples related 
to the interpretation of the spectra of SisHs 
and SisCle. 

For studying the distributions referred to in 
the preceding paragraph the methods already 
referred to® are employed, supplemented by some 
further considerations. When the symmetry is 
lowered to that of a subgroup of the original 
group, those IR of the original group possessing 
the same characters with respect to the operations 
of the subgroup merge to form a single IR of the 
subgroup. Mulliken® has published tables show- 
ing some of these relations. When the symmetry 
is altered so that the new point group is not a 
subgroup of the original one, a correlation be- 
tween the two is then obtained by seeing how the 
IR of each reduce to those of the group possessing 
all elements of symmetry common to the two. For 
our present purposes we are interested in these 
relations for the point groups C;,, D3, Dsa, Dsi, 
and T,; these are shown in Table III. IR desig- 
nated by A are one dimensional, and frequencies 
associated with them are nondegenerate; those 

TABLE IV. Irreducible representations to which the funda- 


mental frequencies of certain models belong, and selection 
rules for Raman and infrared spectra. 























For- Sym- INFRARED 
MULA | METRY FREQUENCIES RAMAN ACTIVE ACTIVE 
XY3 | C3, | 2Ai, 2E 2A,, 2E 2A;, 2E 
XY, | Ta Ai, E, 2T2 Ai, E, 2T2 2T» 
XY3Z} Cz, | 3Ai, 3E 3Ai, 3E 3Ai, 3E 
Xe2Ye| Dea | 3Aig, Aiu, 2Aeu | 3Aig, 3Eg 2Aou, 3Eu 
3E,, 3E, 
X2Ye| Ds, | 3Ai’, Ar’, 2Ae”’| 3A)’, 3E’, 3E” | 2A2”, 3E 
3E’, 3E” 
X2Ye} Ds; | 4Ai1, 2A, 6E 4A,, 6E 2Ao, 6E 











Note: By active is meant only that a frequency is permitted to 
appear, and not that it necessarily does appear. 


® Phys. Rev. 43, 294 (1933). 
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designated with letters E and T are two and three 
dimensional, respectively, and the frequencies 
associated with them are doubly and triply 
degenerate. 

In Table IV are shown the distribution of the 
fundamental vibrations for several models among 
the IR of the appropriate point group, and also 
the selection rules for the fundamentals in their 
infrared and Raman spectra. For those vibrations 
belonging to the completely symmetric IR (A; of 
D;, A; of Czy and of Ta, Ay’ of Ds;, and Ai, of D3.) 
the corresponding Raman lines have a polariza- 
tion ratio p< 6/7; for all others p=6/7. 

In Table V are shown the Raman selection 
rules for first overtones and one-one combina- 
tions for these same models. 

Coordinates belonging to appropriate IR for a 
given model of X2Y. may be constructed readily 
from similar coordinates for the XY3 model. 
(The normal coordinate for a given vibration is a 
linear combination of those symmetry coordin- 
ates belonging to the IR with which that vibra- 
tion is associated.) The method is as follows: For 
each symmetry coordinate of XY; construct two 
for X2Y¢. In both of these each —XY3 radical 
vibrates with the XY; symmetry coordinate, but 
in one the two radicals vibrate in phase, in the 
other 180° out of phase. This will be made clear 
by reference to Fig, 2, in which are shown the 
symmetry coordinates thus'® found from the 
symmetry coordinates formulated by Howard 
and Wilson" for X Y3. The “‘in phase” coordinates 
(those with subscript s) are unchanged by rota- 
tion about a twofold axis of the X2Y¢ model, 


TABLE V. Raman selection rules for first overtones and for 
one-one combinations of fundamental frequencies of X2Ye 
models. All are allowed except those shown. 











SYMMETRY FORBIDDEN RAMAN LINES 
D; Ai+Az2 
D3a (Aig or Aiu)+(Azg or Agu), Aig+Eu, Aiu+E,, 
Ao, +E., Ao +Ey,. 
Ds, (Ay’ or Ai’’)+(Ad)’ or A2’’) 











Slight modification of the XY; coordinate may be 
necessary to eliminate over-all rotation as in forming 7, of 
Fig. 2 from wz, w, of XY3, or to avoid distortion as in 
forming wz, wy of X2Y¢ in Fig. 2 from T,, T, of XY3. We 
a to thank Dr. J. B. Howard for calling our attention 
to this. 

"J. B. Howard and E. B. Wilson, Jr., J. Chem. Phys. 3, 
630 (1934). 


RAMAN SPECTRA OF Si2Cle¢ 
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Fic. 2. Symmetry coordinates for X2Y¢ from those for XY3. 


whereas the same operation changes the “‘out of 
phase”’ coordinates (those with subscript a) by 
180° in phase. The equilibrium configuration for 
X2Y¢ shown in Fig. 2 is that for symmetry D3a, 
but beneath each coordinate is noted the IR of 
the point groups D3, Ds, and D3, respectively, 
to which the suitably drawn coordinate belongs 
when the —XY; radicals are combined to form 
molecules of these symmetries. These coordinates 
have been given designations for future reference. 

In the following discussion we shall refer to 
these coordinates as if they were normal coordin- 
ates. The symmetry coordinates of AB; have been 
chosen so that each one involves primarily a 
bending or a stretching of the bonds. For most 
AB; molecules the large difference in the frequen- 
cies of the two vibrations belonging to either A; 
or E of C3, show that the normal coordinates are 
probably very similar to these symmetry co- 
ordinates. The same is probably true of the co- 
ordinates shown for ABs. This would be ex- 
pected to be true, especially for the degenerate 
modes, in case of free internal rotation. 
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DISCUSSION 
Si.H, 


There is no question as to the presence of in- 
ternal rotation in SigH¢. In ethane the difference 
in stability of the Ds, and D2 configurations is 
considerably less than kT at room temperature 
so that the internal rotation is almost free.!? In 
SisHs, due to the greater distance between the 
two central atoms than in C2H¢, rotation of the 
two halves of the molecule with respect to each 
other must be even less hindered than in ethane. 

Av=2163 cm™ is doubtless due to »,(z) 
because of its great intensity, its sharpness, and 
its numerical proximity to v(7)=2187 cm™ of 
SiH? 

Av=910 to 955 cm™ is probably 6(c). In SiH, 
we have assigned? v=910 cm™ to IR Tz of Ta, 
v=978 cm to E of Ta, and of these two frequen- 
cies only Avy=978 cm™ was observed with ap- 
preciable intensity in its Raman spectrum. On 
the assumption that the intensities of the Raman 
lines of SioHs are related to those of SiH, in 
accordance with Table III we conclude Av=910 
to 955 cm! of SigH¢ probably belongs to E of Ds, 
hence is probably 6(c) rather than é(7). 

Av=434.5 cm™ is very likely 7. The funda- 
mental frequencies of SiH, indicate that all the 
fundamental frequencies of SiegHe are in the 
neighborhood of either 1000 cm or 2000 cm™ 
except 7, w; and r. w; is probably zero as nearly 
free internal rotation is present. » would be ex- 
pected to be sharp, and perhaps of appreciable 
intensity because of its appearance in the Raman 
spectrum of C2H, (g).!* If it is assumed that 434.5 
cm~! is the frequency of vibration of two rigid 
—SiH; groups with respect to each other, a 
stretching force constant of 1.710° dynes/cm 
is found for the Si—Si bond. 


Si.Cl, 


In Si.Cl, the situation as to internal rotation is 
not so clear. The electrostatic interaction due to 
assuming a positive charge of 4.5X10~® e.s.u. on 
each Si atom and a negative charge of 1.5107 
e.s.u. on each Cl atom (corresponding to the high 
value of 3.0 Debyes for the dipole moment of the 

12 Eucken and Weigert, Zeits. f. physik. Chemie B23, 
266 (1933); Eucken and Parts, Géttingen Nachrichten 
(1932), p. 274; Zeits. f. physik. Chemie B18, 61 (1932); 


Eyring, J. Am. Chem. Soc. 54, 3191 (1932). 
13 Teller and Topley, J. Chem. Soc. 885 (1935). 
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Si—Cl bond)" results in stabilizing the D3a con- 
figuration by 475 cal./mole with respect to the 
D;, configuration. The van der Waals attraction 
between the Cl atoms tends to stabilize the D3, 
configuration by an amount of the order of 450 
cal./mole.’* Due to the mutual repulsion of the 
Cl atoms, there is a barrier of the order at 500 
cal./mole between the two configurations." 
Hence it seems probable that resistance to in- 
ternal rotation at room temperature is probably 
of the order of kT, the Dsa configuration being 
somewhat more stable than the Dg, one. 

In accordance with Table IV we see that if 
SieCl¢ possesses a rigid staggered (D3,) structure, 
six distinct fundamental frequencies are _per- 
mitted in the Raman spectrum, three of which 
would havea polarization ratio p< (6/7). If SieCl, 
possesses a rigid structure having a symmetry 
plane perpendicular to the threefold axis, nine 
distinct fundamental frequencies are permitted 
to be Raman active, and three of these would 
have p<(6/7)." If free internal rotation occurs 
the symmetry of the molecule is only D; most of 
the time. For a rigid D; model Table IV shows 
that ten distinct fundamental frequencies are 
permitted to be Raman active. However, if free 
rotation is present, one of these, w;, is zero, and 
intuitively it seems very likely that under such 
conditions the corresponding ‘‘in phase” and 
“out of phase’”’ vibrations not symmetric with 
respect to the threefold axis (the degenerate vi- 
brations) would become indistinguishable, thus 
resulting in only six distinct fundamental fre- 
quencies permitted in the Raman spectrum. 

4 Distances used: Si— Cl =2.00A, the SiCl, value. Si—Si 
= 2.34A the value recently found in SiH, by electron dif- 
fraction. (Dr. Brockway, private communication.) The 
C—Cl dipole moment is about 1.5 Debyes and the Si— Cl 
separation is about twice the C— Cl separation on Pauling’s 
— scale (See J. Am. Chem. Soc. 54, 3570 

1 : 

15 See Slater and Kirkwood, Phys. Rev. 37, 682 (1931). 

16 If for the repulsive potential between two bound Cl 
atoms, the expression (7) = be@¥0-” /? is used, with b= 10 ” 
erg, p=0.2091A, and yo=1.370A, the calculated barrier !s 
about 350 cal. /mole. These values of 5b and p were used to 
give the best fit for neon (see Bleick and Mayer, J. Chem. 
Phys. 2, 252 (1934)) and yo was chosen by comparison with 
vo for neon and for the alkali and halide ions (see Huggins 
and Mayer, J. Chem. Phys. 1, 643 (1933)). The barrier 
(exclusive of Coulomb interaction) calculated if the Cl. 
atoms are assumed to be chloride ions is about 2000 cal. / 
mole. (See Huggings and Mayer, above.) The correct value 
is probably nearer that calculated for bound chlorine atoms. 

17 By rigid here is meant that resistance to internal rota- 


tion is sufficiently great so that only small vibrations about 
the equilibrium position may occur. 
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In addition to the selection rules, we have some 
knowledge of the fundamental frequencies of 
HSiCl;'8 and of BrSiCl,'* to aid us in our analysis 
of the SisCls spectrum. In Table VI are shown 
the observed Raman frequency shifts for these 
compounds and the assignment of the frequen- 
cies. The assignment shown for the HSiCl, fre- 
quencies is that of Mecke.”° The analysis of the 
BrSiCl; frequencies we have made by analogy to 
Mecke’s analysis?® of the HSiCl;, HCX; and 
CH;X frequencies (X=halogen). The designa- 
tions v and 6 refer to vibrations roughly ap- 
proximated by the symmetry coordinates shown 
in Fig. 2 for XY; where X now becomes a rigid 
H —Si (or Br—Si) group, and (7) and n(c) refer, 
respectively, to vibrations in which the motion is 
essentially a stretching or a bending with respect 
to the —SiCl; group of the H—Si (or Br—Si) 
bond. Calculation of (7) for BrSiCl;, assuming 
it to involve only a vibration of the Br atom 
relative to a rigid —SiCl; group and assuming 
the stretching force constant of the Si— Br bond 
to be unchanged from its value in SiBr;, yields 
314 cm“. 

Considering both BrSiCl; and SizCle as a 
—SiCl; group with a rather large mass attached 
along the threefold axis, one would expect the 
v(r) and 6(7) vibrations in one to be of approxi- 
mately the same frequency as in the other. 
Without concluding as to the symmetry of SisCle, 
this prediction coupled with the selection rules 
and the values of p for these lines lead us to 
assign 354 cm~!= (7) and either 124 or 132 cm™ 
= 5(7) in Si2Cl.. 

Since the mass of the H atom is so small rela- 
tive to that of Si or Cl, the frequencies of HSiCl; 


TaBLE VI. Assignment of the fundamental vibration frequen- 
cies of HSiCl; and of BrSiCls. 











HS1Ci3 IR DESIGNATION BrS1Ci3 
489 (6.6) Ai v(m) 362 (Vv. STR.) 
249 (2.6) Ai 5(1) 123 
2256 (8.3) Ai n(m) 325 

587 (1) E v(o) 410 

179 (4) E 5(c) 183 (sTR.) 
797 (1) E n(a) 201 

















‘8 de Hemptinne and Peeters, Bull. Sci. Acad. Roy. Belg. 
i (1931); Urey and Bradley, Phys. Rev. 37, 843 
1). 
*de Hemptinne, Wouters and Fayt, Bull. Sci. Acad. 
Roy. Belg. 19, 318 (1933). 
°° Handbuch und Jahrbuch der 


Chem. Phys. (Leipzig, 
1934), Band 9, Teil II, p. 390. 
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(omitting n(7) and n(c))*! are very nearly those 
of a free —SiCl; group. If nearly free internal 
rotation is present in SieCl¢, the perpendicular vi- 
brations v(o) and 6(¢) would be expected to be of 
nearly the same frequency as those of a free 
—SiCls group, since then the net instantaneous 
momentum in any direction perpendicular to the 
threefold axis would be essentially zero for each 
half of the molecule. This condition need not 
obtain for a rigid” configuration. The appreciable 
intensity, the polarization ratio of 6/7, and the 
numerical value of 590 cm™ in Si,Cl, indicate 
that it is v(o). This proximity to the value of v(c) 
= 587 cm in HSiCl; favors the conclusion that 
there is nearly free internal rotation in SieCl,. 
The appreciable intensity and p value of 212 
cm lead to its assignment as 4(c). This is 
slightly higher than 6(¢) in either HSiCl; (179 
cm~') or BrSiCl; (183 cm). 

We are unable to make definite assignments of 
more than these four lines as fundamental fre- 
quencies of SisCls. The relation of these assign- 
ments to the analyses of the frequencies of 
HSiCl;, BrSiCl; and SiCl,; from the viewpoint of 
Table III is shown in Table VII. 

Several of the weaker lines. observed may be 
assigned as combinations or overtones: 


421 cm": The first overtone of 212 cm™. 
487 cm: Combination 132 cm!+354 cm™!. 
565 cm™!: Combination 212 cm=!+354 em. 
706 cm: First overtone of 354 cm™'. 


The remaining weak lines, Av=377, 387 and 
460 cm~', are also probably due to combinations 
or overtones. Numerically they agree with 124, 
132 and 212 cm“, respectively, combining with 
a frequency of about 252 cm~. As fundamentals 
we can see no feasible assignment for them. 

The component of the 124, 132 cm= doublet 
which is not 6(7) may well be 7 with a polarization 
ratio of 6/7. The observed p=0.72 is an over-all 
value for both lines. 

Two possibilities suggest themselves for 179 
cm: (1) It may be 7; (2) If the correct configu- 
ration is a rigid” one, it may be 6,(¢), its numer- 
ical value being practically the same as 6(c) in 
HSiCl; and BrSiCl;. If » is assumed to be the 
vibration of two rigid —SiCl; groups with respect 
to each other, and if the Si—Si bond stretching 


*1 This is equivalent to calling the symmetry coordinates 
described above for n(z) and (c) normal coordinates. 
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TABLE VII. Correlation of fundamental frequencies of Taste VIII. Assignment of the Raman frequencies of SizCl., 
SiCh, HSiCl;, BrSiCls and SieCls. assuming tnternal rotation.” 
S1Cus (Tg) | HS1Cus (Cs) BrSICLs3 SrCLe (Ds, Dsa)(?) Av cm! INTENSITY ASSIGNMENT (See Fig. 2) 
A 425 (10) 489 (7)) {362 (v. stR.) | 354 (10)) FUNDAMENTALS 
250 (3) tA 123 124 (8) ~* 
or isi 354 10 v,(7) 
ei en | satis 132 (or 124) 5 (or 8) 6,(7) 
T2 220 (5) 179 0.4 n (?) 
799 (1)) (201 a . - 
E 150 (5) 587 (1)}E | 4410 590 (3) E,E g) 
179 (4)} (183 (sTR.) 212 (4)f di 124 (or 132) 8 (or 5) r (?) 




















force constant is that calculated above from SizH., 
its frequency should be 201 cm~. (A similar 
calculation using C2H,'* and C.Cl,” data predicts 
a frequency of 354 cm™ for n of C2Cl¢, and a line 
belonging to the completely symmetric IR is ob- 
served at 342 cm™.) Since the frequency corre- 
sponding to the y vibration in CsH¢, SigH¢ and 
C.Cl¢ in each case appears with appreciable in- 
tensity in its Raman spectrum, we should expect 
it to appear also in SieCls. Therefore it seems to 
us more likely that 179 cm~ is n rather than 6,(¢). 

The observed line at Av=625 cm™ appears 
with considerable intensity, yet it does not seem 
possible to us that it is due to a fundamental fre- 
quency. Its polarization ratio shows that it is not 
due to a degenerate fundamental frequency. 
Reference to Table IV shows that only »,(7), 
5,(7), n, and w; are permitted to be Raman active 
among the nondegenerate fundamentals in any 
of the models considered. Of these only »,(7) 
would be expected to be of relatively high fre- 
quency, but 354 cm~ is definitely assigned as 
vs(7). Avy=625 cm@ must therefore be due to a 
combination or overtone. Its intensity then indi- 
cates that it probably corresponds to a level in 
accidental degeneracy with 590 cm™. Since 590 
cm- shows a p value of 6/7 and 625 cm a p 
value of <0.5, and since interaction between two 
levels may occur only if their wave functions 
belong to representations including the same IR, 
then the wave function corresponding to 625 cm 
must belong to a representation including both 
A, and E of the appropriate point group. Of the 
first overtones or one-one combinations for each 
of the symmetries D3, D3¢ and Ds, only an over- 
tone of a degenerate fundamental or a combina- 
tion of two degenerate fundamentals belonging 
to the same IR satisfy this requirement.”* If this 


22 Heidenreich, Zeits. f. Physik 97, 277 (1935). 
23 See Tisza, Zeits. f. Physik 82, 48 (1933). 








COMBINATIONS AND OVERTONES 


421  —_e 4 25(c) 

487 Cae eI 5,() + v5(7) 

565 < <i 5(0) + y5(7) 

706 <<<il 2v;(7) 

377 <<e“y 124cm1+252cm™ (?) 
387 < <1 132cm?+252cm™ (?) 
460 Cele 5(o)+252cm™ (?) 
625 2 vs(w7) +252cm-'or (?) 


in accidental deg. with 


v(o)  (?) 








interpretation is correct®* either 26(¢), 27 or 
5(¢)-+7 should be approximately 608 cm~, al- 
though this is not the case for the present pro- 
posed assignment. It is interesting to note that 
combination of 354 cm™ with the hypothetical 
252 cm™ frequency gives a level at approxi- 
mately 608 cm, but this level would not satisfy 
the requirements just mentioned. 

The above analysis does not permit a definite 
answer to the question of internal rotation in 
Si,Cls at room temperature, but affords some 
evidence for its existence. In Table VIII is shown 
a tentative assignment of the observed lines 
under the assumption that free rotation does 
exist.” 

In conclusion we wish to thank Professor E. B. 
Wilson, Jr. and Dr. J. B. Howard for valuable 
criticism and discussion. 


Note added in proof. Very recently Kemp and Pitzer 
(J. Chem. Phys. 4, 749 (1936)) have presented evidence for 
the existence of a potential barrier of about 3150 cal. /mole 
between the D;4 and D3, configurations in ethane. Since 
the Si—Si distance (2.34A) in SisHg and Si2Cle is con- 
siderably greater than the C—C distance (1.52A) in ethane, 
their conclusions would not appear to exclude the possibility 
of internal rotation in SigH¢ and SicClg. 


*4In this case the numerical agreement between 587 
cm of HSiCl; and 590 cm= of SizCls loses some of its 
significance as evidence for free rotation. 

2 If instead a rigid (reference 17) Dsu, or even Daa, 
structure exists, subscripts s on the degenerate frequency 
designations in Table VIII would be the only change neces- 
sary for the appropriate new tentative assignments. 
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Photochemical Reduction with X-Rays and Effects of Addition Agents 
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Ceric sulfate solution was reduced by exposure to 
x-radiation. In order to follow this reaction a very accurate 
method was developed for the determination of ceric 
sulfate in dilute solutions by the use of orthophenanthroline 
indicator. The amount of reduction of ceric sulfate by 
x-rays was found to be practically independent of concen- 
tration. Reduction was shown to be directly proportional 
to time of exposure up to 50 percent and slightly less 
efficient beyond that point. Silver perchlorate and silver 
nitrate decreased the amount of reduction of ceric sulfate 
by x-rays. Mercuric nitrate and acetic acid increased the 
amount of reduction caused by x-rays. Curves are also 
presented to show the effect of changing the concentration 
of these addition agents. Indications are that the mercuric 
nitrate was truly a catalyst to the reduction of ceric sulfate 
by x-rays and that the acetic acid increased the amount 
of reduction by direct action on the ceric sulfate under the 
influence of x-rays. Perchloric and nitric acids increased 
the reduction of ceric sulfate appreciably but not to as 
great an extent as did Hg(NOs)o. Sodium nitrate, per- 
chlorate and sulfate showed no effect. Various organic 
substances were tried and found to increase the reduction 
by varying amounts. Experiments were carried out with 


various combinations of the catalysts and inhibitors to 
determine their relative efficiencies. The reduction of 
potassium permanganate by x-rays proceeded to manga- 
nese dioxide for low acid concentrations and to manganous 
ion at higher acidity. The reaction was greater than for 
ceric sulfate but not as nearly proportional to dosage. The 
effect of addition agents was similar to that for ceric 
sulfate. X-rays did not produce the same reaction in 
solutions of uranyl sulfate-oxalic acid as does ultraviolet 
radiation. Solutions of potassium iodate were reduced by 
x-rays but to a less degree than ceric sulfate or potassium 
permanganate. Hydrogen peroxide was formed in samples 
during radiation but did not effect the potassium iodate. 
Solutions of potassium bromate were not reduced by 
x-rays but even larger amounts of hydrogen peroxide were 
formed. Calculations were made for: ions reduced per r 
unit absorbed, r units necessary for reduction of one 
equivalent, ions reduced per ion pair formed by x-rays, 
and calories of absorbed energy necessary for reduction of 
one equivalent, by two independent methods. A mechanism 
was proposed for the formation of hydrogen peroxide by 
x-rays in the presence of dissolved oxygen. 





HOTOCHEMICAL reduction of inorganic 
compounds by x-rays is an interesting field 
of investigation which already has received some 
attention. Chamberlain! has reported the reduc- 
tion of potassium permanganate and potassium 
iodate by x-rays. These results were only qualita- 
tive. Clark and Pickett? made the same observa- 
tion and did quantitative work on the photolysis 
of potassium nitrate and potassium iodide by 
x-rays. Fricke and Brownscombe,* and Fricke 
and Washburn‘ have noted the reduction of 
potassium dichromate by x-rays and attributed 
it to the formation of hydrogen peroxide. These 
workers have also observed large effects of im- 
purities in the irradiated solution. 

The formation and decomposition of hydrogen 
peroxide by x-rays have been the subject of much 
investigation because of the possibility of its 
influence on any x-ray photochemical reaction in 
aqueous medium. Glocker and Risse®: ® studied 
the effect of wave-length and concentration upon 
the decomposition of hydrogen peroxide. Risse 
reported the formation of hydrogen peroxide by 
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x-rays in water containing dissolved oxygen. 
Fricke’? and Brownscombe® confirmed his report 
and showed that no hydrogen peroxide was 
formed in pure air-free water. 

The purpose of the present work was to con- 
tribute to the knowledge of the quantitative 
aspects of reductions by x-rays of familiar in- 
organic reagents in dilute aqueous solution and 
the catalytic or inhibitory effects of various 
substances added voluntarily to these solutions. 


EXPERIMENTAL 


The x-ray tube used for the irradiation of solu- 
tions during this investigation was one designed 
by the authors and constructed in the machine 
shop at the University of Illinois. It was designed 
with a metal body for self-protection from stray 
radiation, a water-cooled copper or tungsten 
anode to allow a high tube current, and a special 
aluminum-beryllium alloy window for the eff- 
cient transmission of x-rays but not light or heat. 
A complete description of this tube will appear 
elsewhere. 








Determination of intensity of radiation 


Intensity measurements were taken with a 
Victoreen R meter and corrections made for 
temperature and atmospheric pressure to give 
values in terms of the internationally accepted 
_ r unit® of x-radiation. The effective wave-length 
of x-rays emitted from both copper and tungsten 
targets operating at 55 kilovolts was determined 
by the method of Duane, Hudson and Sterling’? 
involving the determination of percent trans- 
mission by a 2 mm aluminum plate. Values of 
0.62A for the copper target and 0.60A for the 
tungsten target were obtained. The general ex- 
perimental procedure for irradiation was to ex- 
pose 10 cc of the solution in a 50 cc beaker at a 
distance of 15 or 20 cm from the focal spot, with 
the tube operating at 55 kilovolts. Some ex- 
posures were made with a tube current of 20 and 
others with 30 milliamperes. 

Absorption coefficients calcvlated for water 
for the two effective wave-lengths used (u=0.826 
for \=0.62A and u=0.746 for \=0.60A) agree 
satisfactorily with interpolations from a series 
given by Glocker and Risse.® 


Method for estimation of ceric sulfate 


Since it was desired to determine quantitatively 
the amount of ceric sulfate reduced by reasonably 
short exposures to x-rays a very precise method 
of titration had to be developed. Potentiometric 
titration with 0.001 N ferrous sulfate was tried 
but found to be only partially satisfactory. Since 
orthophenanthroline" has been widely used as an 
oxidation-reduction indicator in solutions of or- 
dinary concentration and found to give a very 
sharp end point, it was thought probable that 
application might be made of this indicator even 
in dilute solutions. The following procedure has 
been used and found very satisfactory: the 
orthophenanthroline ferrous complex was made 
up to one-fifth the usual concentration (0.005 
instead of 0.025 M) and one drop of the solution 
used for a 10 cc sample. Sufficient sulfuric acid 
was added to the sample to make it approxi- 
mately 1 NV and titration carried out with 0.001 NV 
ferrous sulfate. Under these conditions the end 
point was detected as the first appearance of a 
permanent pink color in contrast to the colorless 
solution before the end point. This change was 
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noticeable with 0.01 cc of 0.001 N ferrous sulfate 
and very pronounced with 0.02 cc, particularly 
when illuminated with a strong light on a white 
table and with a white background. It is evident 
that the indicator correction was very large since 
the indicator contained 0.005 NV Fett. However, 
nonirradiated controls were always titrated at 
the same time as the irradiated samples so the 
difference due to x-radiation did not depend on 
indicator correction as long as it was kept con- 
stant. The proper correction was made when 
direct comparison of the two solutions was 
necessary. 

The procedure finally adopted for the prepara- 
tion of these ceric sulfate solutions was to re- 
crystallize the ceric ammonium sulfate until a 
very pure product was assured; then it was 
ignited in a muffle furnace at 900 to 1000°C for 
eighteen hours to convert it to the oxide. De- 
composition took place at a lower temperature 
but also some reduction was caused by the am- 
monia liberated. The prolonged heating in the 
presence of oxygen of the air reoxidized the 
product to ceric oxide. This oxide was then 
treated with a small amount of concentrated 
sulfuric acid (9 parts HeSO, to 1 part H2O) to 
convert it to the sulfate. The purified ceric sulfate 
was finally dissolved in water and diluted to con- 
venient stock concentration. The solutions were 
always strongly acidic with the excess sulfuric 
acid necessary to react with the oxide. 

These solutions, as well as all others used for 
irradiation in this investigation, were prepared 
from “conductivity”? water redistilled from 
alkaline permanganate. 


Reduction of ceric sulfate solutions by x- 
radiation 


Experiments showed that the amount of reduc- 
tion was practically independent of the concen- 
tration of ceric sulfate for dilute solutions. This 
is consistent with most photochemical effects of 
x-rays. The amount of reduction did increase 
slightly for more concentrated solutions of ceric 
sulfate (above 0.001 /). 

Since the ceric sulfate stock solution was quite 
acidic with sulfuric acid the dependence of reduc- 
tion upon acid concentration was determined. 
Increasing the sulfuric acid concentration from 
0.05 N to 1.2 N decreased the amount of reduc- 
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tion by a given dosage of x-rays in 0.0002 M 
Ce(SO,)2 by approximately 10 percent. 

The next step was to determine the function of 
reduction with time of radiation. This series of 
experiments was carried out with 10 cc samples 
of 0.0002 M and 0.0004 M Ce(SO,)e in 0.1 N 
H.SO,. Fig. 1 is a curve showing the micro- 
equivalents reduced against minutes irradiation 
(20 milliamperes, 590 r/min. incident intensity). 
It can be seen that the amount of reduction is 
very nearly proportional to dosage up to 50 
percent. The sixty minute value for 0.0002 M@ 
Ce(SO,4)2 which corresponds to approximately 75 
percent reduction is somewhat low indicating 
that the rate of reaction has decreased slightly. 
The one hundred and twenty minute value for 
0.0004 M Ce(SOx4)2 shows a similar behavior. 


Effect of added substances upon reduction of 
ceric sulfate 


During the investigation of the effect of various 
dosages upon amount of reduction one sample 
was exposed after a drop of the orthophenan- 
throline indicator had been added. The amount 
of reduction in that sample was almost double 
that of samples given the same irradiation with- 
out the indicator present. Further experiments 
with Fe** added to the concentration present in 
the indicator failed to cause any irregularity in 
results. The increased reduction, then, must have 
been due to the organic part of the indicator. 

It was thought of interest to determine what 
effect small amounts of Ag+ might have on the 
amount of the reduction of ceric sulfate by 
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x-rays. Inasmuch as the Ag* would cause the 
solution to absorb a larger percentage of x-rays 
it might be predicted that a greater amount of 
reaction would take place. On the contrary Ag* 
was actually found to be a very efficient inhibitor 
for the reaction. Fig. 2 is a curve showing amount 
of reduction in 10 cc of 0.0002 M Ce(SO,)s, 
0.1 N H.SO, against the log of the concentrations 
of AgClO, present. AgNO; was also used and the 
values agree within experimental error with 
those for AgClO,. These values were taken with 
copper target radiation with the tube operating 
at 55 peak and 30 milliamperes for fifteen minute 
exposures. The solutions were placed 15 cm from 
the focal spot and the incident intensity of radia- 
tion was 775 r/min. From the curve it is seen that 
even 107° M AgClO, has a definite inhibitory 
effect upon the course of the reduction and 
0.1 M7 AgClO, decreases the amount of reduction 
to approximately one-third of its ordinary 
amount. 

Since silver salts were found to have such a 
marked effect, other salts of heavy metals were 
tried to see if they might behave in a similar way. 
However, mercuric nitrate was found to show a 
very strong behavior in the opposite direction, 
that is, increasing the amount of reduction of 
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TABLE I 








ADDITION AGENT 
(CONCN. APPROX. 0.02 M) 


MICROEQUIVALENTS 
CE(SO4)2 REDUCED 





None 0.48 
Ethyl alcohol 61 
Diethyl ether i 
Dibutyl ether .73 
Acetone 51 
Carbon tetrachloride 57 
Ethyl acetate 51 
Dioxane 55 
Benzene 1.12 








ceric sulfate. One curve of Fig. 3 shows the 
amount of reduction in 0.0002 M@ Ce(SO,)s, 
0.1 N H.SO, plotted against the logarithm of 
molar concentration of Hg(NO3)e for the same 
conditions of irradiation as used for sample repre- 
sented by Fig. 2. A curve is also shown for acetic 
acid which has been found to have an effect 
similar to that of mercuric nitrate. It is well to 
mention again that the values given for amount 
of reduction are differences in amount of standard 
FeSO, necessary to titrate the irradiated sample 
and a similar nonirradiated control. The values 
recorded are then due to the action of x-rays and 
independent of any direct action beween ceric 
sulfate and the added substance. 

The effects of HCIO, and HNO; on the reduc- 
tion of Ce(SO,4)2 by x-rays were also determined. 
In each case the amount of reduction was in- 
creased but not to as great an extent as for 
HC.H;02 and Hg(NOs3)2. Whereas, 0.48 micro- 
equivalents of 0.0002 M Ce(SOx)2, 0.1 N H2SO, 
were reduced under conditions stated above, 
0.1 M HCIO, increased this value to 0.57 micro- 
equivalents and 0.1 4 HNO; to 0.59 micro- 
equivalents. Bismuth nitrate was also found to 
increase the amount of reduction considerably 
but the separation of bismuthy] sulfate prevented 
an exact quantitative estimate under the same 
conditions as used for the other substances. On 
the other hand, 0.05 14 CuSO, acted more like 
AgClO, and decreased the amount of reduction 
to 0.38 microequivalents. NaNO;, NaClO,, and 
Naz2SO, were also tried but found to give no 
appreciable effect. 

Various organic substances were also tried. 
Some were found to increase the amount of 
reduction very markedly and others to a lesser 
extent. Table I shows the substances studied and 
gives the microequivalents reduced under the 
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same conditions of irradiation as _ preceding 
experiments. 

Since both catalytic and inhibitory effects were 
noticed for single substances added to ceric 
sulfate solutions, it was of interest to see if these 
behaviors would be additive when two substances 
were present in the same solution. Experiments 
with 0.06 Hg(NOs)z, 0.06 HC:H30s., and 
0.02 M Heg(NOs)o, 0.02 M HC.H;0.2 indicated 
that the reduction was increased to a greater 
extent than by either single substance but not as 
much as to be directly additive. The next com- 
bination tried was one catalyst and one inhibitor 
such as HC2H;O2 and AgClO,, or Hg(NOs3)2 and 
AgClO,. Curves in Fig. 4 show the amount of 
reduction for a constant concentration of the 
catalysts (0.02 M@ HC.H;02; 0.06 M Hg(NOs3).) 
in 0.0002 M Ce(SOx4)e, 0.1 MN HeSO, for varying 
concentrations of the inhibitor AgClO,. In both 
cases AgClO, was found to be much more efficient 
as an inhibitor than the other substance as a 
catalyst. This is easily seen on the graph since the 
points where the curves intersect the line drawn 
to indicate reduction in pure Ce(SO,)2 represent 
the condition of balance between catalyst and 
inhibitor. Thus, 0.0001 4 AgClO, (log=—4) 
neutralizes the effect of 0.02 MM HC2H;O2 and 
0.002 M AgClO, (log=—2.7) corresponds to 


Experiments with potassium permanganate and 
added substances 


The reduction of potassium permanganate has 
been noted by other investigators ; Chamberlain! 
reported reduction to Mnt*t* and Clark and 
Pickett? reported the formation of MnQOz. These 
results were only qualitative, probably because a 
satisfactory method for estimation of dilute 
solutions of potassium permanganate was not 
available. Preliminary qualitative investigation 
confirmed the observations cited and showed 
that MnO, was formed by irradiation of solutions 
of potassium permanganate less than 0.1 N with 
respect to H2SO,. Solutions of higher acidity 
(1.0 N HeSO,) were reduced directly to Mn** 
and varying amounts of MnO, and Mn** were 
formed for intermediate concentrations. 

Quantitative determination of reduction was 
found possible by a procedure similar to that used 
for ceric sulfate. However, the orthophenan- 
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throline end point with 0.001 N FeSO, was 
neither as sharp nor as permanent as for ceric 
sulfate and a titration error of approximately 
+0.03 cc was thus introduced. 

A study of the amount of reduction of KM nO, 
by different dosages of x-rays was carried out in 
a manner similar to that for Ce(SO,4)2. Ten cubic 
centimeter samples of 0.0004 VN KMnQ,, 1.0 NV 
H.SO, were used. Fig. 5 is a curve showing 
microequivalents reduction plotted against min- 
utes radiation, 20 cm from tungsten target, 
operating at 55 peak and 20 milliamperes (inci- 
dent intensity 590 r/min.). 

From the curve it is seen that the relation be- 
tween reduction of potassium permanganate and 
dosage is not so nearly linear as was the case for 
ceric sulfate. The curve shows a sharp rise for the 
first few minutes; this may indicate that x-rays 
catalyze the reduction of KMnQ, by slight 
traces of organic materials remaining in the 
solution in spite of precautions taken to prevent 
such contamination. This is further borne out by 
the fact that slight amounts of added organic 
substances were found to increase the reduction 
to a great extent. Orthophenanthroline as dilute 
as 0.00002 M increased the reduction for a ten 
minute exposure to three times its normal amount. 

Mercuric nitrate and acetic acid increased the 
reduction of potassium permanganate in the same 
manner as for ceric sulfate. Whereas, approxi- 
mately 0.49 microequivalents were reduced by a 
ten minute exposure of 0.0004 M@ KMnO,, 1 NV 
H.SO, solution, the presence of 0.05rM Hg(NOs)2 
or 0.05 M HC.H;02 augmented this value to 1.09 
and 1.12, respectively. 





(INCIDENT inTENSITY 590 P/oIN,) 
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The reduction of KMnQ, was inhibited by 
AgClO,,and AgNO; as would be expected by 
comparison with their behavior with Ce(SO,)>. 
This reduction was decreased from 0.49 micro- 


equivalents to 0.18 in the presence of 0.05 
AgClO, and to 0.22 for 0.05 M AgNOs. 


Experiments with uranyl sulfate-oxalic acid 
solution 


Urany] sulfate-oxalic acid solutions have been 
widely used as an actinometer for sunlight or 
ultraviolet radiation.” * The procedure usually 
employed is to titrate radiated samples and con- 
trols with standard potassium permanganate to a 
faint pink color, the difference in amount of 
permanganate being a measure of the amount of 
oxidation of H2C2O, by radiation, catalyzed by 
uranyl sulfate. Since the reaction is quite large 
with ultraviolet, it was thought that some reac- 
tion might be noted by x-rays. 

A solution of 0.005 M H.C.0,, 0.002 M@ 
UO.SO, was prepared and samples exposed to 
x-radiation. No reaction could be detected, so 
again it was found that behavior with ultraviolet 
was not safe evidence for prediction of the action 
of x-rays. Further experiments using a very large 
dosage of x-rays (17 cm, Cu target, 30 milli- 
amperes, 55 kv, 2 hours) even showed that ap- 
proximately 0.10 cc more 0.01 N KMnO, was 
necessary for the titration of the radiated samples 
than for the controls. However, the titration 
reaction was very slow near the end point for the 
radiated samples; this was probably due to the 
fact that a small amount of hydrogen peroxide 
had been formed by x-radiation and was re- 
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sponsible for the extra amount of _KMnQ, 
necessary. 


Experiments with potassium iodate and po- 
tassium bromate 


Since this investigation was primarily con- 
cerned with reduction effects of x-rays, it was 
desirable to examine the effect upon potassium 
iodate and potassium bromate. Stock solutions 
were carefully prepared from reagent quality 
chemicals and redistilled water. Solutions for 
exposure to x-rays were diluted to 0.0006 N and 
enough sulfuric acid added to make the solutions 
0.1 N. Ten cubic centimeter samples were placed 
15 cm from the copper target tube operated at 
55 peak and 30 milliamperes (incident intensity 
720 r/min.). 

The samples were examined by adding excess 
potassium iodide and titrating almost to the dis- 
appearance of the liberated iodine with 0.001 NV 
Na2S2O0;. Starch was then added and the titration 
continued to the disappearance of the blue color. 
This end point is capable of detection to +0.01 cc 
of 0.001 NV Na2S.O3. It was necessary to add more 
sulfuric acid for the titration of the bromate 
samples. 

Preliminary experiments with potassium io- 
date showed that reduction did take place but to 
a less extent than for ceric sulfate and potassium 
permanganate. The time of exposure was in- 
creased to one hour and several measurements 
taken on the amount of reduction for that time. 
An average of 0.55 microequivalents was re- 
duced. Titration of these samples seemed to 
proceed to a sharp and stable end point but upon 
standing thirty minutes or longer the radiated 
samples become blue again while the controls 
remained unchanged. Further titration with 
sodium thiosulfate solution after two hours gave 
another sharp end point which was stable over a 


TABLE II. 








FIRST SECOND 
TITRATION TITRATION 
cc 0.001 N cc 0.001 N 

NaA2S203 NaA2S203 


5.91 


SAMPLES 
10 cc, 0.0006 N KIOs, 
0.1 N H2SO« 





Controls 
Radiated one hour 


0.42 microequiv. H2O2 added 
(not radiated) 
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long period except as atmospheric oxidation 
caused a gradual change for the acidity used 
(0.1 N H2SO,). Repeated experiments revealed 
that the second titration of the radiated samples 
depended on time of radiation in the same manner 
as the reduction noted by the first titration. This 
second titration averaged 0.42 microequivalents 
for one hour exposures. 

Since hydrogen peroxide might have been 
formed in the radiated samples the next pro- 
cedure was to see if hydrogen peroxide added to 
the control samples would cause them to behave 
in the same manner as those that had been 
radiated. This was at first thought doubtful 
because the H,O2 might be oxidized by the KIO; 
or at least be destroyed during the first titration. 
However, a series of samples was prepared and 
the results are indicated in Table II. 

It is noted that H2Oz did not affect the KIO; 
and did not disturb the titration since the con- 
trols necessitated the same amount of Na2S.0; 
as did those containing H2Oe. There is a differ- 
ence of approximately 0.54 microequivalents be- 
tween the controls and radiated samples cor- 
responding to the reduction by x-rays. After the 
first titration the solutions were allowed to stand. 
After ten minutes the radiated samples and those 
to which H2O:2 had been added became slightly 
blue due to the slow liberation of Iz as the HO 
reacted with HI formed by previous addition of 
excess KI to the 0.1°N H2SO,. The blue color 
continued to develop in both solutions at the 
same rate proving beyond a doubt that H.O2 
had been formed in the radiated samples. The 
second titration was carried out after about two 
hours and corresponded to the amount of H.Oz 
formed or added. The amount of oxidation of the 
control samples by air was equivalent to 0.02 cc 
of NazS:O3. In more concentrated H2SO, solu- 
tions the reaction between HI and H2O2 was 
more rapid and consequently the end point for 
the first titration could not be detected as 
accurately. 

The effect of x-rays upon solutions of potas- 
sium bromate was tried in a manner similar to 
that for potassium iodate but mo reduction could 
be detected. However, the second titration 
showed that even larger amounts of H2O2 had 
been formed than in the case of potassium iodate. 
Approximately 1.01 microequivalents of H2O2 
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TABLE III. 




















Ce(SOx)2 KMnOx KIO8 H202 
Reduced Reduced Reduced Formed 
(0.1. N H2SO.) | (1.0N H2SOs) | (0.1.N H2SOs) | (0.1 N H2SO,) 
‘‘Equivalions” converted per 7 unit absorbed 4.4210” 7.45 X 10” 1.2710" 2.30 10" 
r units absorbed per equivalent converted 1.37104 0.81 10" 4.78 X10" 2.63 X10" 
Ions converted per ion pair 2.0 i 0.62 1.1 
Calories of absorbed energy per equivalent reduced* 379,500 205,135 1,224,193 690,000 
Calories electronic energy for conversion of one equivalent 344,000 204,000 1,095,000 602,000 




















Note: ‘‘Equivalions” equal gram equivalents times 6.06 < 1023, 


* Calculated from ions per ion pair, electron volts per ion pair, and heat equivalence of electron volts. 


were formed for a one hour exposure. Further 
experiments revealed that about the same amount 
of H:O2. was formed by same radiation of 
0.1 N H2SO, (not air-free). 

An attempt was made to study the effect of 
x-rays upon solutions of potassium chlorate but 
a satisfactory method for the accurate determi- 
nation of dilute solutions could not be found. 


Energy relationships 


Ceric ions reduced per r unit.—Values are given 
in Fig. 1 for amount of reduction of 0.0002 
and 0.0004 M ceric sulfate in 0.1 N HeSO, by 
tungsten radiation (effective \= 0.60) for varying 
periods of exposure. From the straight portion of 
the curve it can be seen that forty-one minutes 
were necessary for the reduction of 1.0 micro- 
equivalent weight or 6.0610" ceric ions. Since 
the incident intensity was 590 r/min. per sq.cm, 
590-11.94=7050 per min. incident radiation on 
the total surface. Of this amount 47.5 percent 
was absorbed by 0.1 N H2SO,; the absorption 
due to the small amount of Ce(SO,)2 present is 
negligible. From this it is calculated that 137,000 
r units were absorbed in forty-one minutes. 
Therefore, 6.06 X10'” ions per 137,000 7 units or 
4.42 X10" ions per r unit absorbed were reduced ; 
1.3710" r units were necessary for complete 
reduction of one mole of ceric ions. 

By definition one r unit corresponds to an 
ionization such that a current of 1 e.s.u. flows in 
1.0 cc of air under saturation potential. The 
number of ion pairs in air is therefore: 


total energy 1 





charge on electron 4774 «10-19 


= 2.1 10° ion pairs/cc 


+ Calculated from total heat equivalence of r unit, and total heat dissipated as electronic energy. 


in air and 2.1 10°-795*=1.7 X10" ion pairs per 
cc per r unit incident radiation on the water 
surface. These calculations assume a very thin 
layer of fluid so that the incident intensity is the 
same on all portions. Since this was not the case 
for the experimental conditions used, it is neces- 
sary to calculate the average incident intensity 
compared to that at the surface. Using the ab- 
sorption equation J = Je~** and allowing J» to be 
100 percent, the area under the curve of intensity 
against depth of solution can be found as 
JS.s4 100 e~#*dx. Dividing this value by 0.84 gives 
an average for all portions of the solution of 73.3 
percent of the intensity incident upon the surface. 
There were, then, 590 7/min.-0.733-1.7 10" 
=7.35X10" ion pairs per cc per min. in the 
solution under the conditions used. But, one 
microequivalent (6.0610! ions) of ceric ions 
was reduced for a forty-one minute exposure in 
10 cc; this gives: 


6.06 x 10" 
7 35X10"- 10-41 





= 2.01 ions reduced per ion pair. 


Energy of reduction in calories per equivalent of 
ceric tons.—This calculation comes out directly 
from the ions reduced per ion pair (2.01), the 
energy in electron volts to produce an ion pair 
(33 volts) and the heat equivalence per mole of 
the electron volt (23,000). The result is 379,000 
calories per mole of ceric sulfate. By an entirely 
independent method based on the best data 
available on total heat equivalence of the r unit 
(148 ergs per r unit at the given wave-length) and 
the proportion of total heat dissipated as elec- 
tronic energy (0.71), the value 344,000 calories is 


* The best value correlating the number of ion pairs in 
water and in air in terms of density difference. 
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obtained. A check within 10 percent for these 
two different calculations would seem to be 
highly satisfactory. 


Discussion of results 


Table III shows the results for ceric sulfate in 
tabular form and also values for KMnQO,, KIOs, 
and H.O:z as computed from similar calculations. 
It should be noted that all values are in terms of 
oxidation-reduction equivalents instead of moles. 
The values would be the same for ceric ion but 
would differ by a factor of 5 for KMnQ,, 6 for 
KIO; and 2 for H2Oe. This method of tabulation 
allows direct comparison of the oxidation- 
reduction change brought about by x-rays in the 
different substances. An appropriate term, ‘‘equi- 
valions,”’ has been used to indicate a single 
“equivalent” ion and bears the same relation to 
gram equivalent as does the molecule to the gram 
mole. The values given for Ce(SO4)2, KMnQO,, and 
KIO; are for amount of reduction by x-rays and 
the values for H2O2 designate the amount formed 
in 0.1. N HeSO, in the presence of dissolved 
oxygen. 

The reduction of ceric sulfate and potassium 
permanganate by x-rays is undoubtedly a com- 
plicated process involving both the formation of 
H.O2 and direct reduction by activated water. 
The importance of the activated water molecule 
in x-ray photochemistry has been repeatedly 
shown by Fricke“ and co-workers. This accounts 
for the fact that most effects of x-rays are inde- 
pendent of concentration of reactant. 

Since all solutions irradiated were open to the 
air, dissolved oxygen was present in solution and 
the production of H,O, by x-rays proceeded 
according to the over-all equation: 


2H:0+02=2H:0>. (1) 


The mechanism for this reaction has never been 
established although various ones have been sug- 
gested. Risse postulated two mechanisms for 
the decomposition of H,O by x-rays and sub- 
sequent formation of H2:O2. These mechanisms 
with added values for heats of reaction (H=heat 
evolved) as calculated from data given by Jnter- 
national Critical Tables and Lewis and von Elbe'® 
are as follows: 


CLARK AND W. S. 


COE 


2(H.O (act.)=H+OH) A= —228,800 cal. 
20H =H.Orz AH= = 34,100 cal. 
2H+0, =H.O- AH= 147,700 cal. 


Total 2H.O0+0,= 2H202 AH=-— 47,000 cal. 
or 2(H.O (act.)=2H+O) AH=-—501,400 cal. 
20 =O. AH= 163,000 cal. 
2(2H+0,. =H:202) AH= 291,400 cal. 


Total 2H2O+0,=2H.0,2 AH=-— 47,000 cal. 


These mechanisms require very large activation 
energies and also have the disadvantage of re- 
quiring collision between two groups that could 
not be present in large amounts, such as 2 OH’s 
in the first case and 2 O’s in the second. In view 
of these facts a more logical mechanism would 
seem to be: 


H;0 (act.) +O2.=H.0.+0O AH= —105,000 cal. 
H.O+0 = HO. AH= 58,000 cal. 


Total 2H.O+0,=2H.0.2 AH=-— 47,000 cal. 


This is also consistent with the fact that the 
formation of H2O, has not been observed in pure 
(Oz free) water. The activation energy involved 
is moderate and the water molecule itself is a 
reactant in each step. 

Since hydrogen peroxide reduces ceric sulfate 
in acid solution, any H2O2 formed by x-rays 
reacted immediately, according to: 


2Cet++++H2O02=2Cet**++2H*+O>2. (2) 


However, it can be seen from Table III that such 
a reaction could not account for the amount of 
reduction noted. Therefore, an additional direct 
reaction must have been caused by other activa- 
ted water molecules. The actual mechanism of 
reduction was probably a direct combination 
with secondary electrons liberated from water 
molecules by x-ray photons, but the over-all 
reaction would be: 


2Cet+++-+H,0 (act.) 
=2Ce++++2Ht+1/202. (3) 











It is recognized, of course, that the production of 
H.O; is determined in cerium-free solutions and 
that in this case part of the H2O2 produced is 
again decomposed by x-rays.!7 From the data of 
Fricke under the conditions of the experiment, 
nearly half of the HO. would be decomposed. 
Hence, the reduction of ceric ions can more nearly 
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be accounted for by allowance for the initial 
production of H,O2. Further proof for the type 
of direct reaction mentioned above is furnished by 
the fact that KIO; solutions were reduced by 
x-rays but not affected by H2Oz2 added or formed 
by irradiation. This is also consistent with the 
smaller amount of reductior noted for KIO; than 
for Ce(SO4)2. KMnQO, was reduced to an even 
greater extent than Ce(SO,)2, probably due to a 
higher efficiency of the direct reduction by 
activated water or secondary electrons. 

From Eq. (3) it is to be noted that the libera- 
tion of oxygen accompanies the process of reduc- 
tion in solutions irradiated by x-rays. Therefore, 
conditions for formation of hydrogen peroxide 
are soon established even if irradiation is executed 
in air free solutions in closed containers. A few 
exposures with freshly boiled solutions showed 
that the reaction proceeded a Jittle more slowly 
at the beginning. 

The reduction of IO;~ apparently proceeded 
only to the stage of IO~ because no free I, was 
detected by the addition of starch for the titra- 
tion until excess KI was added. The over-all 
equation would be: 


1O;-=H,0 (act.)=IO-+O02+H,0. (4) 


It is impossible to state whether the Oz liberated 
came from H,0 or IO;~. The answer may be that 
one O came from each. 

As yet, no logical explanation has been found 
to interpret the fact that BrO;~ was not reduced 
by x-rays in a manner similar to that for IO;-. It 
is even more puzzling in view of the fact that the 
reaction, BrO;- = BrO-+-Oxg, is slightly exother- 
mic while the similar one for IO3~ is endothermic. 

The decrease in amount of reduction of ceric 
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sulfate by x-rays in the presence of Ag* indicates 
that the heavy silver ion absorbed the activation 
energy from H.O before it was effective in pro- 
ducing reduction. On the other hand, Hgt+ 
seemed to act as an efficient acceptor for the 
activation energy from the water or for secondary 
electrons and brought about increased reduction 
of ceric sulfate by intermediate reduction to 
Hg2** and reoxidation by Cet++*+. 

Acetic acid also increased the amount of re- 
duction of ceric sulfate by x-rays, but as pointed 
out previously, the difference in the shape of the 
curves of reduction versus concentration of 
Hg(NOs3)2 and HC;H;O, (Fig. 3) would seem to 
indicate that x-rays produced the necessary 
activation energy to promote direct reaction be- 
tween the acetic acid and ceric sulfate. Since no 
method was available for the accurate determina- 
tion of the amount of acetic acid remaining, this 
could not be proven definitely: 
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The absorption coefficient of oxygen in the continuous band having a maximum at 1450A 
is not altered by the addition of 80 percent nitrogen at pressures less than 15 mm of mercury, 
the range covered by these experiments. The conclusion is drawn that, at least at low pressures, 
the absorption of oxygen in the Schumann region is not effected by collisions with foreign 


molecules. 





BSORPTION of ultraviolet light by oxygen 
in a band having a maximum near 1450A 
has been used as a means of inducing photo- 
chemical reactions involving oxygen. According 
to Herzberg, the absorption of a light quantum 
of wave-length shorter than about 1800A results 
in the production of a neutral and an excited 
oxygen atom. Since the interpretation of some of 
these photochemical studies depends on the as- 
sumption that the absorption of oxygen in this 
band is not altered by the presence of foreign 
molecules, an investigation of the validity of this 
assumption seemed advisable. 

The apparatus for measuring the absorption of 
oxygen consisted of a hydrogen discharge tube, a 
vacuum spectrograph and two cells of length 2.27 
and 4.75 cm equipped with lithium fluoride 
windows. The cells were introduced, one at a 
time, between the light source and the slit of the 
spectrograph. The absorption of the gas was 
measured quantitatively by comparing the in- 
tensity of the light passing through the tube when 
evacuated to a pressure less than 5X10-> mm of 
mercury with the intensity when gas was intro- 
duced at a known pressure. Further details of the 
method of photographic photometry used in 
these experiments have been given in a previous 
publication.? 

A mixture of 20 percent oxygen and 80 percent 
nitrogen was selected for investigation because 
this combination was used as a standard by G. B. 
Kistiakowsky’ in studying the reaction of oxygen 
and hydrogen induced by absorption of extreme 
ultraviolet light. Tank nitrogen when purified by 
passing over hot copper and dried with phos- 

1 Herzberg, Zeits. f. physik. Chemie 4B, 223 (1929). 


* Schneider, Phys. Rev. 45, 152 (1934). 
5 Kistiakowsky, J. Am. Chem. Soc. 52, 1868 (1930). 


phorous pentoxide was completely transparent in 
the wave-length region covered in this investiga- 
tion provided the pressure was lower than 10 cm 
in the longer cell. Since the pressure was kept 
below this value for all of the measurements, the 
absorption of the nitrogen may be neglected. 
Furthermore, since Warburg‘ has found that the 
yield of ozone per quantum of absorbed ultra- 
violet light is not changed by the addition of 
nitrogen, the nitrogen may be considered as an 
inert gas. 

Because the oxygen disappeared slowly during 
prolonged exposure to ultraviolet light, it was 
necessary to keep the gas flowing through the 
cell. The flow was sufficient to completely change 
the gas in the absorption tube five times a minute. 

Since there is negligible change in reflection of 
the light at the surface of the absorption tube 
when the pressure of the gas is altered a few 
millimeters, the intensity of the light, I, trans- 
mitted by the gas column may be expressed in the 
following form: 


[=I ¢-*r*, 


where I, is the intensity through the cell when 
highly evacuated, a is the absorption coefficient 
per unit length per unit pressure, p is the pres- 
sure, and d is the length of the gas column. 
This equation shows that a plot of log. Jo/J 
against pd should be a straight line if the ab- 
serption coefficient of oxygen is not altered by an 
increase in pressure of the gas. Fig. 1 shows such 
a plot for three wave-lengths within the con- 
tinuous absorption band of oxygen. The points 
may be seen to lie along a straight line to within 
experimental error for total gas pressures ranging 


‘Warburg, Sitz. Akad. Wiss. Wien. 216 (1912); 872 
(1914). 
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from 2.3 mm to 15.0 mm. Because of the high 
absorption at greater pressures, it was impossible 
to make accurate measurements beyond 15 mm 
of mercury. Since an increase in pressure in- 
creases the number of collisions between mole- 
cules, these experiments show that at low pres- 
sures an increase in the number of collisions has 
no effect on the absorption of oxygen in the 
Schumann region. 

The slopes of the lines in Fig. 1 give the ab- 
sorption coefficient for the nitrogen-oxygen mix- 
ture. When expressed in terms of the partial 
pressure of oxygen, the absorption coefficient for 
the maximum of the band, 1450A, is 0.55/cm/mm 
of mercury. Since the measurements were made 
at a temperature of 25° centigrade, the absorp- 
tion coefficient reduced to normal temperature 
and pressure is 460 as compared with the value 
of 490 per cm at N.T.P. obtained by Ladenburg 
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Fic. 1. The logarithm of the opacity of a mixture of 20 
percent oxygen and 80 percent nitrogen as a function of 
the pressure in mm of mercury times the thickness. (0) 
represents a layer of gas 4.75 cm thick and (+) 2.27 cm. 


and Van Voorhis® for pure oxygen. The fact that 
these two values agree to within the experimental 
error shows that the addition of nitrogen had no 
influence on the absorption of the oxygen. 


5 Ladenburg and Van Voorhis, Phys. Rev. 43, 315 (1933). 
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The purpose of this paper is to present a consistent thermodynamic treatment of solutions 


based on activity and activity coefficient. On the basis of the argument presented interrelations 
of these quantities are obtained for dissociated and undissociated components. Expressions are 


N 1878 Gibbs deduced an equation for dilute 
solutions. If this equation were to be redefined 
in terms of molar units, one of Gibbs’ constants 
identified with the gas constant, and the re- 
sulting equation be assumed applicable over the 
whole range of concentration we would have the 
ordinary ideal solution definition. In application 
it was found in general that some solutions 
approximated very closely to this equation, 
others only when very dilute, and still others 
not at all. The latter type, namely solutions of 
electrolytes, were classified separately and a 
limiting law obtained by formally patching the 
existing one. 

To extend the usefulness of this type of 
equation so that it might be applicable to all 


also derived for the fugacity and activity of a component in a gas mixture. 





solutions and reduce to it when “infinitely” 
dilute and when ‘“‘ideal,’’ Lewis and Randall 
replaced the molar concentration factor by a 
function of the independent variables and named 
it activity, which was a truly inspired idea. The 
ratio of the activity to the molar concentration 
of the species, which they named activity 
coefficient, thus automatically cataloged all 
solutions in terms of deviations from the “‘ideal 
solution.’’ This general pattern has, however, 
been partially concealed by their phraseology 
and some apparent ambiguity in defining their 
functions; as a result there have appeared in 
text books a variety of interpretations some of 
which are justifiable whereas others lead to 
mutually inconsistent interrelations of the ac- 















tivities, in particular for dissociated and undis- 
sociated solutes. 

In order to make the relations derived herein 
readily apparent the usual method of partially 
defining activity, a;, of a component “2” in a 
binary system as a function of absolute tempera- 
ture T, pressure ~, and molar concentration X; 
(such as mole fraction or molality), by the 
equation 


a;=elFi-Fi)/RT = geF i RT (1) 


will be used. Here R denotes the gas constant, 
F; (T, p, X;) the molar chemical potential (u;M; 
of Gibbs) for the state (JT, p,X;) and F,® a 
quantity that makes the mole fraction activity 
equal to the mole fraction of ‘‘z’’ if the solution 
satisfy the equation of an “‘ideal”’ solution. The 
state for which F;=F;°, i.e., a;=1, has been 
named the ‘‘standard state,’’ and so F;° could be 
called the molar chemical potential for the 
standard state. Considered from this point of 
view F,° is, by definition, a function of both 
temperature and pressure and the ‘standard 
state’ is free to vary in the (t, p) coordinate 
plane. 

If the components of a solution were miscible 
in all proportions then a single expression would 
suffice to evaluate F;° of each component but 
because complete miscibility is not in general 
obtainable in a solution F;° (or ¢;) is evaluated 
separately for solvent and solute, denoted here 
by subscripts 1 and 2, respectively. The two 
equations are, however, limiting values of a 
single relation, namely the ordinary ideal solution 
definition. 

In order to define ¢2 and thus F;° of a solute 
the dilute solution law is assumed as a physical 
hypothesis, namely! 


Xe 





e#/AT=4(T, p), (2) 
1—x2 


where x2 denotes the weight fraction of solute, 
u2 is Gibbs’ chemical potential per gram. A, 
(=R/M2), is a constant for the particular solute 
(where M, denotes the mole weight of the 
species), and ¢ (T, p) is a function of T and p 

1 This is equivalent to Gibbs’ Eq. (215). Cf. R. W. 
Goranson, Thermodynamic Relations in Multicomponent 


systems, Carnegie Inst. Washington Publ. No. 408 (1930), 
pp. 320-322. 








108 ROY W. GORANSON 


characteristic of the solute. The equation is to 
be considered as exact for the general case only 
when treated as a limiting law for x2—0. In terms 
of molar units Eq. (2), when treated as a limiting 
law, becomes 


Noe~F2/RT = $4(T, p), (3) 


where ¢2 (7, p) is a function of T and p charac- 
teristic of the substance. If Eq. (3) be extrapo- 
lated over the entire concentration range it 
becomes the usual definition of an ideal solution. 
If it be written in the logarithmic form its 
relation to the definition of activity will be 
evident. But Eq. (3) presupposes that a mo- 
lecular weight has been chosen for the component 
to which it is ascribed, the molar concentration 
and molar chemical potential being expressed in 
terms of this mole weight. Therefore since F,° 
of the solute activity will be defined in terms of 
this or a similar equation the solute activity so 
defined refers to the mole species chosen. 

Empirically it is known that for undissociated 
solutes Eq. (3) expressed as a limiting law for 
N20 exists and here the corresponding mole 
fraction activity is defined such that 


RT In g2= — F,?. (4) 


Other molar concentration units could also have 
been used, for example a common unit is mo- 
lality, m, the number of solute moles in a kilo- 
gram of solvent. The molality activity is defined 
correspondingly such that 


RT In n¢o2= — F?°, (4a) 


where 7 denotes the number of solvent moles per 
kilogram. 

F,° is therefore a function of T and p and, 
for different concentration units, will vary in 
absolute values.” 

For strong electrolytes it is known that a limit 
does not exist if the formula weight of the undis- 
sociated solute be taken as the effective mole 
weight of the species. If the electrolyte is 
assumed, however, to be completely dissociated 
into its ion parts the mean effective mole weight 
of the average ion will be not the formula weight 

2 Falkenhagen (Electrolytes, translated by R. P. Bell, 
(Oxford, 1934), p. 55) equates the molality, mole fraction, 


and moles per liter of solution activities; his activities are 
therefore not equivalent to those of Lewis and Randall. 
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M:2 of the solute but M2/v if v denotes the number 
of dissociated parts or ions. 

Let a solute RS, if it dissociates, separate into 
v parts or ions such that 


RS@TR ALS). (5) 
1 jut 


(= 


One could then formally generalize the ordinary 
dilute solution equation and write the expression 


$i (T, p) 


Nye-F kT = 


=$:(T, p) (6) 


for the 7th species. 

Herein N2=(n:/[mi+m:2]) denotes the mole 
fraction of solute where m; and m2 are the number 
of moles of solvent and solute, respectively; 
N./ denotes the mole fraction of undissociated 
solute and NV+ that of a mean species of dissocia- 
tion, or 


(1—a)N2 akNo 
” Sale tilt 1 4a(v—1)N2 








where a denotes the degree of dissociation and k 
is in general a constant for the particular solute. 
The participating or effective parts may be 
considered as (1) all distinct, k=1, (2) two 
distinct species, k=[(v,)’+(v_)’- ]}"’”, or (3) one 
mean species, k=v. It is immaterial, however, 
which viewpoint be taken because Eq. (6) may 
be divided through by & and thus all three cases 
become equivalent in treatment so long as the 
relative proportions of the species remain con- 
stant. The relationships are simplified for mo- 
lality concentrations, namely m, m’=(1—a)m, 
and ms=akm; here too the three cases will be 
made equivalent by dividing the corresponding 
equation through by k. 

In many thermodynamic problems dissociated 
and undissociated states are not differentiated 
because to do so would serve no useful purpose. 
In such cases a quantity can be taken to denote 
the average molar chemical potential of a dis- 
sociated component RS, defined here by {Fe}, 
and which is equivalent to the sum of the molar 
chemical potentials of the parts which make up 
an average mole weight of solute in the solution. 
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Thus if a denotes the degree of dissociation 


(F(T, p, N2)) =F: +X F,)= Fs for a=1, (7a) 


i=1 7=1 
=(1—a)F.,(T, p, Ne’) +avFs 
for 0<a<i, (7b) 


= F,(T, p, Ne) fora=0. (7c) 


The brace is used here to distinguish this 
quantity from the chemical potential of the un- 
dissociated part of the solute with which it might 
otherwise be confused. 

Thermodynamically one cannot distinguish 
between the participating parts but obtains only 
the mean effect of the composite. Thus, when the 
variable is the concentration of the solute, a 
measure will be obtained of {Fe} or of some 
value proportional to {F2} and not of any one 
constituent species; therefore, if there are v 
parts and @ remains unity, { F2}/v could be and 
is equivalent to the quantity generally assumed 
as the molar chemical potential of a mean ion or 
participating species. 

If the condition given by Eq. (5) exists, such 
that 0<a<1, we have F,=vFs and if this be 
combined with Eq. (7b) we have 


{ F.(T, Dp, N2)} = F.(T, p, N2’) 
=vFs(T, p, N+). (8) 


It has just been demonstrated that the molar 
chemical potential {F2} of the solute in the 
solution is equivalent to the molar chemical 
potential F; of the undissociated part of the 
solute and this again equal to the sum of the 
molar chemical potentials of the species of 
dissociation or vFs, assuming coexistence of 
these species. But unless the degree of dissocia- 
tion and therefore the number of dissociated and 
undissociated moles of solute can be determined 
the concentrations N,’ and Nz will not be 
known. Because ordinarily a is not known it is 
generally assumed that strong electrolytes are 
completely dissociated or that 


NY =0 and Ns=kN2/[1+(v—1)Ne2]. (8a) 


The ideal solution law will accordingly be 
written as 


Nae7tF2(7, p, N2))RT = 9. (T, p) (9) 
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and for solutions sufficiently dilute so that Ne 
and m can be considered as proportional to one 
another, in the form 


mse lPaRT=ngs(T, p)=o4'(T, p), (10) 


where 7 denotes the number of solvent moles per 
kilogram. 

If the assumption (8a) be made then changes 
in chemical potential which are the result of 
association will be included with other deviations 
from these solution ‘‘laws’’ and manifested as 
departures from ‘‘ideality.” 

It is known empirically that, for a solution of 
a strong electrolyte which is generally assumed 
to be completely dissociated, Eq. (9) or (10) 
exists only in the limit as Ne or m—0. The mean 
ionic activity is then given by the expression 
az(T, p, N+) =oxe!F2)'"®? and the activity of 
the undissociated solute by the expression 
ao(t, p, No’) = deel ¥2}/®T so that 


(a+)” [o+(T, p) 1’ 
ae o2(T, p) 


and for an ideal solution defined as one satis- 
fying Eqs. (3) and (9) then 


(N)’ 
—=K(T, p), 





=K(T, p) for 0<a<1, (11) 


(11a) 


which is the equivalent of the law of mass 
action. If a=1, as is ordinarily assumed for 
dilute solutions of strong electrolytes, then az is, 
by definition, zero and Eggs. (5), (7) and therefore 
(11) do not exist. Thus (a+)”’ is not equal to the 
activity d2 of the undissociated solute as has 
sometimes been stated but may be correlated 
with a quantity defined as the average activity of 
the solute only if the latter be defined by the 
expression {d2} = (a+)’. 

The mole fraction and molality® ‘‘mean ionic’”’ 
activity coefficients are defined by 


[1+(v—1)N2]ax(T, p, Nx) 
Ne 
a+(T, p, m+) 


y+(T, Pp, m+) = ’ 
m 





y+(T, p, Nz) = 





and 


3 This molality coefficient is equivalent to that of Lewis 
and Randall since here a is equivalent to their ‘‘a/k” 
where k=[(v,)’+(v_)’- }'”. (See reference 7.) 
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respectively. If the solutions furthermore behave 
according to the ideal law as defined by Eq. (9) 


1+(v—1)Ne 
1+a(v—1)Ne 


if a<1. 





y+(T, p, N+)=a 


Thus for very dilute ideal solutions the value of 
y+(T, p, N+) may be regarded as equivalent to 
the degree of dissociation. The corresponding 
molality activity coefficient is given by 


a(1— Ne) 
1+a(v—1)Ne 





vy+(T, p, m+) = ifa<l. 


It is therefore a linearly decreasing function of 
the concentration Nz and unity only in the 
limit as m—0. 
The activity of the undissociated solute in 
solution, if ideal, would be given by 
a(T, P, N2’) 
1 +a(v = 1)Ne 


=1-—y+(T, p, Nz), 


1—a 





y(T, p. N2’) _— 


ova 


which is unity for a=0 and zero for a=1; 
and by 
(1—a)(1—Ne) 


[1 +a(v— 1)N2] 





v(T, Pp, m') = 


which, for a=0, is equal to (1—N2). 

It is now a simple matter to evaluate the 
derivatives of F2°, the molar chemical potential 
of Lewis and Randall’s ‘‘standard state”’ which 
has been defined by Eqs. (1) and (4), (4a), (9), 
or (10) as a function of T and p—a simple 
method is as follows. Differentiate Eq. (1) with 
respect to the independent variables T and #, 
then let the molar concentration m converge to 
zero as a limit and we obtain‘ 


4 Lewis and Randall do not write a pressure derivative of 
solute activity, in fact their only pressure derivative (refer- 
ence 7, p. 255) relates to the activity of a simple substance. 
They do, however, say (reference 7, p. 259) of their standard 
solute state ‘‘the partial molal volume, heat content, and 
heat capacity are those of the solute at infinite dilution.” 
Their intent seems therefore evident. Some writers (e.g., 
H. S. Taylor, Ed., A Treatise on Physical Chemistry, | 
(Harned, New York, 1931), p. 761), have assumed that 
Lewis and Randall’s pressure derivative for a simple 
substance was also applicable for all activities, in which 
case F,° would not be a function of pressure. Such an 
activity could readily be written by defining F,° so that 


lim a,(T, Pp, m) a 
m—»0 a ss 


P yo* 
exp iRTv 














" 0(F2"/RT) ve(7', p,m)] v2*(T, p) 
— =tim | |-" 






THERMODYNAMIC 





0(F2°/RT) —11,(T, p, m) 
— =lim E 
oT m—0 RT? 


H2*(T, p) 
= ———_———-,_ (12a) 
RT? 


, (12b) 
ap mov0 RT RT 
where //, denotes the partial molar heat content 
and v2 the partial molar volume (=0H/dn2 and 
dV/dn2, respectively, where nz denotes the num- 
ber of solute moles, H the total heat content 
(Gibbs’ x) and V the total volume). 

The activity and activity coefficient can now 
be written as the line integral® 


y+(T, p, Ns) = 





= 


T, p, Na 11,*—I1, Vo—V9* 
=exp f 0 +———_¢p 


Fo eo. Nant pRT? vRT 
Nz 0{ Fe} 
+|- —-1 In Nz, (13) 
vRT d@Nz 


where Na=Ne/[1+(v—1)N2] and expA de- 
notes e4. This and the following line integrals 
represent the most general expressions and may 
be integrated over any path between the two 
limits. It can be broken up into three parts and 
because the limit of y+ as N20 is, by definition, 
unity for all temperatures and pressures, it is 
possible to proceed from (7 o, fo, N+.) to 
(T, p, N+) along constant T and p paths in 
which case the first two terms under the integral 
sign vanish; for example, in order to evaluate 
y+(T, p, Na) —v4(To, Po, Nao) one may inte- 
grate from (To, po, N+.) to (To, fo, N+=0) and 
then from (T, p, N#=0) to (T, p, Nx). If the 
temperature and pressure derivatives of y+ can 
be expressed as known functions of (T, p, N+) 


if the ‘‘standard state’’ were assumed to be at the fixed 
pressure p=1. L. H. Adams (Chem. Rev. 19, 1 (1936)) 
has interpreted their definition of activity as the latter 
alternative. (See reference 6.) 

° If the differential form be preferred it can be obtained 
as follows: (a) for ‘‘d In y ,.”’ delete the ‘‘exp”’ and integral 
sign. (b) for ‘‘d In a” delete the ‘‘exp”’ and integral sign 
and cancel out ‘‘d In N” or ‘‘d In m,” as the case may be, 
trom each side of the equation. This method can be used 
tor all the line integrals expressed in this paper but is not, 
nowever, applicable to line integrals in general. 
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over the region considered then the difference 
y+(T, p, Na) —v+(To, po, Nao) may be obtained 
by integrating Eq. (13) from (7, po, N+.) to 
(T, p, N+) without the necessity of evaluating 
either y1(To, po, N+,) or y+(T, p, N+). 

The molality activity coefficient may be ob- 
tained by replacing N+ by mz in Eq. (13) or by 
multiplying the right hand side of this equation 
by (1—Ne2)/[1+(v—1) Ne]. 

If the solute is not dissociated then y is re- 
placed by unity, Ns by No, asz(T, p, N+) by 
ao(T, p, Ne) and y+(T, p, N+) by y2(T, p, Ne). 
The equation then conforms with Lewis and 
Randall’s definition of the activity of an undis- 
sociated solute. 


SOLVENT 


The solvent activity is defined so that 
F;°= F(T, p) where F(T, p) denotes the molar 
free energy of the pure solvent substance. *® 

The mole fraction activity and activity coeffi- 
cient can thus be evaluated as the line integral 


a(T, p, Ni) 





vi(T, p, Ni) _ 


1 


T.pMi JI—TJI, v,;—v 
=exp [ —_——dT +——-dp 


To. we Nxt RI? RT 
N, OF, 
+] — 1 fain N,. (14) 
RT dN, 


This integral may be subdivided as was Eq. 
(13) and the same reasoning applied to it since, 
by definition, y: has the same limit unity as 
N,->1 for all values of T and p. 


6 Lewis and. Randall state ‘‘a/N=1 when N=1” and 
since a is in general a function of (7, p, Ni) this would 
restrict F;° to be a function of 7 and p only. However, they 
also state ‘‘if a pure liquid or a pure solid X at atmospheric 
pressure be chosen as the standard state of unit activity 
we may change the activity” by a change of concentration. 
This is presumably a slip but if not they should have said 

. aT, p, N1) PY 
a” rd. 
where v denotes the molar volume of the pure substance. 
It would be inconvenient for correlation purposes to have 
the solute and solvent defined unsymmetrically in this way 
and it is assumed here that the first quotation represents 
the intent of Lewis and Randall. This agrees with the 
pressure derivative as written by Schottky (W. Schottky, 
Thermodynamik (Berlin, 1929), p. 284). L. H. Adams has 
interpreted the definition of solvent activity as the second 
alternative which is symmetrical with his definition of 
solute activity. 
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If we wished to treat the solvent as a dis- 
sociated substance and consider the activity and 
activity coefficient of a mean resulting con- 
tributory part of the dissociated solvent, we 
could proceed as for a solute or by purely formal 
means substitute the symbol + for the subscript 
1’ in Eq. (14) where + denotes a mean solvent 
ion, and divide the right hand side of the equa- 
tion through by »v with the exception of unity 
in the bracket. 

We note then that, for systems miscible in all 
proportions, it is immaterial which of the two 
definitions we assign to the components because 
it is evident that the two definitions will express 
the same relation; here, however, the mole 
fraction will be the preferred unit of concentra- 
tion because as N2 approaches unity m becomes 
infinite. 


Gas MIXTURES 


Gaseous constituents are miscible in all pro- 
portions and therefore a single expression will 
suffice to evaluate the activities of these com- 
ponents. An activity and activity coefficient of a 
gas component could then be defined by Eq. (13) 
or (14). It has been customary, however, to 
define for gases a quantity called fugacity. 

An ideal binary gas mixture, by definition, 
satisfies the relations p;= Nip and F;(T, p, Ni) 
= F(T, p:), where p denotes the total pressure 
and p; the pressure constituent 7 would have if 
it existed by itself at the same temperature and 
with the same value of its chemical potential. 
F(T, pi) denotes the free energy per mole of 7 
at T and p;. These relations, however, become 
applicable in general to real gas mixtures only 
for small values of ». The physical hypothesis 
will, accordingly, be written as 


lim pNye-F kT =9,(T). 


p-0 


(15) 
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The fugacity, f;, of component i will be defined 
by Eq. (1), in which f; replaces a;, and the 
relation RT In ¢;(T) = — F,°, and may be written 
as 


f(T, P, N;) T, p, Ni 
1 mexp f 


H*—H; 
——_iT 


o p=0, 


pr; N; OF; 
+{ > —1 | In p+|— esceoel ke In N;, (16) 
RT RT dN; 


where /H° is a function of T only. 

When N;=1 then p denotes the pressure of the 
pure constituent 7 and the equation reduces to 
Lewis and Randall's’ definition of fugacity for a 
simple gas; H; and v; become // and 2, the heat 
content and volume per mole, respectively, and 
the last term vanishes. 

If the activity, a;, of a gas component 7 be 
defined by 


aT, Pp, N3) =f,(T, p, Ni) /fiT, p), 


fi(T, p) denoting the fugacity of the pure con- 
stituent, it will be evident that this activity may 
be expressed by Eq. (14). When Eq. (3) is satis- 
fied over the entire concentration range then 
a;/N;=1 and fi(T, p, N)=N.fi(T, p). When 
Eq. (15) is satisfied over the entire pressure 
range then 


fAT, p, N:)=Nip. 


In conclusion the writer wishes to acknowledge 
the helpful criticism received from L. H. Adams. 


7G. N. Lewis and M. Randall, (Thermodynamics and the 
Free Energy of Chemical Substances (New York, 1923), 
p. 192), for their definition of F° ‘‘complete the definition 
by making the fugacity of a perfect gas equal to the pres- 
sure.”’ For the argument it is not necessary to assume 
|»y—RT/p| converges to zero but merely that it remains 
finite as p converges to zero because we are concerned here 
with ratios, not differences. 
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The limitations of the simple collision theory for reactions are discussed. The Menschutkin 


and related reactions are shown to involve the formation of an activated complex with an 
electric moment of a magnitude which is calculated for the various cases. The free energy of 
activation for proton transfer is found from the experiments on conductance in water and the 
dispersion effects on polarization in ice. A partition function for liquid HO and D,O based 
upon experimental quantities is given. The abnormal rapidity of the denaturation of proteins 
in spite of an enormous activation energy is typical of reactions, involving the simultaneous 


breaking of many bonds. 


INTRODUCTION 


E shall use as the general equation for the 
specific reaction rate constant! 


KkT F,! 
Pe oe exp (—Eo/kT) (1) 
t F, 
or the equivalent forms? 


b! =——Kt = —K"( — 


h h at 


kkT 4012" * * Ay 
= —_—_g— (ANG /RT) pAS%/RE 
h at 


KkT 102° * *An 
= edie 
h at 


KkT 


KkT 
= ——e~AHT/RTeASt/R —___e-(AFt/RT), (2) 


h h 





kkT ’ kkT 109" * =") 





The superscript zero indicates the standard state 
for the rate process. The other symbols are 
largely self-explanatory and in any case have 
been repeatedly defined in the literature so that 
we shall only define them here as use is made of 
them. If leakage through the barrier is im- 
portant (but not too large) an additional factor 
1+1/24(hv,/kT)*) should be introduced into 
the above equations, where v, is the imaginary 
frequency associated with crossing the barrier.* 

‘H. Eyring, J. Chem. Phys. 3, 107 (1935). 

2a. Wynne-Jones and Eyring, J. Chem. Phys. 3, 492 
nS (7405) Evans and Polanyi, Trans. Faraday Soc. 21, 

Wigner, Zeits. f. physik, Chemie B19, 203 (1932). 
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If leakage through the potential barrier is very 
great, as for example in the rapid oscillation of 
ammonia between the two mirror image forms, 
we can best start from the specific rate expression 
in the form proposed by Tolman,‘ 


Dkip: exp (—E;/kT) 
Dip: exp (—E,/kT) 


where k;, p; and E; are specific rate, statistical 
weight and energy respectively for the 7th level. 
This can be written in the more convenient form 


EK; 
k'= Pet( Zhen exp (- ) 
kT 


kT 
+ —— exp (-E,/k7)) ’ (4) 


U 





(3) 





where we have added together the rates for the 
discrete levels and the continuum on our po- 
tential energy surface. 

An interesting example of the solvent effect in 
preventing leakage is supplied by ~CR,R2R; 
which does not racemize in solution® but in the 
gas phase would probably racemize in a time not 
very different from the ~10~'° seconds required 
for ammonia. Thus we know that this radical 
does not exist in the free state, and should not be 
thought of as uncombined with solvent or other 

4Tolman, Statistical Mechanics. (The Chemical Catalog 
Co., New York, 1927), p. 261. LaMer has used this equation 
to develop an expression for the temperature coefficient of 
the entropy of activation and related quantities; J. Chem. 


Phys. 1, 289 (1933). 
5 Jones and Wallis, J. Am. Chem. Soc. 48, 169 (1926). 
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molecules even for such very short times as 107"? 
seconds. 


The collision theory 


It has been customary to interpret the specific 
reaction rate constant for bimolecular reactions 
by the equation 


k'=PZ exp (—E/RT) (5) 


where Z is the number of collisions calculated by 
using diameters more or less appropriate for 
momentum transfer, E is RT*d 1n k’/dT, and P 
is whatever is left over to make (5) true. Clearly 
there can be no objection to this on the basis of 
its correctness. The equilibrium constant for a 
reaction of the type 


A+B=C+D (6) 
can be written as 


=—= exp ((—£;+£;)/RT), (7 
ky P;Zy ? 


PZ; 
=exp (AS/R). (8) 
f4s 


i.€., 


Here the subscript f refers to the final products 
and i to the initial reactants. The E’s are defined 
as in (5). Since, in any actual case, 


= 
as= [ AC AT/T 
0 


always depends on the temperature while Z;/Z; 
does not, then P;/P; must. Equating K of (7) to 
the statistical expression, one obtains 


K=(F;/F;) exp (—Eo/RT), (9) 


where F; and F; are the partition functions for 
the initial and final states and Ep is the excess 
energy of the final over the initial states at the 
absolute zero of temperature. The point to be 
emphasized is that (9) shows that K is inversely 
proportional to all the properties of the initial 
molecules, such as their diameters, frequencies, 
etc., which enter into the partition function Fj, 
while in (8), since the momentum transfer 
diameters are used, K is calculated as. directly 
proportional to the collision diameters of the 
initial molecules. The two equations are to this 
extent contradictory and of course (9) is correct. 
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If one tries to remedy this by taking the reaction 
rate proportional to the actual diameters in the 
activated state and inversely proportional to 
those in the normal state, as one must as indi- 
cated in (9), one can get Eq. (1), but there then 
isn’t much left of the collision theory in its 
original form. 

In spite of these logical difficulties the collision 
theory is useful if it provides a way of estimating 
reaction rates with reasonable accuracy. It cer- 
tainly will not give the rate correctly in both 
directions, however, unless it leads to the right 
value for the equilibrium constant. This is by no 
means a sufficient condition for a reaction theory 
to satisfy, but it isa necessary one. As an example 
we calculate the equilibrium constant for the 
reaction H.+J,.=2HI by both methods. Sta- 
tistical mechanics gives 


4r;4 Me 2S83—S\ — So 
K =——_- ——— exp ( ———) 
1°72? (myme)? R 


All 
Xexp (-—). (10) 


where the entropy of vibration 


S;=R 





hy;/kT —hv; 
—In (1 —exp ——) | 
exp (hv;/kT)—1 kT 


The collision theory gives 


201, 2” M3 
K= -exp (—A/I/RT). 
a3" (mym2)? 


(11) 





Here the subscripts 1, 2 and 3 refer to JJ, Js 
and HI, respectively, and 7, o, m, v, and All 
are the distance between atoms in the mole- 
cule, the collision diameter, molecular mass, 
molecular frequency and change in heat con- 
tent, respectively. At O°A Eq. (10) yields 
K=37.4 exp (—AH/RT), while, taking o1, 2°=0;° 
(which seems as good a guess as any) we get from 
equation (11) K=11.35exp(—AH/RT). At 
440°C the collision value has not changed, while 
the correct expression has been decreased by the 
term exp (—S2/R) to K=5.88 exp (—AH/RT7). 
This is not bad agreement with the collision 
theory value, and would be better at a higher 
temperature and of course much worse at lower 
ones. At constant volume AE=AH. 








10) 
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At 
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Since the activated complex is the same for 
the forward and the reverse reactions, it is 
possible to choose 1, 2?=0;? equal to some 
value rv? such that one gets correct rates in 
one direction and almost the correct ones in 
the other, i.e., within a factor of 3 at 440°C. 
At least for some reactions 7 is about the kinetic 
theory collision diameter at reaction tempera- 
tures. For other reactions such a value of 7” 
will not be even approximately correct. The 
equilibrium N2+O,.=2N0O illustrates one other 
shortcoming of the simple collision theory.®* In 
this case the entropy corresponding to the addi- 
tional factor (2+4 exp (—343.8/RT7))? must be 
introduced into (10) for the electronic states in the 
NOmolecules. Thus, at 27°C statistical mechanics 
gives the value K =297 exp (—AH/RT), while 
the collision theory gives K =2 exp (—AJI/RT). 
Certainly there is a better way of estimating 
(PZ) than the collision theory for Z and no 
theory for P. Even if we had to guess the 
properties of the activated complex, just as we 
must guess collision diameters in the kinetic 
theory, at least by using (1) we avoid the incon- 
sistencies of the collision theory. Eq. (1) neces- 
sarily leads to correct equilibrium constants for 
reactions. For further examples of large P’s in 
gas reactions see a paper by Bawn.* 

In view of these considerations it must be 
clear that the procedure of taking a standard Z 
factor for gas reactions is of very doubtful 
significance. In the liquid phase there are added 
complications. It has been recognized that for 
reactions where the reactants are ionized while 
the products are not the same Z for both reac- 
tions is not to be expected. However, it has often 
been assumed that where two reactions are of 
the same type with respect to charge the Z’s and 
the P’s should be at least nearly the same, i.e., 
for such an equilibrium there should be very 
little change in entropy (cf. Eq. (8)). This is 
occasionally, but by no means always, true as 
Table I for proton migration reactions shows. 
The data in Table I were calculated from the 
work of White and Jones.’ Many more are 
available but the above suffice for our purpose. 


The first ten of these reactions simply involve 
sean nots and Eyring, J. Am. Chem. Soc. 57, 985 
3D). 


* Bawn, Trans. Faraday Soc. 32, 178 (1936). 
‘White and Jones, Am. Chem. J. 44, 159 (1910). 
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RATES 


the transfer of a proton from some negative ion 
to one of its isomers, and should therefore fulfill 
the criteria for “similar reactions’’ as well as 
possible. AS = 4.6 e.u. corresponds to a ratio of 10 
for the factor P;Z;/P;Z;. The assumption that 
this factor should be unity (nearly) for similar 
reactions is satisfactorily fulfilled in the case of 
reaction 9, but fails by factors of from 3 to 30 
for the other ten reactions. In each of the first six 
reactions there is the alternative of hydrogen 
bonds forming between a carboxylate ion and 
some other group such as NH2, COOH, OH, 
inside the molecule or forming this bond with the 
surrounding water molecules. The first alter- 
native, if a satisfactory bond can be formed, will 
lead to the higher entropy for the dissolved acid 
molecule. This explains the higher entropies of 
the ortho over the corresponding meta ions, and 
that of the meta over the para, and of the cis 
over the trans. In the toluic acids where no 
internal hydrogen bond is possible the situation 
is reversed. Reaction (11) is about as serious a 
violation of the simple collision theory as we 
have found in this type of reaction. Here 
PZ;/P;Z;=2.5X10~. Thusif this reaction obeys 
the collision theory in one direction it is off by 
almost a million times in the other. If the same 
collision diameter be taken for reactants as for 
products then the calculated value of Z;/Z; will 
be practically unity. It is not possible here to 
throw the responsibility for the small P on some 
nonthermodynamic quantity (such, for example, 
as failure to dissipate the activation energy). The 
different rates arise, rather from the fact that 
when butyric acid ionizes, the surrounding sol- 
vent molecules have the forces holding them 
increased by much more than do the water 


TABLE I. For reactions of the type: 
acid, +ion,= ion; +acid. 








ACID ION AS | 4H 








No. K 
1 | o-amino-ben- | p-amino-benzoate 6.29 | 1912 0.94 
zoic F 
2 | metanilic sulfanilate 0.87 969 0.304 
3 | maleic fumarate 4.84 —179| 15.25 
4 | citraconic mesaconate 3.23 42 4.70 
5 | salicylic m-hydroxy-ben- 5.87 208} 13.28 
zoate 
6 | salicylic p-hydroxy-ben- 6.79 —128| 37.19 
zoate 
7 | o-toluic m-toluate —3.19| —1484 2.45 
8 | o-toluic p-toluate —3.89| —1849 3.16 
9 | n-butyric isobutyrate 0.16 26 1.034 
10 | p-toluic phenylacetate 2.28 818 0.794 
11 | n-butyric metanilate —24.51| —5777 0.077 
12 | sulfanilic hydroxide 17.62} —9775 | 6.55 X10 
13 | acetic hydroxide —3.57 | —13750 | 1.8 K 10° 
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molecules surrounding the metanilic acid when 
its ionization occurs. Thus an ice-like sheath 
gathers round the butyrate ion as it forms. This 
sheath is at a lower “‘effective’’ temperature and 
is thicker than that which forms when metanilic 
acid ionizes. This is because in the latter ion the 
charge can distribute itself over a larger volume 
through a reorganization of the bonds in the 
aromatic residue, i.e., resonance between the 
different bond structures. The larger the sphere 
over which a charge is distributed the less the 
force with which surrounding molecules are held 
and the less the decrease in entropy as the ion is 
formed. We can then state the general propo- 
sition. Any ion which can resonate between a 
variety of bond structures is provided with a 
mechanism which distributes the charge over a 
larger volume thus reducing the potential at the 
surface of the ion and therefore the binding with the 
solvent. As a result such ions exhibit a compara- 
tively low entropy of ionization. Conversely the 
entropy of ionization provides a criterion of the 
amount of internal bonding of a type not 
seriously restricting the freedom of relative 
motion of the parts which form the ion. This last 
qualification is necessary, since, for example, a 
hydrogen bond between two points on a long 
chain decreases the entropy more than if the two 
points made bonds with the solvent. 

Frequently a reaction shows marked changes 
in the “‘collision factor” with a change in solvent. 
In general this is accompanied by a compensating 
change in the activation energy with the result 
that the over all change in rate is comparatively 
small. Using the language of Eq. (2) we can say 
that an increase (or decrease) in the entropy of 
activation, ASt, is in general accompanied by a 
more or less completely compensating change in 
the heat of activation, AH?, which leaves the free 
energy of activation, AF?=AH?—TaAS}, but 
slightly changed. The explanation of this may be 
seen as follows: Suppose first we have the case in 
which two ions with unlike charge (or two polar 
molecules) combine to form a neutral molecule 
(or one with a smaller dipole moment). In polar 
solvents the ions will, of course, tend to be more 
solvated than the activated complex. This will 
result in a large value for the absolute entropy of 
activation. The point to be explained is why the 
changing of solvation by changing to a new 
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solvent will affect AF? much less than it affects 
AS*t. When a liquid solidifies by freezing there is a 
large decrease in both the entropy and the heat 
content without any change in the free energy, 
the one exactly compensating the other. Now if 
all the forces between the molecules could be 
increased, the liquid would freeze at some higher 
temperature with, however, this same compensa- 
tion between heat content and entropy. This is 
just what an ion or dipole does when it is 
introduced into a polar solvent. It increases the 
forces acting on neighboring solvent molecules, 
and if this increase is great enough, and the 
temperature not too far above the melting point 
of pure solvent, a certain number of solvent 
molecules ‘‘freeze,” thus decreasing the entropy 
of the system but giving off a compensating heat 
of solidification which more or less completely 
balances the entropy change. Thus in trans- 
ferring a mole of ions from one solvent to another 
where the forces between ions and solvent mole- 
cules will in general be different and where 
solvent melting points are different, there will, on 
the whole, be a different amount of solvent 
frozen before and after the transfer. This results 
in comparatively large changes in the partial 
molal heat content and entropy and a compara- 
tively small change in the partial molal free 
energy. Thus, if S;; is the partial molal entropy 
associated with taking a mole of reactant 7 from 
the gas phase and introducing it into solvent j, 
and if S,* is the same quantity for the activated 
complex, then the entropy of activation, A.S, in 
solvent j exceeds that in solvent k by the quantity 


AS} —AS,t=8,*—LS,;— (Si*- IS) 
=8*-Sit—L(8j—Su)-_ (12) 
Further 
P,Z;/PZu=— exp (R-'(AS—AS;4)). (13) 


Kk 


Here x; and «x, are the transmission coefficients in 
solvents j and k, respectively. If a reaction 
involves changes in multiplicity, the paramag- 
netism of the two solvents would have an im- 
portant effect on the « values. Conceivably, also, 
the x values may vary considerably due to a 
difference in the facility with which solvent 
molecules supply or absorb the energy of activa- 
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tion, but in the great majority of reactions in 
solution the important factor influencing this 
ratio seems to be exp (R-!(AS;t —AS;4)). Under- 
standing this factor, as we have seen, reduces 
simply to understanding the change in partial 
molal entropy with change in solvent for the 
activated complex and for the reactants. Equa- 
tions exactly analogous to (12) can be written for 
the difference in heat content and free energy in 
two solvents. 


AIT} —All,4 =H;* —-H.— 2 (Hi;— Hix), (14) 
AF}t —AF,A=F,*—F,*— Qo (Fij— Fix). (15) 


For the relative specific rates in the two solvents 
we obtain 


kj’ /Rx' = (x;/xx) exp ((RT)~'(— AF; + AF;*)). (16) 


It is of interest to examine the experimental 
work of Muchin, Ginsberg and Moissejewa* who 
measured the rate of the addition reactions 
pyridine to benzyl bromide and triethyl amine to 
benzyl bromide to give the two corresponding 
salts. These reactions were each measured 
in benzene-alcohol mixtures and also in ben- 
zene-nitrobenzene mixtures, where the mole 
fractions of the solvent components were varied 
from 1 to 0. The free energy of activation 


AFt=—[RTInk’—RTInkT/h] (17) 


of these are plotted against the reciprocal of the 
dielectric constant, 1/D, in Fig. 1. Curve III of 
Fig. 1 is based on the work of de Hemptinne 
and Bekaert® on the Menschutkin reaction 
((CsH;);N+C2H;Br—) measured in a benzene- 
acetone solvent, the mole fraction of both com- 
ponents being varied from 1. to 0. All the values 
are for 27°C, and the transmission coefficient, 
x, is taken as unity. The free energies of activa- 
tion are plotted against 1/D since the Debye 
Hiickel theory gives the expression 
Ne? Ne? /47e? 


3 
Ena) 
2Da 3DX\DkRT i 





(18) 


for the electrical work required to charge a mole 
of ions N. Here a, e, €, ;, 2; are the radius of the 

8’ Muchin, Ginsberg and Moissejewa, Chem. Zentr. II, 
2376 (1926). 


*de Hemptinne and Bekaert, Zeits. f. physik. Chemie 
28, 225 (1899). 





ion, its charge in electrostatic units, the charge 
on an electron, the number of ions of type 7 per 
cc, and the charge on the 7th ion in multiples of 
the unit charge. At the beginning, at least, of the 
above reactions the last term of (18) is negligible, 
there being no ions. If we extrapolate curve IV 
(Fig. 1) to 1/D=0 we get 22.07 kcal. for the free 
energy of activation, while at 1/D=1 we get 
25.65 kcal., an increase in AF? of 3.58 kcal. for 
the electrical part of the free energy of activation 
if it were carried out at unit dielectric constant, 
i.e., in the gas phase, as compared with its 
value if carried out in a solvent of infinite 
dielectric constant. We compare this value with 
Ne?/a=329.5/a kcal., the energy of charging up 
two ions with radii of a angstroms in a medium of 
unit dielectric constant. Since, if ions were 
formed, the ionic radii would not exceed a few 
angstroms and so give a very large energy of 
charging, the process of activation clearly does 
not involve the formation of ions, but does 
involve a considerable increase in the electric 
moment of the reactants. We next modify 
Kirkwood’s Eq. (26) slightly’® to obtain the 
electrical part of the free energy for a multipole, 
and combine such expressions for the reactants 
and the activated complex to obtain the free 
energy of activation, AFp;* in a solution of 
dielectric constant D, where the value AF}! is the 
value in a solution of dielectric constant 1. The 
relation is 


o (n+1)(D—1) 


AFpt = —N/2 
. /2 2 (n+1)D+n 








x| C. ES |ec+an, (19) 


b*2n+1 F b2ntl 


where the term G is the part of the free energy 
change (in going from D=1 to D=D) not 
depending on the electrostatic forces. Thus the 
relative solubilities of the reactants and the 
activated complex will change with composition 
of solvent due to changes in the van der Waals’ 
forces and chemical combination with the solvent, 
and so contribute to G. Such effects are discussed 
in detail by Hildebrand." 

10 Kirkwood, J. Chem. Phys. 2, 351 (1934). See also R. 
E. Bell, Trans. Faraday Soc. 27, 797 (1931). 


"Hildebrand, Solubility (Reinhold Publishing Co., 
New York, 1936). 
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Fic. 1. Adding ethyl alcohol or acetone to benzene 
changes the rate of the above reactions only by changing 
the dielectric constant. Nitrobenzene has an additional 
effect. 


In (19) N, D, b; and b* are Avogadro’s number, 
dielectric constant, the radius of the hypothetical 
sphere from which a reactant molecule excludes 
solvent molecules, and the same radius for the 
activated complex. The subscript i is to be 
summed over all the reactants and the asterisk 
indicates the activated complex. 

For a molecule with two charges, +e and —e, 
separated by a distance R and situated at equal 
distances, 7, from the center of a spherical 
molecule with radius 6, Kirkwood finds: 


Qo= (dei)? = 0, 

O:==ER’, 

Q2= 3u?(r?— R?/4), 

O3 = 6u?(r4 — 57?R?2/8+5.R4/48). 
mw and « are the electric moment of the molecule 
and the charge on the electron, respectively. 


(AF pt) ND? wo (2n+1)(n+1) 
a(1/D) 2 »=0((n+1)D+n)? 
G.* Qni 
ag Se 


From (20) we see that when dG/0(1/D) can be 
neglected and when D is reasonably large com- 
pared with the highest » that need be con- 


Frecipr 








S 








Ko 





Now 





(20) 


STEARN AND H. 








EYRING 


sidered, then AFp? is a linear function of 1/D. 
The reaction of benzyl bromide with triethyl 
amine and with pyridine in solvents varying 
from pure benzene to pure alcohol satisfy this 
linear relationship quite well. The same reactions 
are approximately linear up to moderate concen- 
trations of nitrobenzene in benzene. We shall 
therefore use the following modification of (20) to 
estimate u* of the activated complex from the 
slope of AF? against 1/D at the point D=5. All 
Q’s except Q; have been neglected since, in the 
reactions under consideration this is probably a 
good approximation. 


t 2772 7 #2 2 ot 
(AF pt) _ 24N*D fe M1 <] (21) 





a(1/D) (2D+1)LV* Vi Ved 


Here b, the hypothetical molecular radius, has 
been estimated by means of the equation 
8b°N = V. Using (21) we calculate the values for 
u* given in Table II. Here the subscript 1 refers 
to pyridine in I and II, and to triethyl amine in 
the other three. The subscript 2 indicates ethyl 
bromide in III and benzyl bromide in the 
remaining reactions. V* is taken as the sum of the 
molal volumes of the reactants. Without a 
knowledge of the effect of pressure on reaction 
rate, or a calculation of the potential surface, 
this is the best we can do as an estimate of |’* 
and is probably not far from right. The fifth 
column gives the amount by which the free 
energy of activation in the gas phase would 
exceed that in a solvent of infinite dielectric 
constant if 0F?/d(1/D) were constant and had 
the value possessed at D=5. If the dielectric 
constant of the solvent and the moments of the 
reactants and activated complex were the whole 
story, then simply changing the solvent should 
make no difference in the values in column 5, i.e., 
I is to be compared with II and IV with V. At 
least, at D=5 the change in dielectric constant 
seems to be the principal effect. At high concen- 
trations of nitrobenzene (see Fig. 1) some ad- 
ditional factor enters which seems best inter- 
preted as compound formation between the 
activated complex and the nitrobenzene. This 
might have been anticipated since the plot of 
dielectric constant against mole fraction of 
nitrobenzene when mixed with benzene shows a 
rise considerably more rapid than linear, indi- 
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cating association of the nitrobenzene. Such 
association may well be considerably more pro- 
nounced with a very polar substance such as the 
activated complex. Compound formation should 
be signalized by AS*t becoming rapidly more 
negative with the added nitrobenzene and AH? 
decreasing and over compensating the entropy 
change. Fig. 2 shows that this is what happens in 
the case of pyridine reacting with benzyl bromide 
in a benzene-nitrobenzene solvent. Similar effects 
are present in the other cases. The significance of 
this for the ‘‘collision number’”’ has already been 
discussed. 

The values found for yu* as given in the 
last column of Table II are near those to be 
expected for the resultant molecule. Thus the 
activated complex must have very nearly the 
same configuration as the final molecule. That 
this is true in a closely similar reaction can 
be shown by comparing the entropy of activation 
and the entropy of reaction. Thus Essex and 
Gelormini’s results’? for the reaction 


Ce6HsN(CHs3)2+CHsI =CsH;N(CHs)s1 


give ASt= —36.1 and AS= —37.06 at 60°C." 

This also indicates that the activated configu- 
ration is very like the final one. 

If ions are added to the reacting system we 
must add another term to (19) depending on the 
ionic strength. We then obtain 

N x (n+1)(D—-1) G,* Oni 
an ss 
2 n=0 (n+1)D+n | pont 7 p2ntt 





_* 
+AFi+G+RT In —, (22) 
Ily; 


where we refer back to Kirkwood’s Eq. (21) for 
the explicit form of y. 


REACTIONS INVOLVING THE BREAKING OF 
HYDROGEN Bonps!4 


Bernal and Fowler have given arguments for 
thinking of ice as possessing a structure such that 


Me Essex and Gelormini, J. Am. Chem. Soc. 48, 882 
1926). 

18 Wynne-Jones and Eyring (J. Chem. Phys. 3, 492 
(1935)) give too large a value for AS because of an incor- 
rect interpretation of the equilibrium results of Essex 
and Gelormini. 

4 W. M. Latimer and W. H. Rodebush, J. Am. Chem. 
Soc. 42, 1419 (1920). 

Bernal and Fowler, J. Chem. Phys. 1, 515 (1933). 
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RATES 119 
each oxygen is surrounded by four hydrogen 
atoms at approximately tetrahedral angles. Two 
of these hydrogen atoms are to be thought of as 
belonging to a particular oxygen atom and are 
about 0.95A from it—the normal distance in 
steam—and the other two are about 0.86A 
further away and belong to neighboring oxygens." 
This same structure is assumed to persist in 
water near the melting point and to be super- 
seded by a state of increasing disorder with 
rising temperature. Ionization results when an 
oxygen gains a third near hydrogen by one of its 
two far neighbors coming 0.86A closer, leaving an 
O-H behind. The H;*O can now lose an H+ to 
each of three neighbors. The return migration of 
the one proton simply reverses the ionization. 
However, migration of the other two corresponds 
to migration of the H;*O and by repetition of the 
process the positive charge can travel to a distant 
electrode. In a similar way O-H by acquiring 
neighboring protons can transfer the negative 
charge to any distance. This process of a negative 
ion acquiring a proton from a neutral molecule 
is the type reaction considered in Table I, so 
that the free energy of activation AF? is a 
matter of considerable interest. The migration of 
a proton from H;*O to a neutral molecule is also 
of very wide occurrence. It is therefore important 
to ascertain what value the conductance gives 
for the free energy of activation. 

If direct current is to be carried through water 
by proton transfer then since protons are liber- 
ated at one electrode and oxygen at the other, the 
process will come to a standstill unless some 
mechanism exists whereby protons after having 
passed from one oxygen to a neighboring one are 
then able to pass this oxygen. The rotation of 
molecules provides such a mechanism. Thus a 
combination of proton transfer with rotation of 














16 Pauling, J. Am. Chem. Soc. 57, 2680 (1935). 


TABLE II. 
a(aFt) 

8(1/D) Mi ue ue 
REACTION Vi Ve v* CAL. X1018 | K1018 | «1018 
II 80.5} 118.9} 200 | 5650] 2.1 1.81 | 6.52 
I 80.5} 118.9} 200 | 5050] 2.1 1.81 | 6.30 
IV 139.9 | 118.9| 259 | 3650/| 0.8 1.81 | 5.48 
V 139.9 | 118.9] 259 | 4075| 0.8 181 i 3.9% 
Ill 139.9} 76.2} 216 | 3150] 0.8 1.83 | 5.11 
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molecules constitutes one mechanism of con- 
ductance in water while the usual bodily trans- 
port of the ion is another. 

Huckel,!” Bernal and Fowler’ and Wannier!” 
have discussed this process. If the barrier which 
the proton passes over is flat a classical treatment 
is sufficient for the reaction rate. If it were high, 
since it is certainly thin, quantum-mechanical 
leakage would be important and as Bernal and 
Fowler point out the abnormally large mobility 
of the hydrogen ion would in that case not be 
present for deuterium ion. Subsequently, Lewis 
and Doody" have shown that Kt, Cl- and acid 
ion have the mobilities 64.2, 65.2 and 315.2, 
respectively, in ordinary water and 54.5, 55.3 and 
213.7 in heavy water. Thus the D;+O ion has 
nearly as great a mobility relative to K+ and Cl- 
as H;tO has so that we are dealing with a very 
flat barrier. We therefore calculate the height of 
the potential barrier neglecting quantum me- 
chanical leakage by a method similar to that of 
Huckel. Our treatment differs slightly from his 
in that we regard rotation of the molecules and 
proton transfer as two successive steps of the 
same process. The slower of the two steps is the 
rate determining one. 

We first consider the frequency k’ at which a 
proton crosses a barrier under an effective 
potential gradient of B. Using the theory of 
~ -17E, Huckel, Zeits. f. Electrochemie 34, 546 (1928). 


17a C, Wannier, Ann. d. Physik 24, 545 and 569 (1935). 
18 Lewis and Doody, J. Am. Chem. Soc. 55, 2616 (1933). 
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absolute reaction rates we obtain: 


k’ =(xkT/h) exp (—(AFt— B23,060 cos 60.43 
10-8) /RT). 


Here 0.43 X 10-° is the distance to the top of the 
energy barrier, i.e., half the distance the proton 
must travel in crossing a symmetrical barrier. 6 
is the angle which the particular hydrogen bond 
in question makes with the potential gradient. 
The factor 23,060 converts electron volts to 
calories. Effectively the distance the charge 
travels after the H;tO has rotated is from one 
oxygen center to the next, i.e. [18(6.06 x 1078)" }} 
=3.110-8A. This distance must be multiplied 
by cos @ to give the distance traveled along the 
potential gradient. The average distance per 
second along the gradient which the charge is 
carried by each proton on H;+tO is (averaging 
over all angles @) then 


(xkT/h) exp (-aF/RT) [ exp (B23,060 cos 0 
0 
x 0.43 X10-*/RT) cos 6 sin 6d0 


x( sin wi?) 
vo 


exp (—AFt/RT)3.1 10-8 23,060B 
0.43 X10-8(RT)“}. 


1 
= («kT /h) 


The three protons give three times this value for 
H;tO. The molal conductance due to proton 
transfer is then 


A’ =96,500(xkT/h) exp (—AF*/RT) 
3.10.43 X 23,060 X 10~-"B 
RT 
= 3.10610? kB exp (—AF!/RT). (23) 


x 





We take the transmission coefficient x=1. By 
using (23) with Huckel’s estimate for (Au— Ana) 
as the molal conductance due to proton transfer 
A’ where Aq and Ay, are the measured molal 
conductance of H;+O and Nat, respectively, we 
obtain Table III. 

There is considerable uncertainty as to what 
value should be taken for B the effective potential 
gradient along the potential barrier when the 
overall potential gradient is one volt per centi- 
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meter.!® Huckel assumed B = (D+ 2)/3 where D is 
the dielectric constant which varies from 88.2 at 
273°K to 47.1 at 429°K. Using this relationship 
at 298°K where D=78.8 we obtain AF? = 3350 
cal. Now if the slow process is the rotation of the 
molecules our formal expression for A’ would be 
but slightly changed and consequently we would 
obtain nearly the same value for AF?. Since we 
are measuring the slow process it seems clear 
that the free energy of activation for the reaction 
H;+O+H,O=H,.0+H;7*0 is less than 3.5 kcal. 
and lies between this value and 1 kcal. The same 
treatment of the conductance of O-H will give a 
AF? for the reaction HxyO+O-H-—O-H+H,0 
lying in this same range. The next to last row in 
Table III gives the heat of activation calculated 
from successive pairs of points if B could be taken 
equal to unity. Otherwise it must be corrected 
for the temperature dependence of B. 

We next consider the relaxation time in ice as 
revealed by the dispersion of the dielectric 
constant.2° Murphy” found that if 7 expresses the 
relaxation time of the polarized condition of the 
dielectric that the dispersion formulae obeyed in 
ice are 


€g—~ Ex 


1+w?r? 


(€,— €,,) WT 


1+w?r? 


7 


e=e,+ and ¢’= 


where ¢’ and e” are the real and the imaginary 
parts of the dielectric constant respectively. 
These are the formulae Debye showed” apply to 
rotation of molecules if the relaxation time of the 
individual molecules is given by 


€o+2 
= T. 
és +2 


Here w is 27 times the frequency of the applied 
potential e, and ¢,, are the dielectric constants for 
static and infinite frequency respectively. Murphy 
points out that the same formulae apply to a 
variety of types of polarization. 





(24) 


T1 





1 J. G. Kirkwood, J. Chem. Phys. 4, 539 (1936); Lars 
Onsager, J. Am. Chem. Soc. 58, 1486 (1936) ; Hans Mueller, 
Phys. Rev. 50, 547 (1936). 

” Errera, J. de phys. 5, 304 (1924); Smyth and Hitch- 
cock, J. Am. Chem. Soc. 54, 4631 (1932). We also want 
to express our appreciation for the opportunity of discus- 
sion with Professor Smyth. 

(1938) J. Murphy, Trans. Electrochem. Soc. 65, 133 

® Debye, Polar Molecules (Chemical Catalog Co., 1929). 
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We shall now consider what Debye’s model of 
dipoles orienting in an electric field enables us to 
say about the relaxation time. We write 


K 
ryt =k! =—e- (AF U/RT) 


h 


(25) 


the specific reaction rate for rotation. From (25) 
and Murphy’s value of 7 calculated from 
d(e’’/(e’ —e,,))/dw we obtain Table IV. 

In calculating 7; we used ¢,,=2.2 and used 
Murphy’s Table II, column 5, to get e,. There is 
the same kind of uncertainty in the formula to be 
used for calculating 7; from 7 that existed in the 
B for conductance. We, however, use the formula 
(24) proposed by Debye. Larger values of 7; of 
course would mean larger AF? values. The values 
calculated from Smyth and Hitchcock’s experi- 
ments agree quite closely with Murphy’s in spite 
of a considerable difference in the purity of the 
water used and in methods of méasurement. 

The free energy of ionization of water at the 
freezing point is AF=18.68 kcal. Hence 


(H;+O)(O-H) 
(H20)? 
and if (H;+O)=(O-H) then (H;+tO)/(H,O) 


=e (4F/2RT), Now the concentration of ions in ice 
at the freezing point is probably not very differ- 
ent from that in water. Hence AF/2 is approxi- 
mately equal to the free energy of activation AF? 
for the relaxation process. The concentration of 
activated complexes is therefore about the same 
as the concentration of H;+O and of O-H ions, 
i.e., about 10~-* to 10~’ times the concentration of 
water molecules. Our results for water show that 
proton conductance should set in for somewhat 
lower frequencies, i.e., such relaxation times will 
correspond to free energies of activation of about 
1.5 to 3.5 kcal. higher than that for the relaxation 
times measured in ice. Debye has pointed out 
that the observed polarization with an applied 





= g—(AF/RT) 








TABLE III. 

T (abs.) 273 | 291 | 298 323| 348] 373) 401) 429 
H 222.5| 314 | 350 465| 565| 644] 722] 777 
Na 26 | 43.5] 50.9] 82) 116| 155| 203) 249 

H—Na 196.5} 270.5} 299 383| 449] 489] 519] 528 
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potential of 1 volt/cm could be accounted for by 
one molecule in 510° orienting. If the ions 
present in H.O migrated toward the electrodes 
until space charge stopped them (even if there 
were no discharge) the apparent polarization and 
conductance would far exceed that observed. 
This is clear since if pairs of ions moved one space 
further apart under unit potential gradient it 
would approximately account for the observed 
polarization. 

Since the observed polarization in ice is prob- 
ably due to rotation, and therefore faster than 
proton transfer, our assumption that in water 
proton migration is slower than rotation and 
therefore is rate determining seems reasonable. 
This argument is far from conclusive, however. 
Thus the fact that HO in ice rotates with a lower 
free energy than the observed value for con- 
ductance in water does not, of course, prove that 
H,;0* ions in water also rotate with a lower free 
energy of activation. As a matter of fact, 
Wannier!”* from a comparison of his quantum- 
mechanical calculations with the experiments on 
conductance concludes that rotation is slowest. 
Whichever process is the slow one our conclusion 
that the free energy of activation for proton 
transfer is not in excess of about 3.5 kcal. is a 
necessary one. , 

Having considered the simplest reaction in- 
volving the breaking of a hydrogen bond we next 
consider another type—the vaporization of 
water. Since the activated complex is simply the 
vapor molecule, the rate of vaporization from a 
square centimeter of surface in moles per second 
is given by the expression xC,(2rmkT)', where x 
is the transmission coefficient, C, is the equi- 
librium vapor concentration in moles per cc, and 
(2xmkT)! is the average velocity normal to the 
surface.22* The only novel feature in this calcula- 
tion is the partition function used for the liquid 


2K. F. Herzfeld, J. Chem. Phys. 3, 319 (1935). 
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water, 


(2xm,kT)} 
Fi=————V (1 —exp (—500hc/kT))— 


xX (1—exp (—167hc/kT))—*Fyin. (26) 


Here m, is the effective translational free mass, 
V; is the free volume of the liquid and c is the 
velocity of light. The frequency 500 cm™ is 
Cartwright’s observed value” for the oscillation 
in which only hydrogen atoms move, interpreting 
his results in accordance with the model of Bernal 
and Fowler. He finds the corresponding fre- 
quency in D.O to be nearly 500 X 2-!. The oscilla- 
tion 167 cm is also his observed value for the 
other pair of approximately equal oscillational 
frequencies in which the effective mass is about 
that of the oxygen. He found this frequency to be 
not appreciably changed in heavy water. In our 
partition function for heavy water we take these 
frequencies to be 167 (18/20)! cm. Because 
of the way the water is constrained to oscil- 
late by its neighbors the effective transla- 
tional mass is considerably less than 18. The 
free volume is taken to be V;=b7", where 
n is found to be 5.029 from the condition that 
(d/dT)(RT°d |n F,/dT) = 18.0, the molal heat ca- 
pacity of water. For D.O we use 18.32 for the heat 
capacity per mole. This gives »=5.024. The 
value 18.32 was obtained by finding AC, for 
liquid over vapor from the vapor pressure meas- 
urements of Miles and Menzies* and adding this 
to C, for gaseous D,O assumed equal to Lewis 
and Randall’s value for H.O vapor.*> Brown, 
Barnes and Maass” find by direct measurement 
a heat capacity for D,O of 20. We are unable to 
reconcile these two results. The partition function 


23 Cartwright, Phys. Rev. 49, 470 (1936). 

% Miles and Menzies, J. Am. Chem. Soc. 58, 1067 
(1936). 

% Lewis and Randall, Thermodynamics (McGraw-Hill, 
New York, 1923), p. 80. 

26 Brown, Barnes and Maass, Can. J. Research 12, 699 
(1925); 13B, 167 (1935). 


TABLE IV. 








T (abs.) 272.3 | 270.5 | 269.3 | 265.7 | 264.7 
7X 10° 2.06 2.31 2.69 3.74 3.90 
Ti X 107 8.65 9.49) 10.5 15.5 16.1 
AFt 8390 | 8390 | 8400 | 8490 | 8470 











261.3 | 252.6) 240.6] 227.3| 266.0 | 261.1 | 227.3 
6.27; 18.9 66.4} 434 3.47 5.45} 208 
25.3 76.2 | 268 | 1750 19.3 23.2 | 964 
8620 | 8850 | 9010 | 9340 | 8610 | 8350 | 9070 
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for the gas is 
(2rmkT)*V 82?(82°ABC)*(RT)! 
Me 2h8 
X Fyin exp (—Eo/kT). (27) 





F, = F, so that we obtain for the concentration of 
light or heavy water in the gas: 


1 p 
Y ar 


uv 


m'8n2(842A BC) (kT)! 


3 
2h?(m,3b)T™ il (1—exp (—hyi/kT)—“e 





Xexp (—Eo/kT), (28) 


where p is the vapor pressure. We use the mo- 
ments for light and heavy water tabulated by 
Topley and Eyring.?” Two quantities remain to 
be determined from the vapor pressures. These 
are Ey and (m,3b). For HO these values are, 
respectively, 13042. calories and 8.43 x 10-", and 
for D,O 13426.calories and 1.35210-". Here 
m, and therefore m are in atomic units. Newton 
and Eyring?® find that for a variety of gases the 
intermediate value of } is about 10-", so that an 
effective translational mass for HzO and D.O of 
around 10 and 14, respectively, is reasonable. 
Table V shows the agreement between the ex- 
perimental and the calculated vapor pressures. 
The experimental values for D,O are those given 
by Miles and Menzies.*4 Since we have not taken 
account of holes as seems necessary to understand 
critical phenomena our equation must neces- 
sarily give too low vapor pressures at high 
temperatures.”® ; 

The abnormally large value of the entropy of 
vaporization of water arises from the fact that, 
unlike many liquids, the molecules are unable to 
rotate freely because of hydrogen bonds. The 
vapor pressure is about six times as large as for a 
normal liquid with the same heat of vaporization, 
and of course this appears as a correspondingly 
large value for AS? in the expression for the rate 
of evaporation. 





*" Topley and Eyring, J. Chem. Phys. 2, 217 (1934). 
** Newton and Eyring, Trans. Faraday Soc. In press. 
*? Eyring, J. Chem. Phys. 4, 283 (1936). 
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TABLE V. 








F ia 0 20 40 60 80 100 110 140 





4.58 | 17.54) 55.3 | 149.4) 355.1) 760 | 1074.6 | 2710.9 





obs. 
H20 
calc. | 4.596} 17.58] 55.65] 150.3) 356.9| 760.4} 1070.5 | 2708.3 


obs. 15.06} 49.13) 136.4) 331.6] 722.3) 1028.9 | 2640.4 
20 
cale. 15.10} 49.36} 136.9] 332.5| 723.9| 1026.4 | 2633.4 
































We now consider as our final example of a 
reaction with a really abnormal entropy of ac- 
tivation the denaturation of proteins. Such reac- 
tions have been discussed recently by Anson,” 
Mirsky and Pauling,*' and by Wrinch.*!¢ 

The results for the rate of denaturation of 
hemoglobin and of egg albumin, both of which 
proceed unimolecularly are presented in a com- 
prehensive survey of the subject by Chick and 
Martin. Both are first order with respect to the 
proteins and from their results the following 
tables have been prepared. 

Here the specific reaction rate is expressed in 
reciprocal seconds. 


AFt=—RT |n k'+RT In (kT/h) 


d (AF 
Allt = — r—_(—) 
dT\ T 


(Al/* — AF?) 
and ASt =—_—_____-. 


The steric factor is exp AS?/R. Clearly the ac- 
tivated complex possesses enormously greater 
randomness than does the original molecule. This 
must be due to the breaking of many bonds or the 
opening of rings or both. 

The denatured molecule as well as the acti- 
vated complex has a much greater entropy than 
the native protein as is shown by the equilibrium 
results of Anson and Mirsky on tripsin.** Thus 
their experimental results yield for the reaction 
native tripsin—denatured at 45° the values 
AF=-—-—190 cal.; AH=—67,600; AS=213 e.u. 


39 Anson, The Chemistry of Proteins, Chap. 1X, edited by 
C. L. A. Schmidt. 

31 Mirsky and Pauling, Proc. Nat. Acad. Sci. 22, 439 
(1936). 

3a FP), M. Wrinch, Nature 137, 411; 138, 241, 651 
(1936). 

82 Chick and Martin, Kolloid chemische Bethefte V, 49 
(1913). 
33 Anson and Mirsky, J. Gen. Physiology 17, 393 (1934). 
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TABLE VI. Hemoglobin. 








STERIC 
T (abs.) Rk! aFt AH? | ast FACTOR 


340.7 1.2310-% | 24,700 | 58,000} 98} 1.89107 
338.7 7.33 10-4 | 24,900 | 58,700 | 100} 5.14 10?! 
335.7 3.17X10-4 | 25,300 | 60,300 | 105 | 6.25 10” 
333.1 1.55X10-4 | 25,600 | 60,400 | 105; 6.25 x10” 


























Egg Albumin 





344.2 | 1.18 10-4 | 26,600 |135,800 | 317 | 6.66 10% 
346.1 4.41<10-4 | 25,850 |140,900 | 332} 1.20107 
349.4 | 25.4. X10-4 | 24,880 |/128,100 | 295| 1.11 10% 











Here the entropy change for the equilibrium falls 
in between the two previous values for the acti- 
vated complexes. Unfortunately there are no data 
on the equilibrium and the reaction rate for the 
same reaction. One other well established fact is 
that acid and alkali affect both the equilibria and 
rates profoundly and this effect is usually to 
lessen AF, AH and AS each in much the same 
way as it affects the three corresponding rate 
quantities. 

Mirsky and Pauling’s suggestion that dena- 
turation involves the breaking of hydrogen bonds 
provides a reasonable picture of the rate process. 
Let us assume that the 60,000 calories for AH in 
hemoglobin comes from breaking 12 hydrogen 
bonds of 5000 calories each. If these are between 
an amino and a carboxyl group then since there 
is only an oscillation of these groups before the 
bond rupture while there is rotation afterwards 
we should expect an increase in entropy. The 
rotational entropies of an amino and carboxyl 
group are 4.1 and 7.4 e.u. respectively or 11.5 for 
the sum. Twelve such pairs would give an en- 
tropy of 138 e.u. Before comparing this with the 
observed 100 e.u. of activation we should sub- 
tract something for the oscillational entropy of 
the groups before rupture and add something for 
the increased freedom given the other groups. 
All we can say then is that the increase in entropy 
seems quite reasonable on the basis of the hydro- 
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gen bond mechanism. The same calculation for 
egg albumin gives 312 e.u. to be compared with 
317 at 344.2° (abs.) with similar results at the 
other two temperatures. In estimating the en- 
tropy we have assumed the amino and carboxy] 
groups in the activated complex stage are unable 
to form new hydrogen bonds with water. In 
solutions containing acid, alkali or urea, however, 
such bonds probably do form causing A//? and 
AS? to drop greatly in such a way as to largely 
compensate each other and so give only moderate 
decreases in AF? as is observed. 

The results on the heat and entropy of acti- 
vation do not provide a basis for a decision be- 
tween the hydrogen bond theory and the cycliza- 
tion theories of Dr. Wrinch providing the latter 
bonds are not much stronger than the hydrogen 
bonds. The only requirement is that the bonds 
broken in denaturation are such as to provide a 
degree of randomness unusual in chemical reac- 
tions and more nearly comparable with the 
breaking of bonds three or four times as strong 
as those broken in the melting of ice. 

In previous papers a steric factor of about 10° 
for the unimolecular reaction by which a hexyl 
radical closes a five-membered ring simultane- 
ously expelling a CH; group has been calculated 
and compared with experiment.?* (See also refer- 
ence 6 for additional examples.) In denaturation 
we have a factor acting in the opposite direction 
making the reaction abnormally rapid by a factor 
of about 10”? and 10”, respectively. The activated 
complex theory provides a means of understand- 
ing these and of estimating the factor to be 
expected in other cases. The problems of reaction 
rates are, in fact, principally those already met in 
equilibrium problems, so that methods or in- 
formation useful in one field have an obvious 
application in the other. Each peculiar type of 
equilibrium is apt to be exhibited somewhere as 
an equally peculiar reaction rate. 

We want to thank Dr. Anson for many helpful 
discussions on denaturation. 
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The Homogeneous Unimolecular Decomposition of Gaseous Alkyl Nitrites 


VIII. The Decomposition of Ethyl and n-Propyl Nitrites at Low Pressures Together with a General 
Discussion of the Results for the Entire Series 
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The kinetics of the decomposition of ethyl and n-propyl 
nitrites have been investigated at low pressures. The 
results, together with those for the other members of the 
nitrite series, are discussed from the point of view of the 
Kassel theory. Attention is focused particularly on the 
effect of increasing molecular complexity as one ascends 
an homologous series, since the nitrites furnish an ideal 
series for such a discussion. It is concluded that the Kassel 


theory furnishes a complete explanation of the facts insofar 
as the effect of pressure on the rate is concerned. On the 
other hand, there is no satisfactory quantitative explana- 
tion of the increase in the high pressure rate with increasing 
molecular complexity. The factors which might contribute 
to such an increase are discussed, and it is concluded that 
the main factor is probably a progressive change in binding 
energy as one ascends an homologous series. 





INTRODUCTION 


ONSIDERABLE attention has been di- 
rected in recent years to the study of 
unimolecular reactions. It now appears that the 
initial step in nearly all organic decompositions 
is unimolecular in character, and that the whole 
question of the stability of gaseous organic sub- 
stances is virtually reduced to the determination 
of the unimolecular rate constants. 
The variation of the constants A and E in the 
expression for the rate constant, 


k=Ae-#!RT, 


as one passes from one member to the next in an 
homologous series has been the subject of fre- 
quent investigation during the past few years. 
In general a number of criticisms may be directed 
against the series of compounds previously in- 
vestigated. In the case of the aliphatic ethers 
the mode of decomposition varies from member 
to member of the series, and in any case recent 
work seems to indicate that the decompositions 
are much more complex than appeared at first 
sight.' The azo compounds’? in addition to 
possessing a somewhat complex decomposition 
mechanism, are further complicated by the 
presence of a side reaction of increasing promi- 
hence as one goes up the series. The esters and 





_} Rice and Sickman, J. Am. Chem. Soc. 56, 1444 (1934); 
Steacie and Solomon, J. Chem. Phys. 2, 503 (1934); 
Steacie, Hatcher and Rosenberg, ibid. 4, 220 (1936). 

* Ramsperger, J. Am. Chem. Soc. 49, 912, 1495 (1927); 
50, 714 (1928); 51, 2134 (1929). 
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paraldehydes*® are rather more satisfactory for 
comparative purposes. Here, too, however, while 
the reactions are definitely comparable, they 
involve a mechanism in which two bonds are 
broken and one other bond is converted from 
single to double. The physical significance of 
the energy of activation is therefore somewhat 
obscured. 

Recently a series of investigations on the 
decomposition of the gaseous alkyl nitrites has 
been made in this laboratory.‘ It appears that 
the primary and rate determining step in all 
cases is 


RONO=RO+NO.5 


Hence only one bond is broken in the reaction, 
and it is the same bond for every member of 
the series. The reactions have also been shown 
to be homogeneous and unimolecular over a 
wide pressure range. 

The present communication deals with an 
investigation of the decomposition of ethyl and 
n-propyl nitrites at low pressures, together with 
a general discussion of the results for the entire 
series. 


* Coffin, Can. J. Research 5, 636 (1931); 6, 417 (1932): 
7, 75 (1932); 9, 603 (1933). 

‘Steacie and Shaw, J. Chem. Phys. 2, 345 (1934); 3, 
344 (1935); Proc. Roy. Soc. (London) A146, 388 (1934); 
A151, 685 (1935); Steacie and Calder, J. Chem. Phys. 4, 
96 (1936); Steacie and Rosenberg, ibid. 4, 223 (1936): 
Steacie and Smith, ibid. 4, 504 (1936). 

(935) and Rodowskas, J. Am. Chem. Soc. 57, 350 
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EXPERIMENTAL 


Ethyl nitrite was prepared and purified under 
reduced pressure as previously described. It was 
stored as a liquid in a bulb at — 80°C during the 
investigation. 

The reaction was carried out in a constant 
volume system, and was followed by observing 
changes in pressure. The reaction vessels were 
two Pyrex bulbs of about one liter capacity. 
One of them possessed a surface volume ratio of 
about 0.6 cm~!: the other was packed with short 
lengths of Pyrex tubing, giving a surface volume 
ratio of about 3.0 cm~'. The reaction vessel was 
immersed in an electrically heated oil bath, the 
temperature of which could be maintained con- 
stant to within 0.2° by manual control. Tempera- 
tures were measured by means of standard 
mercury thermometers. 

In the higher pressure range the apparatus 
was constructed of 3/16-inch tubing, with stop- 
cock connections lubricated with Apiezon N 
grease. Pressure changes were followed with a 
differential mercury manometer, and readings 
could be taken at very frequent intervals. (See 
Table I.) At lower pressures wider tubing was 
employed, all stopcocks were replaced by mer- 
cury seals, and a set of McLeod gauges was 
substituted for the manometer. These gauges 

TABLE I. Complete data for typical runs with ethyl nitrite. 


Run No. 9. 494.6°K. Initial pressure =10.9 cm. Pressure 
measurements made with a manometer. 























Time AP k Time AP 
(sec.) (cm) (sec.~!) (sec.) (cm) (sec.71) 
0 | 0.0 —- 300 5.2 0.0029 
15 | 0.3 0.0023 360 5.7 0.0028 
30 | 0.6 0.0025 420 6.15 0.0028 
45 | 1.0 0.0026 480 6.6 0.0028 
60 | 1.35 | 0.0028 600 7.2 0.0028 
90 | 2.0 0.0029 840 8.0 0.0028 
120 | 2.65 | 0.0028 1500- 8.65 0.0024 
180 | 3.60 | 0.0029 2100 8.85 0.0025 
240 | 4.5 0.0029 3000 8.9 — 








Graphically T2s = 102 sec., Ts9 =235 sec. 


Run No. 84. 494.6°K. Initial pressure =0.0422 cm. 
Pressure measurements made with a McLeod gauge. 





0 |0.0 — 1800 0.0304 | 0.00131 
300 |0.0130)} 0.00142 2400 0.0376 —- 
600 {0.0210} 0.00136 3000 0.0376 — 


900 {0.0269} 0.00139 3600 0.0376 — 
1200 |0.0304; 0.00138 . 








Graphically T25=210 sec., 750 =510 sec. 
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were of very small volume to reduce to a mini- 
mum dead space fluctuations when they were 
cut in and out of the system. Actually each 
gauge comprised less than one percent of the 
total volume of the reacting system and about 
ten percent of the total dead space. 

A number of control runs were made to check 
the time necessary for the establishment of 
pressure equilibrium, etc. The system could be 
consistently evacuated to about 10-° cm by 
means of a two-stage steel Leybold diffusion 
pump, and no appreciable pressure change 
occurred when the system was cut off from the 
pump and allowed to stand for a time com- 
parable with the duration of a normal experi- 
ment. On a few occasions the system was flushed 
out with nitrogen before making a run, but this 
procedure had no effect upon the result. During 
a number of runs a piece of gold foil was inserted 
between the reaction vessel and the mercury 
surfaces of the gauge and cutoff, but this also 
had no observable effect on the progress of the 
reaction. It was also shown that heating the 
gauges and connecting tubing to 90°C during 
runs had no effect on the rate of reaction or 
upon the pressure increase at completion. 
Various samples of ethyl nitrite were used, and 
these gave identical results. 

In carrying out an experiment the reaction 
system was first cut off from the pump. A suitable 
amount of ethyl nitrite vapor was then expanded 
into the reaction system by opening the stopcock 
or mercury cutoff leading to the reservoir for 
seven seconds. The initial pressure reading was 
made eight seconds later, and further readings 
were taken at intervals until the pressure ceased 
to change. In runs made with a McLeod gauge, 
the gauge was kept connected to the reaction 
system between observations, and was cut off 
only when readings were being taken. It was, of 
course, impossible to make readings as frequently 
with a McLeod gauge as with a manometer. 

Propyl nitrite was prepared and stored in a 
manner analogous to that of ethyl nitrite. 
Several samples were prepared at normal and at 
reduced pressures, and their purity was checked 
by boiling point determinations. All runs were 
made at pressures below 3 cm, using McLeod 
gauges. To eliminate the possibility of the con- 
densation of the products of reaction, all experi- 
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ments were made with the gauges thermostated 
at 100°C and with the connecting tubing heated 
electrically to the same temperature. Runs made 
with the packed reaction vessel did not deviate 
appreciably from those with the empty vessel. 
In spite of all precautions, the experiments with 
propyl nitrite were much less reproducible than 
those with ethyl nitrite. There seemed to be no 
apparent explanation of this fact. 


EXPERIMENTAL RESULTS 


All runs were carried to completion. The 
times for 50 percent, and 25 percent, of the total 
pressure increase to occur were obtained graphi- 
cally, and were used as a measure of the rate of 
reaction. 


I. Ethyl nitrite 


In all about 100 runs were made at three 
different temperatures. Complete data for two 
runs in the empty reaction vessel are given in 
Table I. The data for all runs are given in 
Figs. 1 and 2, in the form of log P—log 104/T so 
and log P—log 10*/T»; graphs, and the data for 
a series of representative runs are given in 
Table II. 

The mean pressure increase for all runs was 
84.8 percent, in excellent agreement with Steacie 
and Shaw’s high pressure value of 86 percent. 
The final pressure was quite sharply defined, 
and showed no tendency to fall off slowly as 
was reported in the case of methy] nitrite. 

For the sake of comparison, Steacie and Shaw’s 
high pressure rates are plotted in Fig. 1 (indi- 
cated by double circles). The agreement is 
remarkably good. It may be observed that 
Steacie and Calder reported rates slightly lower 
than those of Steacie and Shaw in the case of 
methyl nitrite, and this is also the case here at 
the one temperature at which the difference was 
appreciable. 

Runs made with the packed reaction vessel are 
indicated by solid circles in Fig. 1. It is obvious 
that the reaction is not appreciably hetero- 
geneous down to pressures of 10-? cm. 


II. n-Propyl nitrite 


The experiments with propyl nitrite were all 
made at 483.6°K. Considerable difficulty was 
experienced in obtaining reproducible data, and, 
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Fic. 1. The variation of 7s» with pressure for ethyl nitrite. 


as Table III and Fig. 3 indicate, some uncer- 
tainty attaches to this work. The results are, 
however, accurate enough to establish fairly 
definitely the pressure at which the rate begins 
to fall off appreciably. As may be seen from 
Fig. 3, the agreement with the results of Steacie 
and Shaw is well within the experimental error. 
Packing the reaction vessel had no apparent 
effect on the results. 


DISCUSSION 


The most important characteristic of a uni- 
molecular reaction is the falling off in rate at 
low pressures. As predicted by the theory, ethyl 
nitrite maintains its high pressure rate down to 
lower pressures than does methy! nitrite, and 
the still more complex propy! nitrite maintains 
it to even lower pressures. 

We shall now proceed to a discussion of the 
effect of pressure on the rate for the alkyl 
nitrites from the standpoint of the Kassel 
theory.® Kassel assumes that if a molecule has j 
quanta divided among s oscillators it will react 
when m quanta are localized in a particular 
bond. He then gets directly for the rate of 
reaction at high pressures 


o /j+s—1 
pee adnewnry (HAN) 


=m s—1 





j\G—m+s-—1)! 
Xx . . 
(j-—m)'(j+s—1)! 


emake kT (1) 


6 Kassel, Kinetics of Homogeneous Gas Reactions (New 
York, 1932), p. 96 et seq. 
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Fic. 2. The variation of T25 with pressure for ethyl nitrite. 


which, on carrying out the summation reduces to 
k= Ae mir, (2) 


Here A has the general character of the reciprocal 
of a relaxation time, and y is the frequency of 
the oscillators, all of which are, for simplicity, 
assumed to be identical. 

At low pressures we must consider the number 
of collisions, etc., and for k, the rate at low 


pressures, we get 
—~ ') 
e~ (i-—m) hv/ kT 


j-m 





k=k,(1 —e~WlkT)s 7 
j=m A (j-m+s—1)!j! 


aN (j—m)\(s—1)! (3) 





where a, the collision factor, =4(rkT/mo)?-o°, 
mo is the molecular weight, and o the molecular 
diameter. 

In applying this relationship, we put mhvNo 
=E, where E is the energy of activation, and 
No is Avogadro’s number. The value of v is 
adjusted so as to give agreement between the 
observed and calculated values of k/k,, at 
various pressures, although, of course, it must be 
roughly consistent with known vibrational fre- 
quencies of the molecule. The number of effective 
oscillators in the molecule, s, is taken as the total 
number (37—6, where n is the number of 
atoms) less the number of C—H bonds, since 
these possess too high a frequency to be thermally 
excited at the temperatures under consideration. 

As a result of the series of investigations on 
the decomposition of the alkyl nitrites, we now 
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have available for the first time data for an 
homologous series of compounds in which (a) we 
are sure that exactly the same primary reaction 
is taking place in all cases, and (b) the sub- 
stances are simple enough to enable us to 
determine the falling off in rate at low pressures 
for several members of the series. 

For what we might call an ideal series of 
homologues, the same bond would be broken in 
every case, and the bond strength would be the 
same throughout the series. A methyl derivative 
possessing definite values of m and s should be 
followed by an ethyl derivative with met= Mie, 
and Set=Sme+7, and a propyl derivative with 
Myr=Mme, ANd Spr=Smet 14, etc. Furthermore E£, 
since it is equal to Nomhv, should be constant 
throughout the series. Hence, from Eq. (2), 
k,, should be the same for all members of the 


TABLE II. Summarized data for representative runs with 
ethyl niirite. 








Po Final T2s5 Ts0 Po Final “e5 ) 
(cm) % AP | (sec.) | (sec.) (cm) % AP | (sec.) (sec.) 





Manometer 494.6°K 





21.25 | 35.1 96 | 230 ||4.15 87.9 126 | 312 























15.3 83.7 | 102 | 235 ||2.05 110.0; 108) 285 
10.6 83.9 | 102 | 235 |/0.95 106.0 90; 300 
7.3 84.9 | 108 | 255 | 





McLeod Gauge 494.6°K 





1.12 | 74.2 | 135 | 300 |lo.0153 | 80.4] 180| 630 
0.205 | 78.5 | 180 | 390 |l0.0093 | 89.1| 420) 930 
0.146 | 74.4 | 120 | 270 |l0.0061 | 81.8! 390) 1080 
0.0835| 73.9 | 165 | 420 |l0.00207| 86.6| 1020 | 1650 


0.0283) 84.1 | 240 | 570 | 





Packed Bulb 494.6°K 









































17.1 84.8 | 126 | 300 ||0.0332 | 90.3} 180| 480 
6.9 | 92.1 | 135 | 330 ||0.00536| 84.0} 420| 930 
4.45 | 96.6 | 126 | 300 ||0.00343| 84.0| 540 | 1020 
0.344 | 83.7 | 180 | 360 

483.6°K 
16.3 | 87.1 | 258 | 642 10.321 | 79.0} 276| 734 
8.5 | 88.2 | 252 | 645 ||0.0506 | 72.7| 390) 1050 
4.85 | 94.8 | 258 | 642 ||0.0032 | 103 | 1530 | 3510 
0.996 | 78.9 | 285 | 720 

505.0°K 
15.5 86.1 | — | 114 0.146 70.0 72 
4.8 | 82.3 | — | 123 110.0257 | 75.1) — | 125 
0.906 | 71.2} — 84.3 228 








90 | 0.0083 
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series, provided that A does not vary. This 
latter probability is discussed later. 

The results obtained for the high pressure rate 
constants of the nitrite decompositions are: 


Rmethy) = 1.84 10"3 e-36400/RT cee. —1, 
Rethy! = 1.39 10"4 e-37700/ RT 
Rocenyt =2.75X10"4 e~37650/RT 


Pise-prepyt _ 1.26X 10" ree 


Ryutyi = 3.02 K 10 e—37000/ RT, 


The variations in E are within the experimental 
error, but their interdependence with the non- 
exponential term obscures any regularities in 
the rates. For the moment let us arbitrarily put 
E=37000 calories per gram mole for all the 
reactions, and recalculate A. We get 


Substance A A/Ame 
Methyl nitrite 0.33 XK 10" (1.0) 


Ethyl nitrite 0.65 X 10" 2.0 
n-propyl nitrite 1.3510" 4.1 
iso-propyl nitrite 1.2610" 3.8 
n-butyl nitrite 3.02 X 104 9.1 


In other words the rate approximately doubles 
for each added CHe group. This is in general 
agreement with Coffin’s results. It may also be 
noted that within the experimental error the 
rates for n-propyl and iso-propyl nitrites are 
identical. 

We now come to the consideration of the low 
pressure rates. Fig. 4 gives the experimental 
curve for methyl nitrite together with the pre- 
dictions of the Kassel theory for ¢= 3.0 10-* cm, 
s=12, m=13 and 14, and a temperature of 
494.3°K. The value used for A in these calcula- 
tions was one-half the experimental value, since 
it appears certain that in.the nitrite decom- 
positions two molecules disappear for each 


TasBLe III. Summarized data for experiments with propyl 
nitrite. 


Temperature = 483.6°K. 























Po Final Tso | Po Final Tso 
m | YAP (sec.) | (cm) % AP (sec.) 
1.066 | 70.5 240 || 0.121 77 270 
0.709 | 68 285 || 0.116 78 315 
0314 | 72.5 360 || 0.0920 90 270 
0.275 | 94 240 =|! 0.0420 83 555 
0.210 | 66 270 || 0.0415 95 270 
0.182 | 89 270 ‘|| 0.0395 97 360 
0.125 | 88 270 | 0.0320! 106 420 
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Fic. 3. The variation of T;» with pressure for propyl nitrite. 


elementary act.’ The agreement of theory with 
experiment seems quite satisfactory with .m 
equal to either 13 or 14. 

Fig. 5 gives similar results for ethyl nitrite, 
s being taken as 19, and other factors as shown in 
the figure. It is seen that excellent agreement is 
obtained if we use the same values for m and o 
as were used with methy] nitrite. 

A similar calculation for propyl nitrite, using 
s=26 and other factors as with methyl nitrite, 
shows that the falling off in rate should not be 
appreciable above 0.1 cm. This is consistent with 
the experimental results given in Fig. 3. 

It may therefore be concluded that the Kassel 
theory gives excellent agreement with the experi- 
mental facts. In this connection it should be 
emphasized that all adjustable parameters are 
exhausted in fitting the theory to the results 
for methyl nitrite, and that there are no fui ther 
opportunities for adjustment of the theoretical 
values for ethyl and propyl nitrites. The nitrites, 
therefore, approximate very closely to the ideal 
homologous series discussed above. The theory, 
however, offers no explanation of the slight but 
regular increase in the high pressure rates as one 
ascends the series. 

In looking for an explanation of this effect, it is 
obvious in the first place that the diameter and 
mass of the molecule, steric effects, etc. can have 
no effect since they all vanish in the expression 
for the high pressure rate. A variation in the 
proportionality factor, A, would account for the 
observed effect. While the significance of A is 
rather hazy,’ it has the general meaning of the 


7 See Kassel, Chem. Rev. 10, 11 (1932). 
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reciprocal of a relaxation time. It would therefore 
be expected that it would diminish with in- 
creasing molecular complexity, which is the 
reverse of the experimental observations. It 
therefore seems unlikely that the explanation of 
the facts is to be found here. 

There is evidence that bond strengths may 
diminish as we ascend an homologous series. 
Thus Ramsperger has reported a progressive 
decrease in the activation energies of the azo 
compounds. Rice and Johnson? in discussing the 
strength of the C—C bond in saturated hydro- 
carbons report a progressive decrease in the 
energy of activation for the split into free 
radicals with increasing molecular weight. There 
is, of course, considerable doubt regarding the 
relationship between the low temperature de- 
composition and the free radical decomposition 
of the saturated hydrocarbons, but it appears 
almost certain that the primary step in the 
nitrite decompositions is a direct split into 
radicals, and is therefore directly comparable 
with Rice and Johnson’s work. 

It seems therefore that the most probable 
explanation of the increase in the high pressure 
rate as one ascends the nitrite series is to ascribe 
it to small changes in E, A probably remaining 
roughly constant. Since E=Nomhyv, it follows 
that either m or v must change. A change in m of 1 
would correspond to a change in E of 2580 
calories, while a change of 600 calories from 
member to member is all that is required to bring 


8 T, Am, Chem, Soc, 56, 214 (1934), 
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the high pressure rates into agreement. Of course 
the ‘“‘quanta” involved are highly artificial in 
character and fractional changes might be per- 
mitted. The most simple and logical assumption, 
however, is that there is a slight progressive 
decrease in v as one ascends the series. Assuming 
m= 14 for methy] nitrite in the above calculations 
is equivalent to putting hy=2580 calories per 
mole, whence v=910 cm. All that would be 
required to bring the high pressure rates into line 
would be a decrease in v of 15 cm™ from member 
to member. Such a decrease is by no means 
unreasonable from the point of view of molecular 
structure. 

Rice and Gershinowitz® on the basis of the 
‘freezing out” of degrees of freedom, have 
recently calculated the A factors for methyl and 
ethyl nitrites. They find the same value for both. 
The calculation, however, is based on the as- 
sumption of the same frequency in both cases for 
the bond which is about to break, and a small 
change in this frequency from member to member 
of the series could materially change the 4 
factor. 

It follows that the whole situation is compli- 
cated in that a small change in the frequency of 
the pertinent bond may be expected to change 
both A and E. At the moment there is no theory 
sufficiently detailed to predict accurately what 
will happen as one ascends an homologous series, 
but qualitatively at least it appears that the 
assumption of small changes in binding strengths 
is a sufficient explanation of the facts. 


9J. Chem. Phys. 3, 479 (1935). 
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A Study of Long Chain Organic Compounds by Electron Diffraction 


K. H. Storks Anp L. H. GeERMER 
Bell Telephone Laboratories, New York City 
(Received November 28, 1936) 


A layer of lard deposited upon a metal surface gives an 
electron diffraction pattern consisting of segments of 
diffuse inclined bands. These arise from line gratings made 
up of long hydrocarbon chains. The mean inclination of 
the axes of these chains to the surface normal is calculated 
to be 27°, and the separation of carbon atoms 2.50A. The 
saturated glyceride, tri-stearin, vields a diffraction pattern 
of diffuse bands parallel to the surface, unlike the diffrac- 
tion pattern from lard, although the latter is predominantly 
a mixture of glycerides. We believe that the inclination of 
the long hydrocarbon chains in lard is due to lack of 
saturation, although this hypothesis has not been tested. 
Diffraction patterns from stearic acid deposited upon a 


EVERAL experimenters have examined long 
chain organic compounds by electron diffrac- 
tion.' The first to report an investigation of this 
kind were Thomson and Murison who studied 
long chain greases of the paraffin type deposited 
on metal surfaces. They, as well as later workers, 
obtained diffraction patterns consisting of several 
diffuse parallel and equally spaced bands on 
which appear a few fuzzy spots. These patterns 
are attributed to extremely imperfect crystals 
containing long hydrocarbon chains standing ap- 
proximately normal to the supporting surface. 
In general the materials studied in this manner 
have been mixtures of saturated single chain 
compounds containing carbon and hydrogen 
only. Although crystallographic data on such 
compounds are not extensive, they indicate that 
these often, if not always, crystallize in the 
orthorhombic system with the axis of the long 
chain parallel with one of the crystallographic 
axes. The observations by electron diffraction 
have been consistent with this type of structure, 
and from these observations it is deduced that 
paraffin type hydrocarbon compounds orient 
themselves upon a metal surface with the two 
short crystal axes lying in the plane of the 
surface. 


‘ The earliest papers are: Thomson and Murison, Nature 
131, 237 (1933); Nelson, Phys. Rev. 44, 717 (1933); 
Murison, Phil. Mag. 17, 201 (1934). 


metal surface also contain inclined bands. These are 
relatively well defined, and have inclinations which differ 
from one sample to another. This variation is attributed 
to the existence of large crystals the c axes of which are 
differently inclined to the plane of incidence. We calculate 
that this axis is always inclined 33° to the surface normal. 
A similar investigation of methyl stearate vields the angle 
25°. For both stearic acid and methyl stearate the separa- 
tion between diffraction bands corresponds to the atomic 
spacing 2.50A. This is in good agreement with the known 
distance between alternate carbon atoms in a zigzag 
hydrocarbon chain. 


LARD 


We decided to study an entirely different type 
of greasy compound. For this we chose com- 
mercial lard as a typical grease mixture of animal 
origin. Samples of lard were prepared on blocks 
of polycrystalline nickel which had been carefully 
polished and freshly cleaned. Diffraction patterns 
were obtained from these by scattering high 
speed electrons at very small glancing angles. 
Essentially similar patterns were found when 
lard was deposited by evaporation from benzene 
solution, when it was simply rubbed upon the 
metal surface, and when it was fused and allowed 
to flow freely. In still another experiment lard 
was deposited upon the naturally formed face of 





Fic. 1. Electron diffraction pattern from a layer of lard 
upon a block of polished nickel. (4 =0.051A. 
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a single crystal of carborundum, and in this case 
also the pattern was like those found earlier. 

A typical diffraction pattern obtained from 
lard is shown in Fig. 1. The diffuse smears of this 
figure are portions of two sets of parallel bands, 
the members of each set being equally spaced and 
inclined at angles of about 27° to the plane of the 
specimen surface. Except for this inclination the 
features of Fig. 1 are rather like those observed 
for paraffin, and we interpret them in an an- 
alogous manner. The two sets of bands corre- 
spond to hydrocarbon chains which are inclined 
to the surface normal by angles within a range 
about 27° but are otherwise unrelated. 

To see how this comes about we consider first 
the diffraction effect of a single chain of carbon 
atoms. We designate by d the separation of atoms 
measured along the direction of the chain, by a 
the angular inclination of the chain to the surface 
normal, and by ¢ the azimuth of the chain 
measured from the plane of incidence. This chain 
of carbon atoms would give rise to diffraction 
cones intersecting the photographic plate in lines 
approximately straight, parallel and equally 
spaced.” These lines are described by the equation 


1X 
y=n———-+ * tan asin ¢ (1) 
d cosa 


where x and y are coordinates in a Cartesian 
system with its origin at the intersection of the 
primary beam with the photographic plate and 
its x axis parallel to the supporting surface. (A is 

2 Except for the very special case of @ near 7/2 and ¢ 
llear zero. 


AND 


Fic. 2. Diffraction patterns from stearic acid. (A=0.051A.) 
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the electron wave-length, 1 the distance from 


specimen to photographic plate and 7 any whole 


number. ) 

The common slope of the lines is dy/dx=tan 
asin ¢, and their intersections with the y axis are 
nLX/d cos a. We note that if ¢ is varied while a 
remains fixed the lines are turned about these 
fixed intersections. The slopes can be varied from 
—tan a to +tan a; the angles which the lines 
make with the x axis are varied from —a to +a. 
As the azimuth plane of the molecule is rotated 
the diffraction lines sweep out symmetrical pairs 
of sectors of angular width 2a. If all the mole- 
cules in the diffracting material have the same 
value of a but are randomly disposed in azimuth 
these sectors will appear on the pattern. They 
will not be uniformly intense across their widths, 
but most intense at the edges and least intense at 
the centers. 

The diffuse bands of Fig. 1 are the edges of 
these predicted sectors. Nevertheless only short 





lic. 3. Diffraction pattern from methyl stearate. 
(A=0.051A.) 
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sections of the bands appear, and their behavior 
at the y axis is not that to be expected from the 
analysis just given. We must assume that a also 
is somewhat variable from one chain to another. 
From Eq. (1) one finds that, for g= 2/2 which is 
the most important case, as a@ varies a line will 
remain at a fixed distance from the origin and its 
inclination will remain equal to a. 

We believe that the observed diffraction 
pattern arises from hydrocarbon chains dis- 
tributed randomly in azimuth and inclined at 
somewhat differing angles to the surface normal. 
Because of the complex nature of lard and the 
extreme fuzziness of the features of Fig. 1 it 
seems unprofitable to attempt to carry the an- 
alysis further. We can, however, conclude that 
most of the hydrocarbon chains in lard are in- 
clined to the surface normal by angles in the 
neighborhood of 27°. The separation of carbon 
atoms, d, can also be calculated. From Eq. (1) 
the distance between adjacent bands is given by 


nL» 
D= (2) 
d(1—cos? ¢ sin? a)! 





As the diffuse bands represent the edges of the 
sectors we have g=7/2, and Eq..(2) allows us to 
calculate d=2.50A. 

We have investigated some of the compounds 
which might be in lard. As it is known to be a 
mixture of glycerides of long chain organic acids 
plus small amounts of free acids we chose to study 
stearic acid and tri-stearin, a glyceride of stearic 
acid. These compounds are not exactly typical of 
the constituents of lard because the chains of 
both are completely saturated whereas the 
natural product is known to contain a con- 
siderable proportion of unsaturated chains. 

We have found that tri-stearin produces a 
diffraction pattern of straight lines lying parallel 
to the supporting surface, entirely unlike the 
pattern obtained from lard. Separations of these 
lines correspond again to the spacing 2.50A. 
Lines of even order are strong and lines of odd 
order weak. From this pattern we predict that 
the structure of tri-stearin will probably be found 
to be orthorhombic. 

Examples of diffraction patterns from stearic 
acid are shown in Fig. 2. Although these exhibit 
inclined lines and resemble the pattern from lard 
in a rather striking fashion, it is known that lard 





contains little free acid. We were then led to test 
a mixture made up of tri-stearin and a small 
percentage of free stearic acid. This mixture gave 
the diffraction pattern of tri-stearin, without any 
observable modification due to the acid. It is 
interesting to note that the mixture is very 
different from lard mechanically. We now believe 
that the inclination of the long hydrocarbon 
chains in lard is due to lack of saturation of these 
chains in the acid radicals in the glycerides. This 
hypothesis has not been tested as we have been 
unable to obtain glycerides of unsaturated long 
chain acids which are solid at room temperature. 


STEARIC ACID 


The diffraction patterns of Fig. 2 were pro- 
duced by stearic acid deposited upon polished 
nickel by evaporation from benzene solution. 
Although these patterns exhibit fairly well de- 
fined spots characteristic of single crystals, they 
also show the inclined bands which are attributed 
to the separate molecules. We conclude that the 
crystals are extremely imperfect, so imperfect 
that diffraction patterns from the independent 
line gratings of separate molecules are still 
evident. 

Microscopic examination shows that layers of 
stearic acid have a strong tendency to form large 
crystals. Differences in disposition and size of 
these crystals undoubtedly account for the differ- 
ences between the patterns of bands in Fig. 2. 

The first of these, Fig. 2A, resembles the lard 
pattern. The equal intensities and inclinations of 
the two sets of bands in this figure indicate that 
the stearic acid crystals are more or less ran- 
domly disposed in azimuth. The actual inclina- 
tions and spacings of the observed bands will 
then correspond to long hydrocarbon chains lying 
in the two azimuths normal to the primary beam. 
That this assumption is true is borne out by the 
fact that most stearic acid diffraction patterns 
exhibit bands having approximately the same 
inclination as those of Fig. 2A. This inclination 
occurs most frequently although substantially 
smaller inclinations are sometimes found, and 
even bands parallel to the surface are observed in 
rare cases. 

The single set of bands of Fig. 2B has the same 
inclination as the two sets of Fig. 2A. This single 
set of bands corresponds to hydrocarbon chains 











lying in one azimuth normal to the primary beam 
direction, or at least in azimuths within about 25° 
(estimated by means of Eq. (1)) of this normal 
position. The lack of sharpness of the features of 
Fig. 2B suggests that the pattern does, in fact, 
arise from many crystals the hydrocarbon chains 
of which are located in different azimuths. The 
bands of Figs. 2A and 2B yield the values 
d=2.50A and a= 33°. 

The sharp features of Fig. 2C suggest that this 
pattern is produced by a single imperfect crystal. 
That this is true is evidenced further by the fact 
that the inclination of the lines of Fig. 2C is only 
24°. From the relation dy/dx=tan a sin ¢ we can 
calculate the azimuth, ¢, which corresponds to 
the slope dy/dx = tan 24° assuming, of course, the 
value a= 33° obtained from Figs. 2A and 2B. 
This calculation yields g=45° (or 135°) as the 
azimuth of the hydrocarbon chains producing the 
pattern of Fig. 2C. If this value of ¢ is substituted 
into Eq. (2) together with a= 33° and d=2.50A 
we obtain D=8.6 mm as the spacing of the lines 
of Fig. 2C. Direct measurement gives 8.5 mm. 
(The separations in Figs. 2A and 2B, from which 
the value 2.50A was calculated, are 8.0 mm.) 

It is interesting to compare the angle a= 33° 
and the spacing d=2.50A, obtained from these 
experiments, with the results of measurements 
upon single crystals of stearic acid. Miiller® gives 
2.51A as the distance between alternate carbon 
atoms in the zigzag chain of the acid and, for the 
68 modification, he has found the angle between 
the c axis and the normal to the basal plane to be 
26° 22’. For the a@ modification’ this angle is 
reported to be 39° 10’. New and, as yet, unpub- 
lished electron diffraction measurements® by us 
upon single crystals of stearic acid formed by 
evaporation from benzene solution have shown 
that these crystals orient themselves upon any 
surface with the two shorter axes in the plane of 
the surface. These measurements show, further- 
more, that the shorter axes are at right angles 
and have the lengths reported by Thibaud and 
Dupré la Tour. (We have not yet succeeded in 
determining a possible length and inclination for 
the ¢ axis.) From these unpublished experiments 
we feel justified in believing that the very imper- 

3 Proc. Roy. Soc. A114, 542-561 (1927). 

4 Thibaud ard Dupré la Tour, J. Chim. Phys. 29, 153- 


167 (1932). 
5 See, however, Storks, Phys. Rev. 50, 396 (1936). 
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fect stearic acid crystals, with which we are con- 
cerned here, have the form of the a modification. 
We are unable to reconcile our angle 33° with the 
reported angle 39° 10’. It is, of course, possible 
that Thibaud and Dupré la Tour chose a unit 
cell in which the c axis is not parallel to the long 
hydrocarbon chain. No other choice would, how- 
ever, have given an angle in agreement with ours 
without making the c axis too short to contain 
two molecules end-to-end, which is apparently 
the way they must be arranged. 

We have given reasons above for believing that 
the pattern of Fig. 2C is formed by a single 
crystal, and we have calculated its orientation. 
We ought then to be able to account for the 
rather well defined spot pattern of this figure. In 
this we have not succeeded. Even if our data 
could be reconciled with those of Thibaud and 
Dupré la Tour, there would not be sufficient 
information at hand regarding relative intensities 
of different reflections from a crystal of the a 
modification. We are now engaged in studies of 
single crystals of this sort by the transmission 
method,* and hope soon to be able to clear up the 
present uncertainties. 

It has seemed important to us to test whether 
or not the inclination of the hydrocarbon chains 
is the same in various stearic esters as in stearic 
acid itself. Microscopical examination of layers of 
methyl stearate deposited by evaporation from 
benzene solution has not disclosed any observable 
crystals. We conclude that they must be ex- 
tremely small. This conclusion is in accord with 
the fact that diffraction patterns from different 
samples of methyl stearate are similar and indi- 
cate hydrocarbon chains randomly distributed in 
azimuth. Such a pattern is shown in Fig. 3. This 
yields d=2.50A and a=25°. 

The inclinations of the hydrocarbon chains of 
stearic acid (33°) and of methyl stearate (25°) 
differ by an amount which we believe greater 
than experimental error. This difference suggests 
that if we were to study higher single chain stearic 
esters we should find hydrocarbon chains stand- 
ing more nearly normal to the supporting surface. 
We propose to carry out such investigations. 

We are indebted to Dr. C. J. Davisson for 
many helpful suggestions and criticisms. 


6 See Fig. 4 in a note by K. H. Storks and L. H. Germer, 
Phys. Rev. 50, 676 (1936), 
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Linear or Edge Energy and Tension as Related to the Energy of Surface Formation 
and of Vaporization 


WitiiAM D. HARKINS 
George Herbert Jones Laboratory, University of Chicago and Baker Chemical Laboratory, Cornell University 


(Received August 17, 1936) 


The free energy of the edge of a lens of a nonvolatile oil 
on water is found by a consideration of the magnitude of 
the molecular energy involved to be very small, and of the 
order of 3X10 erg per cm. Thus the linear tension of 
the edge should be of the order of 3X10-® dyne. The 
smallness of the energy does not indicate a lack of im- 
portance of the edge, since the molecular energy is higher 


than in the surface. Linear energy and the presence of 
edges are of importance in practically all solids, since these 
are in general built up of crystal grains, and are of specfal 
significance in certain biological systems. Equations which 
give the relation between the free, latent and ‘‘total’”’ 
energy of the edge and the latent heat of vaporization, 
are presented. 





1. INTRODUCTION 


HERE two surfaces meet aa edge is formed. 

This may be considered as a simple edge. 

A more complicated or multiple edge is produced 

by the meeting of a greater number of surfaces or 

interfaces and this is the only type which is stable 

for liquids. A surface may be considered to consist 

of a two-dimensional array of molecules: an edge 
by a linear array. 

Where two or more edges meet, a point or 
angle is formed. Phases, surfaces, edges and 
points may all be indicated by the more general 
term regions. The energy of a certain amount of 
material increases in general as it is moved from 
a phase, or three-dimensional region, to a surface, 
and this increase continues when it is moved to 
an edge and then into a point. The increase of 
energy thus produced may be considered as 
potential energy insofar as there is equipartition 
of the kinetic energy. 

The surface energy of any system is not the 
energy present in its surface regions, but only the 
excess of energy introduced by the presence of the 
surfaces. What may be termed the phasal energy 
is not considered as a part of the surface energy. 
Thus, if any surface region is considered to con- 
sist of a definite amount of material, the energy 
which this material would have in the phase is not 
included in the surface energy. The somewhat 
analogous, but more complicated relations of the 
linear energy are treated in Section 3. 

The topics ‘‘linear wetting” and linear adsorp- 
tion have been considered by Talmud and his 
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associates! but they do not consider the magni- 
tude of the linear energy, which is the subject of 
the present paper. 

It is shown later in this paper that the magni- 
tude of the edge energy is very_small, but that 
this extreme smallness is not due to the smallness 
of the energy per molecule in the edge, but is due 
entirely to the extremely small number of mole- 
cules in one centimeter of an edge, as compared 
with the number in one square centimeter of 
surface or one cubic centimeter of a solid or 
liquid phase. 

The extreme smallness of the linear or point 
energy does not mean that the presence of edges 
and of points is of necessity unimportant. It has 
been found, at least for nonpolar solids, that a 
point on a crystal has a higher solubility and a 
greater vapor pressure than an edge and that in 
an edge these are both larger than in a surface. 

The fact that a given amount of material has a 
higher energy in a point or an edge is of particular 
importance in its effects upon the form of the 
minute subdivisions of biological systems, where 
the form of these subdivisions is of such great 
importance. Thus the importance of edges and of 
points may be said to be greatest in crystalline 
and especially in biological systems. 


2. ORIENTATION OF MOLECULES IN EDGES 


The theory of orientation of molecules in sur- 
faces and interfaces may be applied, in the form 


1D. L. Talmud, B. A. Talmud and S. E. Bresler, Physik. 
Zeits. Sowjetunion 4, 854 (1933); S. E. Bresler, W. W. 
Druschinin and D. L. Talmud, Zeits. f. physik. Chemie 
A164, 389 (1933); D. L. Talmud, B. A. Talmud and S. E. 
Bresler, ibid. A163, 91 (1933). 
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in which it was developed by Hardy and Harkins,” 
to the orientation of the molecules in edges. 

Thus, in a simple edge the least polar groups 
are oriented toward the sharpest part of the edge, 
in order to give the least possible value, under the 
conditions, to the free energy of the edge. The 
more exact relation is, obviously, that the free 
energy of the system is a minimum. In a multiple 
edge the orientation is such as to give all of the 
transitions at the edge as little abruptness as 
possible. That is the sum of all of the transitions 
from phase to edge, and from interface to edge is 
as small as possible in terms of the electrical fields 
involved. Thus active groups, such as NO», CN, 
COOH, COOM, COOR, NH, NHCHs, NCS, 
COR, CHO, I, OH, or groups which contain N, 
S, O, I, or double bonds, turn toward the most 
polar phase, which is often water, or toward 
mercury, or toward a hydrophyllic solid, such as 
titanium dioxide, silicon dioxide, stannic oxide, 
oxides of the metals in general, barium sulphate, 
other salts, etc. 

The hydrocarbon chain of an organic molecule 
in the edge in general orients itself preferably 
toward the vapor phase, or turns toward a non- 
polar phase. Thus at an edge between octyl 
alcohol, water and vapor (with or without air) 
the abruptness of the transition between vapor 
and water, vapor and octyl alcohol, and octyl 
alcohol and water must all be as small, in 
terms of the molecular fields, as is possible. This 
means that the last row of octyl alcohol molecules 
exhibits an orientation with the nonpolar group 
toward the vapor, and the polar group between 
a polar hydroxyl group of the alcohol and a 
molecule of water. 


3. LINEAR OR EDGE ENERGY AND TENSION 


A unit area of the surface of a liquid or solid 
has a certain amount of free surface energy (7) 
and also of what is commonly called total surface 
energy (kh), in excess of that which would be 
present if there were no surface. It has been 
found that in ordinary pure liquids most of the 
surface energy is stored in the outermost layer 
one molecule thick, although some of it pene- 


2 W. B. Hardy, Proc. Roy. Soc. A86, 634 (1911-12) ; 88, 
303 (1913); W. D. Harkins, G. L. Brown, E. C. H. Davies, 
J. Am. Chem. Soc. 39, 354 (1917). 
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trates more deeply. In‘a liquid the free surface 
energy is made evident by a surface tension of the 
same magnitude. The relation between the free 
and the surface energy is given by the equation 


y+l=h, 


in which / is the latent heat and h the Gibbs heat 
function or total surface energy. The entropy (s) 
of the surface is //T which may be shown to be 
equal to —(dy/0T)a, p. 

In this paper all terms in the equations are 
taken to represent one molecule. Thus in Eq. (1) 
y represents the free surface energy for an area of 
surface occupied on the average by one molecule. 

The internal latent heat of vaporization (A) is 
given by Eq. (2) 


A=h+j=ytl4+j, (2) 


in which (j) is the “jumping out energy”’ or 
energy of thermal emission from a surface. 

An edge may be considered as a part of the 
surface or the interfacial region where one radius 
of curvature of at least one of the surfaces is 
extremely small. It is well known that the surface 
energy should increase as the surface curvature 
is increased (radius of curvature decreased) but 
even for ordinary small drops the increase in 
energy is negligible. This is, however, by no 
means the case as the drop approaches molecular 
dimensions. 

For an edge the corresponding equations would 
by analogy be written as 


Yet] =he (1’) 
and A=Ae+Aje=Vetle+je. (2’) 


These quantities do not have, as in Eqs. (1) 
and (2), entirely definite values, but may vary 
between maximum and minimum values, accord- 
ing to the conditions specified. Thus (1) since the 
magnitude of the surface energy terms is depend- 
ent upon the curvature, Eqs. (1) and (2) have 
constant valued terms only because they are 
defined as applicable to plane surfaces alone. An 
edge, however, must be considered as having 
different curvatures, according to its sharpness. 

In addition (2) the values of the terms in the 
equations for linear energy are different depend- 
ing upon whether they refer to the excess in 
energy over that in the phase or in the surface. 
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Thus these equations become : 
Yes +1; =h.., (1!) 
N=Nestjes=Vestlestijes (21) 


referred to the surface, and when referred to a 
volume phase 
Vevtlev =hev, (111) 


A=Nevtjev=YevtletJev- (21) 


Eq. (2!!!) may be used to obtain an idea of the 
order of magnitude of the linear energy and also 
of the linear tension. The value of \ is the same in 
Eq. (2) as in this equation. The effect of the edge 
is to increase the free surface energy y and to 
decrease the energy of thermal emission, so 


Ye>Y, Je<j. (3) 
Also, it may be assumed that /, >/. 


4. SIMPLE EDGE OF AN IDEAL Liguip LENS 
IN VAPOR OR AIR 


Since the latent heat of vaporization (A) is 
greater than the total edge energy (h,) per 
molecule, the first step will be to obtain some 
idea of the order of magnitude of \ for the number 
of molecules assumed to be present in 1 cm of a 
simple liquid edge which represents a highly 
idolized system, since such a lens is not stable. 
Now 


A=sT, 


where s is the entropy of vaporization per mole- 
cule. The Hildebrand-Trouton relation shows 
that for normal liquids the entropy is dependent 
upon the concentration of the vapor alone. Even 
with excessively dilute vapor the value is not over 
20 10-'* erg per molecule per degree. Let it be 
supposed that the edge is so thin as to contain 
only a single row of molecules. Then the number 


TABLE I. Order of magnitude for the linear energy of a sharp 
edge of ethyl ether. 








ERGS PER CM 
X 106 





Free linear energy referred to 


phase vey 2 
Free linear energy referred to 

surface 1 
Total linear energy h.» 4.5 
Linear latent heat dA,» 2.5 








of organic molecules in 1 cm of edge would 
commonly be 2.210’. 

On such a basis the value of \ at 300°K or 
ordinary temperature would be 300X2.2x10? 
X20X10—* erg or 1.310-*° erg. Since y=d 
—l,—j., the value of the free linear energy and 
therefore of the linear tension should be less than 
1.3X10- erg cm, 


Vev <Ney <K=1.3 X 10-5 erg cm-. 


In particular for ether the latent heat of 
vaporization per molecule is 9X10-* erg for 
the number of molecules assumed to be present 
in the edge, while the free surface energy in 
a plane surface is 1.1 10~° erg for this number 
of molecules, and the surface latent heat is 
1.5 10-* erg. 

The value of the free linear energy y-,» would 
then be equal to 1.1 10-°++7,., erg. For a plane 
surface the energy of thermal emission (j) is 
6.0 10-* erg for this number of molecules, and 
the total surface energy is 3X10~ erg. 

A crude method of consideration may be 
adopted for an extremely sharp edge. The mole- 
cules may be considered as a linear array of 
cubes, with only a single linear array below them. 
Thus one-half of each cube would be exposed to 
the vapor, and the total linear energy (h,,.) may 
be taken as 3d; or }X9X107! erg or 4.5X10-° 
erg. If this is divided between latent heat and 
free energy in the same proportions as in a surface 
then the free linear energy y.. is about 210~° 
erg per cm. These results are summarized in 
Table I. 

From the values in Table I, since the linear 
tension is given by the linear free energy, the 
linear or edge tension for a sharp edge of ether 
should be of the order of 2 10~-® dyne. 

The free linear energy and linear tension is 
calculated below by a different method. Consider 
the special case of water whose free surface 
energy is 72.75 ergs per sg. cm or 7.275 X10-" erg 
per molecule. Let the edge be one molecule thick 
and assume that three times as much of each 
molecule is exposed as in the surface. Then the 
free energy would be increased to 14.5 X10-" erg 
above what it is in the surface, or 21.7 K 10~ erg 
above what it is in the interior of the liquid. Since 
there are 3.2107 molecules in 1 cm of the edge, 
the free linear energy is y-s=14.5 or y»=21.7 
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x 10-"+3.2X107, which gives y.;=4.6 or Yer 
=7X10-* erg per cm or the linear tension is 
between 4.6 and 7X10-° dyne. 

At 20° the energy of vaporization of water for 
this number of molecules is 2.0 X 10~5 erg per cm, 
and for this number of molecules the edge energy 
must be less than the energy of vaporization. 


5. VARIATION OF LINEAR ENERGY 
WITH TEMPERATURE 


The variation of the magnitude of the various 
terms in Eq. (1) which gives the latent heat of 
vaporization per molecule from a plane surface 
is represented for liquids by Figs. 1, 2, and 3. 
With the normal liquids the total surface energy 
is almost constant over the ordinary range of 
temperatures for which the liquids are used, but 
decreases extremely rapidly as the critical tem- 
perature is approached. 

The free surface energy decreases rapidly and 
almost linearly with temperature, the latent heat 
of the surface rises rapidly and that of evapora- 
tion and of thermal emission falls rapidly with 
temperature. 

It has been pointed out that curves which 
represent linear energies are not so definite. 
However, if the energy is referred to that of the 
interior of the phase as zero, then the curves have 
the same general form as in the figures. Evidently 
the curve for \ is unchanged, but that for 7.» is 
raised (y..>v), so since this falls to zero at the 
critical temperature, its slope is greater. Since 
l= —T(dy/0T) this causes the latent beat to be 
higher (/.,>/), and thus h,, is higher than h 
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(hey >h). This means that the energy of thermal 
emission is lower for the edge than for the surface 


(Jeo <j). 
6. EFFECT OF THICKNESS OF THE EDGE 


The preceding section indicates that the 
general order of magnitude of the linear tension 
for a sharp edge (angle acute) on a lens of an 
ordinary organic liquid in air is between 10~° and 
10-° dyne. This is based upon the assumption 
that the edge energy is contained largely, not 
wholly, in a linear array one molecule in thick- 
ness. The early work of the writer and his col- 
laborators® gives excellent evidence that most, 
though not all, of the surface energy of a liquid 
is contained in the first monomolecular layer, and 
it is well known that molecular forces from or- 
ganic molecules lose almost their whole effect at 
greater than single molecular distances. 

Suppose that the assumption is advanced that 
the edge is not monomolecular but highly poly- 
molecular, say 100 molecules in thickness. Then 
the magnitude of the linear energy and linear 
tension is by no means multiplied by 100, since 
the total change of direction of the surface in the 
edge is dependent upon the angle in the edge 
alone. The thicker a surface the smaller its 
surface energy, since the transition from one 
phase to the other is thus made less abrupt. 
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The existence of an edge tension for an organic 
liquid or water of 100 times the magnitude esti- 
mated above would be extremely difficult to 


_ explain on the basis of a thick edge, yet such a 


factor of multiplication would increase the 
estimate to only a value less than 0.001 dyne. 

Thus, let it be supposed as above that an 
increase in thickness of the edge does increase the 
free energy of the edge. Then this increase by a 
factor of 100 must be caused by a very thick edge, 
since each molecular thickness cannot now give 
as much as 1/100 the effect in a monomolecular 
edge. On this basis the thickness of the edge 
corresponds to a distance of 10,000 molecules 
along the surface of 45 microns, which is entirely 
too great a distance to represent an edge at all. 
On this basis the linear tension should be much 
less rather than more than 0.001 dyne for organic 
liquids even if it is assumed that an edge is many 
molecules in thickness. 


7. EDGE OR A LENS OF OIL ON WATER 


The calculations given above refer to a liquid 
lens in air, a case seldom or never realized. If a 
lens of oil rests upon the surface of water, then 
the magnitudes obtained above are not very 
greatly changed. Thus, while it may be con- 
sidered that there are now two edges, one of oil 
and one of water, the energy of each edge is much 
lessened by the presence of the other. Also, if one 
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edge is very sharp, the other will in general be 
much less sharp. The result of this is that there 
is a lowering of the linear energy and tension in 
somewhat the same way as the surface energy 
and tension of a liquid are lowered by the pres- 
ence of a second condensed phase. Thus water 
has a surface tension of 72.75 dynes per cm at 20°, 
while that of benzene is 28.88. The interfacial 
tension between the two is 35.0 dynes, which is 
lower than the surface tension of the water, but 
higher than that of the oil. For a petrolatum the 
surface tension was found to be 30.8 and the 
interfacial tension 40.6. Similar values for many 
different organic liquids have been published.‘ 

The surface tension of mercury at 20° is 476 
dynes per cm and the interfacial tension is 375 
toward water, 357 toward benzene, and 304 
toward methyl iodide.’ 

When all of these magnitudes are considered it 
may be estimated that the linear tension of a lens 
of tetradecane on water would be about 3 x 10~° 
dyne if the edge were sharp, or of the general 
order of 2 10~* dyne for the actual edge, which 
is not sharp for the oil phase. It is obvious that 
the angle for the water part of the edge must be 
taken into consideration. 


8. EDGE OF A CRYSTAL 


The magnitude of the linear or edge energy for 
an organic crystal may be expected to be of the 
same general order of magnitude as for the liquid. 
However, with inorganic crystals, especially, 
those of metals, the edge energy is undoubtedly 
very much higher. However, the surface energy 
of sodium chloride, while much higher than that 
of water, is not excessively high. Thus the magni- 
tude of the linear energy of the edge of a cubic 
crystal of this salt is found by Lennard-Jones*® to 
be 4.010-* ergs per cm. The theory is more 
definite in this case, since the theory of the energy 
of an ionic lattice is well developed. The agree- 
ment in the order of magnitude between this and 
the value for a liquid edge as obtained from a 


4W. D. Harkins, G. L. Clark and L. E. Roberts, J. Am. 
Chem. Soc. 42, 700 (1920); W. D. Harkins and Y. C. 
Cheng, J. Am. Chem. Soc. 43, 35 (1921). 

5 W. D. Harkins and E. H. Grafton, J. Am. Chem. Soc. 
42, 2534 (1920); W. D. Harkins and W. W. Ewing, J. Am. 
Chem. Soc. 42, 2539 (1920). 

6 J. E. Lennard-Jones, Proc. Phys. Soc. 43, 461 (1931). 
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consideration of molecular energy magnitudes, is 
as good as could be expected, since the value for 
the liquid edge should be several times lower, and 
it has been found to be about half as high. 
Thus, according to both Lennard-Jones and the 
writer the linear tension should be of the general 
order or 2X10~-° dyne. The value of 6.54 dynes 
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determined by Langmuir’ is in violent disagree- 
ment with this. While the writer has not studied 
the mathematics used by Langmuir, it seems 


obvious that the magnitude which he obtains | 


refers to something else than the true linear 
tension treated in the present paper. 


7 Langmuir, J. Chem. Phys. 1, 762 (1933). 
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The factors which may influence the rates of reactions in solution have been classified and 
discussed with particular reference to the effect on A in the expression Ae~¥/*T for the rate 
constant K. Means by which these different effects may be studied experimentally have been 
discussed. An explanation has been given of the tendency of A and E to change in the same 
direction, and of K and E to change in opposite directions, when comparison is made of the 
rate of a given unimolecular reaction in various solvents. ’ 


NUMBER of recent communications? have 

dealt with the factors affecting reaction 
rates in solution. Since the ideas have not been 
developed very systematically, and since (insofar 
as we are aware) not all possibilities have been 
considered, we feel that it would be useful -to 
present a classification and discussion of these 
factors. 

It is customary in discussing a reaction rate 
constant to write it in the form K=Ae-#/*7, 
where E£, the activation energy, is the difference 
between the average energy of all reacting mole- 
cules, pairs or triplets and that of all molecules, 
pairs or triplets. A bears a somewhat similar 
relation to entropy differences, as has been fully 
discussed in recent publications. We shall refer to 
the reacting molecule, pair or triplet, as an 
activated complex, using the term in the general- 
ized sense of Rice and Gershinowitz.* The term 
reactant or reactants will be used to denote the 
average state of the reacting system. We note the 


1 Present address: Department of Chemistry, University 
of North Carolina, Chapel Hill. 

2 Wynne-Jones and Eyring, J. Chem. Phys. 3, 492 (1935) ; 
Evans and Polanyi, Trans. Faraday Soc. 31, 875 (1935); 
Moelwyn-Hughes and Sherman, J. Chem. Soc. 101 (1936); 
Glasstone, ibid. 723 (1936). 

3 Rice and Gershinowitz, J. Chem. Phys. 3, 482 ws 
Rice, ibid. 4, 54 (1936). 


general rule, now well understood, that anything 
which places a restriction upon the activated 
complex decreases A, while anything which 
places a restriction upon the reactant increases A, 
since (excluding any indirect effect due to a 
concomitant change in the enery) it makes the 
configuration corresponding to the activated 
complex relatively more probable. A number of 
authors, in considering the effect of various fac- 
tors on A and K, for reactions in solution, 
have introduced hypothetical intermediate steps, 
which involve the reacting substances or the 
activated complex in the gas phase or pure liquid 
phase. It is our opinion, however, that very little 
of the inherent complexity of the problem is 
avoided in this way and we prefer to restrict out 
considerations to the actual reacting systems, 
without introduction of hypothetical processes. 

With this understanding, then, we may pro- 
ceed with our classification. We may divide all 
the influences which affect reactions in solution 
into two classes, as follows: 


I. FACTORS WHICH ALSO OBTAIN IN GAS PHASE 


Under this heading we include such matters as 
the orientation or freezing out of rotations into 
vibrations in bimolecular associations. The effect 
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of such factors on A has already been extensively 
discussed. They can also affect E, but this is of 
less interest. These factors operate in the same 
way in solution as in gas phase (or, at least, any 
difference is included in the discussion below) so 
that it is unnecessary to go further into the 
matter in this communication. 


II. FACTORS WHICH OBTAIN ONLY IN SOLUTION 
These may again be divided into two classes: 


1. The internal volume effects 


A solute is not free to move through the whole 
volume of the solution,’ since a large proportion 
of this volume is occupied by solvent molecules. 
Actually only a very small fraction of the volume 
is available. This fact can affect (a) the trans- 
lational motion, (b) the rotational motion of 
solute molecules. In the case of a bimolecular 
reaction, or one of higher order, (a) must always 
operate to increase A, since the translational 
motion takes place in a smaller volume and the 
collision frequency is thus increased. The factor 
(b) may operate in either direction, depending 
upon whether the rotational motion of the acti- 
vated complex or that of all the molecules is the 
more restricted. This in turn depends to a large 
extent on the number of rotational degrees of 
freedom involved in each case. The internal 
volume effect may also influence the energy of 
activation, for we can very well imagine a change 
of internal volume to take place on reaction. This 
will require work against or by the internal 
pressure, with corresponding effect on the energy 
of reaction and hence on the energy of activation 
of either the direct or the reverse reaction or both. 
We should, perhaps, add a few words about the 
possible magnitude of the internal volume effect. 
It certainly will not be as important as other 
effects, but estimates from entropies of vaporiza- 
tion indicate that in some cases the collision 
number in a bimolecular reaction can be in- 
creased by a factor of 20, considering the transla- 
tional effect alone.’ If the rotational entropy of 
the activated complex as compared to that of the 
reactant is also affected and in the same direction 


ase Evans and Polanyi, Trans. Faraday Soc. 31, 891 
5). 

* See, e.g., Rice, J. Chem. Phys. 4, 371 (1936). Evans 
and Polanyi (reference 4) estimated a factor of 100. 
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the increase in A may well be by a factor between 
100 and 500. 


2. Solvation effects 


Solvation means less freedom for both solute 
and solvent molecule. That is to say (since the 
entropy is a measure of the freedom of motion of 
the molecules) if the reactant or the activated 
complex is solvated its entropy (in which we 
include also the entropy of the surrounding mole- 
cules of the liquid, which are really part of the 
reacting system) is lowered. Now a state which 
has a low entropy is an improbable state, so 
solvation will not occur unless the energy of the 
solvated state is less than that of the unsolvated 
state, and sufficiently less so that the probability 
of the solvated state is greater. The result of this 
situation, in the light of our general rule regard- 
ing the effect of the entropy and the probability 
of reactant and activated complex on the rate of 
reaction, is that solvation of the reactants in- 
creases A, increases E, and decreases K (since, by 
the general rule, the statement of the preceding 
sentence is equivalent to saying that the increase 
in E must be more than sufficient to counter- 
balance the increase in A), while solvation of the 
activated complex has just the opposite results.* 

We believe that the above list of solution 
effects is, for all practical purposes, exhaustive, 
provided, at least, that we are dealing with a pure 
solvent and one which does not take part in the 
reaction.* Nevertheless, the factors are suffi- 
ciently varied and complicated, so that one may 
well wonder whether it will ever be possible to 
separate all of them experimentally. Bimolecular 
reactions involving complex organic substances, 
which have been so extensively studied in solu- 
tion, are of the least favorable type from a theo- 
retical point of view, since every one of the above 
factors is operative. It is therefore no wonder 
~ *Since submission of this article for publication an 
article by Evans and Polanyi, Trans. Faraday Soc. 32, 
1333 (esp. p. 1358) (1936), containing some similar con- 
clusions, has appeared. 

6 The above factors are all operative in the case of a 
mixed solvent, but certain complications can arise which 
we do not wish to consider in detail in this communication. 
We may, however, mention what some of these complica- 
tions are. In a mixed solvent we may have solvation by any 
of the constituents of the solvent and compound formation 
between constituents of the solvent ‘including hydration 
of ions when the solvent is an electrolyte solution). It is, 
therefore, necessary to consider equilibria involving the 


various constituents of the solvent, as has been pointed 
out by Professor A. R. Olson in discussions with one of us. 
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that the values of A for such reactions vary over 
a wide range, and it would seem to be only a 
coincidence that the maximum number have a 
value of A corresponding to that expected for 
simple reactions in gas phase, neglecting internal 
volume effects. 

It would appear that the only way to obtain 
independent information on the internal volume 
effect would be to study the bimolecular associa- 
tion of two atoms. In this case only translational 
motion would be involved, orientation effects of 
any kind would be eliminated, and with proper 
selection of reactant and solvent, solvation effects 
might be avoided. Iodine atoms would appear to 
offer nearly the only possibility, but even here the 
experimental difficulties would certainly be great. 

It seems possible, however, to obtain informa- 
tion concerning solvation by comparison of the 
velocity constants for a given reaction in various 
pure solvents,’ and also, when possible, by com- 
parison of the reaction in solution with that in the 
gas phase. It seems plausible to assume that the 
internal volume effect will be similar for most 
solvents, and that the differences in reaction rate 
may be referred chiefly to differences in solvation. 
In fact, for the special case of unimolecular 
reactions which obey the hypothesis of exact 
orientation,® the internal volume effect will be 
negligible as far as A is concerned. For it is 
obvious that it can have no effect upon the trans- 
lational part of the entropy of activation. If the 
hypothesis of exact orientation applies, the rota- 
tional degrees of freedom of the complex are 
identical with those of the reactant, and hence 
the internal volume can have no effect on the 
rotational part of the entropy of activation. 

The unimolecular decomposition of N2O; has 
been studied both in solution and in gas phase.° 


7 It should be pointed out that the effect of varying the 
composition of a mixed solvent is not necessarily similar to 
that of using different pure solvents, as will be clear from 
reference 6. The change of equilibria involving solvent 
molecules with the composition of the solvent may be 
such as to mask effects of the sort considered here. In 
addition, cases in which there is chemical reaction with the 
solvent, changing the nature of the reaction under con- 
sideration, must be excluded. Reaction with the solvent 
presumably occurs in the mutarotation of glucose, for 
example; thus the difference in A and E when the reaction 
takes place in H,O or D.O (Hamill and LaMer, J. Chem. 
Phys. 4, 395 (1936)) is not covered in our discussion. 

8 Rice and Gershinowitz, reference 3, p. 481. 

® Busse and Daniels, J. Am. Chem. Soc. 49, 1257 (1927). 
Eyring and Daniels, ibid. 52, 1473 (1930). For a discussion 
see Moelwyn-Hughes, The Kinetics of Reactions in Solution 
(Oxford, 1933), Chap. VI. 
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There is a number of solvents (so-called ‘‘normal”’ 
solvents) for which K, A and E£ are nearly 
identical with values for the gas phase reaction. 
In other solvents K is smaller, while A and E are 
larger than for the gas phase reaction. 

The fact that a considerable number of liquids 
of various types act as normal solvents for this 
reaction suggests that the increase in activation 
energy due to the effect mentioned under II 
above is either negligible or is compensated by an 
increase of van der Waals forces accompanying 
the formation of the activated complex. Eyring 
and Daniels have suggested that the lower rate 
in abnormal solvents is due to the formation of a 
“complex’”’ between the reactant and solvent 
molecules. The increase in E over that in normal 
solvents may reasonably be attributed to solva- 
tion energy of the reactant. By the above theo- 
retical discussion (see II, 2) this should be ac- 
companied by an increase in A and decrease in K, 
as is actually observed. 

We may now extend the above reasoning to 
other examples of unimolecular reactions, assum- 
ing that the change in activation energy can be 
used as a measure of the amount of solvation. In 
the examples given by Moelwyn-Hughes’ it is 
noted that A and £ for a given reaction in various 
solvents generally increase or decrease together, 
and K varies in the opposite direction. That 
these relations should be true without exception 
is not to be expected, because the change in 
activation energy is by no means an exact meas- 
ure of the change in the amount of solvation of 
the reactant or activated complex.!® 

In the case of nitrogen pentoxide it was 
possible to say that the reactant molecules were 
the solvated ones in abnormal solvents, because 
the rates were lower than the gas phase rates. 
Unless some standard is available, it will not be 


10 The parallelism between E and A has been noted by a 
number of writers in the past (see, e.g., Roginsky and 
Rosenkewitsch, Zeits. f. physik. Chemie B10, 47 (1930); 
ibid. B15, 103 (1931) ; Syrkin, Zeits. f. anorg. Chemie 199, 
28 (1931)), but, as far as we are aware, the explanation 
offered in this paper has not been presented before. We 
feel that it is more reasonable than the explanation of 
Roginsky and Rosenkewitsch, which we believe to be 
scarcely consistent with the more recent theories of chemical 
reaction. Insofar as any parallelism between E and A for 
different reactions exists, it seems likely that the suppost- 
tion of Pollisar (J. Am. Chem. Soc. 54, 3105 (1932)); 
that it is merely accidental, is correct. Whether the rela- 
tionships noted in the data we have been considering are 
to some extent accidental, also, is hard to decide, but the 
matter needs further experimental investigation. 
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possible to tell whether the reactant or the acti- 
vated complex is solvated, but we feel that it is 
probable that it will usually be the reactant. 

The method suggested for detecting solvation 
effects in unimolecular reactions, i.e., comparison 
of rates in various solvents, may also be extended 
to bimolecular reactions, and the same relations 
between FE, A, and K, noted above, should 
obtain.'! The available data” is not entirely con- 
clusive. We may call attention to the addition 
reaction of w-Br-acetophenone and aniline in 
various solvents. In this case it was found that the 
value of E is very low in benzene and chloroform 
as compared to various alcohols, acetone, and 
nitrobenzene. Corresponding to the low values of 
E are low values of A; but they are foo low and, 
contrary to expectations, K is very much lower in 
benzene and chloroform than in the other 
solvents. However, this reaction undoubtedly 

11 However, since the internal volume effect is so much 
more important in bimolecular reactions, its difference in 
different solvents may need some consideration. 


2 For a summary see Holzschmidt, Zeits. f. anorg. 
Chemie 200, 82 (1931). 
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involves a labile bromine, and hence specific 
reaction with the polar solvents (possibly ioniza- 
tion) cannot be excluded. 

Another possibility for investigating solvation 
effects in bimolecular reactions appears to lie in 
the comparison of metatheses involving respec- 
tively an atomic and the corresponding ionic 
reactant.'® As an example we may suggest the two 
optical inversions 


I-+RI-IR+I- I+RI-IR+I1. 


and 


The first of these reactions should involve exten- 
sive solvation, while in most solvents the latter 
should not. It is to be expected that orientation 
and internal volume effects should be practically 
identical for the two reactions. 

One of us (O. K. R.) wishes to acknowledge his 
indebtedness to the seminar of Professor A. R. 
Olson, which contributed greatly to the formation 
of his ideas concerning reactions.in solution. 


13 In this connection see the work by Young and Olson, 
J. Am. Chem. Soc. 58, 1157 (1936), on some ionic reactions. 
Note especially p. 1162. 
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Lattice Energies, Equilibrium Distances, Compressibilities and Characteristic 
Frequencies of Alkali Halide Crystals’ 


Maurice L. HuGGins? 
Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 


Constants in the Born and Mayer exponential expression for the repulsion between two 
ions have been recalculated. Using these, the lattice energies and equilibrium distances between 
ions have been recomputed, for alkali halide crystals. Expressions have been derived and 
tested for the calculation of the compressibility and the characteristic maximum infrared 
frequency. The relation between these two quantities is shown to be approximated by an 


empirical equation proposed by Madelung. 


RECALCULATION OF “‘a,”” LATTICE ENERGIES AND 
Basic RapIil 


HE average potential energy per “‘molecule’”’ 
in a crystal composed of univalent ions 
symmetrically surrounded, as in structures of 
the NaCl and CsCl types, can be expressed?: 4 





‘This paper was presented at the Pittsburgh Meeting 
of the American Chemical Society, September, 1936. 
. ,  aaaaes address, Eastman Kodak Company, Rochester, 


"2M. Born and J. E. Mayer, Zeits. f. Physik 75, 1 (1932). 
ae... L. Huggins and J. E. Mayer, J. Chem. Phys. 1, 643 


by the equation 


ae> C D 
@= ——_ — —_ — —+bc,_Met"+t"--"0) 
ro ro° Fo° 
(1) 
bc__M’ bc..M’ 
+ - e2 (2r-—kero) 4 —_ga(2r+—kero) | 
where 7>=the distance between closest ion 


centers, obtainable from x-ray data® or from the 


5 P, P. Ewald and C. Hermann, Strukturbericht 1913- 
1928 (Akad. Verlagsgesellschaft, Leipzig, 1931). 




















density ; 

e=the electron charge (4.77 X 107!” e.s.u.) 
in the first term and the basis of 
natural logarithms elsewhere ; 

C and D=van der Waals attraction con- 
stants ;* 7 

b=an arbitrary constant, taken for con- 
venience as 10” erg; 

C4, c_. and c,,=factors introduced® to 
allow (roughly) for the dependence 
of the repulsion between two ions on 
their charges; 

a= Madelung’s constant; 

M and M’=the number of nearest unlike 
and like neighbors, respectively ; 

ke=the ratio of the shortest distance 
between like neighbors to that (ro) 
between unlike neighbors; 

r, and r_=“‘basic radii’’ of the positive 
and negative ions, respectively; 
(r++7r_) is the equilibrium distance 
for two ions of opposite sign when 
the total attraction energy has the 
value —bc,_; 

a is a quantity assumed to be constant 
(i.e., independent of r) for a given 
crystal. 


Born and Mayer* computed a for the alkali 
halides from compressibility and thermal ex- 
pansion data, finding it to be roughly the same 
for all, the average value being about 2.90 
(Table III, col. 13).° 











TABLE I. 
1 2 3 4 
LiX NaX, KX, RbX, CsF CsCl, CsBr, CsI 

a 1.7476 1.7476 1.7626 
6s. | 1.315 1.000 1.000 
C44 | 2.00 1.25 1.25 
é.. 10.75 0.75 0.75 
M |6 6 8 
M' |8 8 6 
ke 1.414 1.414 1.155 
ky 2.000 2.000 1.540 
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The lattice energy of an alkali halide crystal 
at the absolute zero is obtained by adding to 
the potential energy given by Eq. (1) a term 
€=—{vmax for the average oscillational energy 
and summing over all the ions. 

Using a=2.90, Mayer and Helmholz” calcu- 
lated lattice energies for these compounds 
(Table III, col. 16) and Huggins and Mayer‘ 
computed a set of basic radii (Table II, col. 2) 
from which the interionic equilibrium distances 
could be accurately calculated (Table III, 
col. 11). 

Wishing to use Eq. (1) and expressions derived 
therefrom in some further work, the writer has 
recalculated a and the lattice energies. This 
seemed advisable for the following reasons: 

1. In their calculations for the lithium halides, 
Born and Mayer and Mayer and Helmholz 
used incorrect values of the coefficients c,— and 
C443 

2. For r+ and r_ they used the Goldschmidt 
ionic radii, which are a poor approximation to 
the true basic radii. 

3. More accurate values of compressibilities, 
coefficients of thermal expansion, van der Waals 
attraction constants and other quantities enter- 
ing into the calculations are available than those 
used before. 

4. Most of the data used in the calculation of 
a are for a temperature of 298°K or thereabouts, 
rather than 273°K, the temperature used by 
Born and Mayer. (The difference in results due 


TABLE II. Basic radii. 














1 2 3 

a=2.90 a=3.00 

10-8 cm 10-§ cm 
Lit 0.475 0.570 
Nat 0.875 0.940 
Kt 1.185 1.235 
Rbt 1.320 1.370 
Cst* 1.455 1.510 
= 1.110 1.050 
cs 1.475 1.435 
pr 1.600 1.560 
1.785 1.750 








6 F, London, Zeits. f. physik. Chemie B11, 222 (1930). 

7J. E. Mayer, J. Chem. Phys. 1, 270 (1933). 

8L. Pauling, Zeits. f. Krist. 67, 377 (1928). 

9a is the reciprocal of the p used in references 3 and 4. 
The change from p to a in this paper is made to conform 
with the usual practice in equations for the energy of a 
diatomic molecule.!® "' Throughout this paper, a is given 
in units of 10° cm™. 





10P, M. Morse, Phys. Rev. 34, 57 (1929). 
( as = L. Huggins, J. Chem. Phys. 3, 473 (1935); 4, 308 
1936). 

12 J. E. Mayer and L. Helmholz, Zeits. f. Physik 75, 19 
(1932). ; 

18M. Born and M. Géppert-Mayer, ‘‘Dynamische Git- 
tertheorie der Kristalle,’"’ Handbuch der Physik, Vol. 24, 
2nd part, second edition (1933), p. 722. 
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to this temperature difference is very small, 
however.) 

The constants and experimental quantities 
used in the first recalculation of a are collected 
in Table I, column 2 of Table II, and columns 3, 
5, 6 to 8, and 17 of Table III. The average a 
obtained was about 3.0, instead of the previous 
2.9. Using this new value, the basic radii were 
then recalculated, the new results being given in 
column 3 of Table II. With these, @ was once 
more computed for each compound, the results 
(Table III, column 14) differing only slightly 
from those obtained with the older radii. 

The lattice energies computed, taking a= 3.00 
and using the corresponding basic radii, are 
compared with those of Mayer and Helmholz 
and with experimental values in columns 15 to 17 
of Table III. 

Interionic equilibrium distances, computed 
with the new a and basic radii, are compared 
with the experimental distances and with those 
computed by Huggins and Mayer in columns 10 
to 12 of Table III. The average deviation 
between the experimental values and_ those 
newly calculated is now only 0.006A 


CALCULATION OF COMPRESSIBILITIES 


Born and Mayer’s’ Eq. (11) can be rewritten as 


Qkiry® Fr, P 


B= 
9°(d*@/dr?) 
with 


T / 0B T s0V 0p 
ronte2() sZ(2) (2) 
B oT - B? V oT P oP T 


2T (aV 
+—(—) . (2b) 
3VX\ OTT p 


k, is the molecular volume (V/N) divided by 1,°. 
(See Table I.) The denominator of Eq. (2a) is 
given by the expression 





(2a) 





dh 
ie = (bc, Ma?)roe2 "++" -—70) 
dr (3a) 
2ae2 42C 72D 
x[1 +5(1-+hade-tren] [Oy 


Yo ro® ro° 








with 
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c_.\ {M’ Coy 
(OHO bg 
Gen. 2M a (3b) 


Calculation of Fr, p from the experimental 
quantities listed in columns 6, 7, 8 and 19 of 
Table III show it (column 18) to have a value 
within a few hundredths of 1.00 in nearly all of 
the compounds with which we are concerned. 
Hence we may write as an approximation 


=< ——______— , | 
. ro?(d*b/dr?) 


Although the two quantities enclosed in square 
brackets in Eq. (3a) are by no means negligible, 
the error caused by the neglect of one is to a 
considerable extent offset by that introduced by 
the neglect of the other. Hence, as a further 
approximation, one can use the relatively simple 
equation 

OR iro 
ne (5) 


4 


The values of 8 calculation from Eqs. (2), (4) 
and (5) (using a=3.00 and the corresponding 
radii) are given in columns 20, 21 and 22 of 
Table III for comparison with the experimental 
results of column 19. The calculated values are 
on the whole low, yet the agreement (between 
columns 19 and 20) is not bad except in those 
cases (LiCl, LiBr, Lil, NaI, CsCl, CsBr, CsI) 
for which the repulsion between IJike ions is 
important. Perhaps Born and Mayer’s assump- 
tion that the same value of a can be used for the 
more distant (like) atoms as for the closest 
(unlike) atoms is not a sufficiently good approxi- 
mation, or the c._, c,, and c__ factors introduced 
by Pauling may be too far off. (Certain theo- 
retical considerations lead the writer to believe 
that these factors, when not equal to unity, 
should be considerably closer to it than assumed 
by Pauling.) 


CALCULATION OF THE CHARACTERISTIC 
MAXIMUM FREQUENCIES 


The characteristic maximum infrared fre- 
quency max is customarily assumed, for an 
alkali halide crystal, to be the frequency of 
oscillation of the lattice of alkali ions, moving 
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TABLE III. 
2 3 4 5 6 7 8 9 
van der Waals Constants (Ref. 7) 7 (37.) » : (2 ) - ~- 5(Sp) r € 
' ff < 4 2 ee op he Ref. , 
Formula ag ae aa a a ao 10-5 deg™ 10-4 deg"! ates 10-2 erg 
LiF 0.8 18 0.6 11 9.28 218 15.3 0.14 
LiCl 2.0 113 2.4 104 12.2% 718 19.8 ( .10) 
LiBr 2.5 183 3.3 190 14.0% 818 24.3 ( .10) 
Lil 3.3 363 5.3 470 16.7% (8) 37.3 ( .08) 
NaF 4.5 46 3.8 31 9.8! (2) 17.5 a 
NaCl 11.2 180 13.9 180 11.0% 7 24.6 07 
NaBr 14.0 271 19.0 300 11.5" 818 25.9 .06 
Nal 19.1 482 31.0 630 13.5" (8) 40.0 05 
KF 19.5 167 21.0 150 10.0% 438 20.2 ( .08) 
KCl 48.0 452 73.0 560 10.1% + a 27.2 .06 
KBr 60.0 605 99.0 800 11.0% 618 31.9 .05 
KI 82.0 924 156.0 1420 12.5" 618 36.7 04 
RbF 31.0 278 40.0 290 (9.5) (5) (21.0) ( .07) 
RbCI 79.0 691 134.0 960 9.85" 7.816 (30.0) 05 
RbBr 99.0 898 180.0 1340 10.4% (8) 35.0 .04 
RbI 135.0 1330 280.0 2240 11.9% (8) 43.0 .03 
CsF 52.0 495 78.0 600 9.5% 10.216 28.4 ( .07) 
CsCl 129.0 1530 250.0 2600 13.6517 8.816 34.0 04 
CsBr 163.0 2070 340.0 3600 13.917 9,116 42.7 .03 
CsI 224.0 2970 520.0 5800 14.617 6.816 49.0 ( .03) 
10 11 12 13 14 15 16 17 
ro a = Lattice Energy per Mole 
exp. H&M H B&M I = M&H H 
Formula | (Ref.'s) | (Ref. 4) | soem | (Ref. 3). | gs tm-1 fe (Ref. 12) m- 
LiF 2.010 2.002 2.002 2.80 3.17 240.1 243.6 
LiCl 2.572 2.557 2.564 2.62 2.87 199.2 200.2 
LiBr 2.745 2.748 2.741 2.62 2.74 188.3 189.5 
Lil 3.000 3.020 3.013 2.17 2.41 174.1 176.1 
NaF 2.310 2.318 2.314 3.09 213.4 215.4 
NaCl 2.814 2.810 2.817 3.07 3.17 181.37° 183.1 183.5 
NaBr 2.981 2.981 2.981 2.99 3.15 174.6 175.5 
Nal 3.231 3.229 3.230 2.60 2.74 163.9 164.3 
KF 2.665 2.685 2.673 3.14 2.98 189.7 192.5 
KCI 3.139 3.127 3.131 3.16 3.14 165.4 167.9 
KBr 3.293 3.286 3.287 3.07 3.09 159.3 161.3 
KI 3.526 3.521 3.525 2.85 2.93 153.87! 150.8 152.4 
RbF 2.815 2.845 2.837 3.09 181.6 183.0 
RbCl 3.270 3.267 3.273 2.81 3.15 160.7 162.0 
RbBr 3.427 3.422 3.427 2.94 3.03 151.3*° 153.5 156.1 
RbI 3.663 3.653 3.663 2.85 3.01 145.3 148.0 
CsF 3.004 2.996 3.004 3.68 173.7 175.7 
CsCl 3.56619 3.572 3.580 3.23 3.11 152.2 153.1 
CsBr 3.713 3.709 3.712 3.06 2.94 146.3 149.6 
CsI 3.95 3.931 3.936 2.82 2.76 141.5?! 139.1 142.5 
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as a unit, relative to the lattice of halide ions. 
Assuming this, we can calculate this frequency 
from the relation 
Vinax = (k/4y) - (6a) 

where yu is related to the masses of the individual 
ions by the equation 

w=m,m_/m,+m_ (6b) 
and k is the coefficient of A (the displacement of 
one lattice relative to the other) in the expansion 
of d@/dA as a power series in A. 

For a crystal of the NaCl type, with the 
oscillations of the two lattices in one of the 
directions of the center lines between closest 
atoms, we can write (neglecting terms which do 
not vary with A and considering only the 6 
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1 1 4 
d= -el + + 
rTo—A rotA (ro? +A?)! 








1 1 + : 
-4 > a 
(ro—A)® (ro+A)® (79? +A?)? J 








1 1 4 7 

—d + + 
(ro—A)$ (ro+A)$ (79? +A?)4) 
+ bc,_[e*” ++r——ro+A) +e (r4+r——ro—A) 
+ gente str—— (ro? 3%) 47), (7) 
(In this equation, c and d are the van der Waals 
attraction constants for a single pair of unlike 
ions.) On differentiating, expanding in powers of 

A, and collecting terms, one obtains 


db/dA =[2bc,-a2(1 — 2a-'rp-"en (+4770) 
— 60cro~§—112dro- 9 JA 











AE RT AONE 











closest ions) +terms in higher powers of A. (8) 
TABLE III.—(Continued.) 
18 19 20 21 22 23 24 25 26 
Fr, P ¥max Cc’ 

ia Ea. Eq. Eq. 5 a Ea. Ea. Eq. 17 
Pe | Oe Se | eee 10-2 barve | 10" barye-t | 10°18 barye? = 10sec | 10sec | 10%em sec? 
LiF 0.72 1.17 0.88 1.22 1.56 9.21 12.18 14.42 125 
LiCl 1.02 3.41 2.36 2.31 3.41 9.37 9.48 1.18 
LiBr 1.03 4.31 2.64 2.57 4.20 8.94 8.33 1.17 
Lil 0.96 6.01 2.75 2.86 5.58 7.57 7.37 1.16 
NaF 0.84 ym | 1.82 2.17 2.01 7.40 8.13 8.30 1.21 
NaCl 1.06 4.26 3.95 3.72 3.58 4,91 6.17 5.76 1.17 
NaBr 1.08 5.08 4.43 4.09 4.12 4.02 5.20 4.85 1.16 
Nal 1.04 7.07 5.06 4.87 5.45 3.51 4.61 4.14 1.15 
KF 0.87 3.30 2.69 3.10 2.77 6.70 6.61 1.18 
KCl 1.03 5.63 5.67 5.53 4,25 4.25 5.06 4.66 1.15 
KBr 1.05 6.70 6.58 6.29 4.89 3.40 4.12 3.71 1.14 
KI 1.04 8.54 7.65 7.32 5.94 2.94 3.67 3.21 1.13 
RbF 1.02 4.1 3.32 3.25 3.06 5.96 5.57 1.17 
RbCl 1.12 6.65 6.55 5.85 4.38 3.54 4.35 3.79 1.14 
RbBr 1.12 7.94 7.64 6.81 5.05 2.63 3.27 2.79 1.14 
RbI 1.10 9.57 9.13 8.27 6.21 2.32 2.78 2.36 1.13 
CsF 1.13 4.25 4.50 3.97 3.78 5.34 5.51 1.16 
CsCl 1.06 5.95 4.80 4.54 4.39 2.94 3.11 3.03 1.43 
CsBr 1.05 7.06 5.72 5.46 4.90 2.24 2.27 2.14 1.42 
CsI 0.98 8.57 6.55 6.67 6.00 1.89 1.77 1.41 



































Note: Values in parenthesis are estimated. 
1 F. A. Henglein, Zeits. f. Elektrochem. 31, 424 (1925). 
16 P, Bridgman, Proc. Am. Acad. 67, 345 (1932). 
17 G. P. Baxter and C. F. Hawkins, J. Am. Chem. Soc. 38, 266 (1916). 
18 International Critical Tables. 
19 E, Brook, I. Oftedahl and A. Pabst, Zeits. f. physik. Chemie B3, 209 (1929). 
*®L. Helmholz and J. E. Mayer, J. Chem. Phys. 2, 245 (1934). 


21 J. E. Mayer, Zeits. f. Physik. 61, 798 (1930). See also reference 11. 
2 R. B. Barnes, Zeits. f. Physik 75, 723 (1932). 
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(The Coulomb terms exactly cancel out.) 
Precisely the same equation is obtained for 
oscillation of the two lattices in other directions 
relative to each other, e.g., in the directions of a 
face diagonal or a body diagonal of the unit cube. 
Substituting the coefficient of A in Eq. (8) 
for k in Eq. (6a), 


Vmax = (272u)—[bc,_a7(1 - 2a! ro- 1) 


X erlr+tr——r0) — 30 cro-§ — 56dro-!9 }}. (9a) 


Correspondingly, for a CsCl type crystal, 


Vnax = (202) bc4_a?(1 — 2a 79-1) ea ++ -—10) 


—40crp-8— (224/3)dro-°]!. (9b) 
The frequencies calculated by means of these 
equations: are compared with the experimental 
frequencies obtained by Barnes” in Table III, 
columns 23 and 24. The calculated results 
average about 20 percent higher than the 
experimental. 

The van der Waals terms in Eqs. (9a) and 
(9b) are relatively small (in most cases less than 
10 percent of the quantity from which they are 
subtracted), hence 


a 


(1 — 2a-— 179-1) ‘a 
row =|— —-~(be,_a°e* (r —wy| : (10) 


2r7u 


Combining this with Eq. (5) for the com- 
pressibility gives the simple relation 


|(— —) “| 
Vmax © - -_* 

27°M up 
This has been used to compute the frequency 
values of column 25, Table III. The calculated 
quantities are about 9 percent higher than the 
experimental, on the average. 


Madelung" has proposed, from dimensional 
considerations, the following relationship : 





(11) 


4 E,. Madelung, Géttingen Nachr. (1909), 100; (1910), 
43 ; Physik. Zeits. 11, 898 (1910) ; M. Born and M. Géppert- 
Mayer, reference 13, p. 656. 
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(M,M_)} 7 
min = | | ons (12) 
(M,+M_)? 


Here Amin is the wave-length corresponding to 
the maximum frequency ymax, M, and M_ are 
the ordinary atomic weights, p is the density, 
and C’ is a constant determined empirically. 
If we substitute 


Amin =Co/Vmaxs (13) 
co being the velocity of light, 
(14) 
(15a) 
 (15b) 


p= m,+m_/kiro’, 
M,. = Nm,, 
and M_=wNm._, 


N being Avogadro’s number, we obtain 


al 
Vmax > — 
C”"?Ni up 


(16) 


This is obviously of the same form as Eq. (11), 
C’ being given by 


2 arco M} 
c~( )( Y= 20-4) (17) 
3N} ki} 


Although not strictly the same for different alkali 
halide crystals, C’ does not vary greatly for those 
of the same structure type, as may be seen from 
column 26 of Table III. 


CONCLUSION 


The results of the lattice energy and interionic 
distance calculations, using Eq. (1) with the 
new values of a, 7, and r_, are believed to be 
even better than those computed earlier. Taking 
the interactions between ions to be as assumed 
in the derivation of that equation, it has been 
shown that the compressibility and the charac- 
teristic maximum infrared frequency can be 
calculated with reasonable accuracy. An approxi- 
mate equation relating the compressibility and 
this maximum frequency has been derived, 
tested, and shown to agree roughly with Made- 
lung’s empirical relationship. 
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A Test of the Third Law of Thermodynamics by Means of Two Crystalline Forms of 
Phosphine. The Heat Capacity, Heat of Vaporization and Vapor Pressure 
of Phosphine. Entropy of the Gas 


C. C. STEPHENSON! AND W. F. GIAUQUE 
Department of Chemistry, University of California, Berkeley, California 


(Received November 16, 1936) 


The heat capacity of solid and liquid phosphine has been 
measured from 15°K to its boiling point, 185.38°K. 
Transitions were found at 30.29, 49.43 and 88.10°K. The 
heats of transition are 19.6, 185.7 and 415.8 cal./mole, 
respectively. The transition at 49.43°K had not previously 
been reported. In order to obtain complete conversion to 
the form of phosphine which is stable below 49.43°K it 
was necessary to keep the unstable form near 40°K for 
several days. A graph showing the rate of conversion to 
the stable form as a function of temperature has been 
given. The transition at 30.29°K occurs in the form of 
phosphine, unstable below 49.43°K. There is also a region 
of abnormally high heat capacity in the supercooled form 
between 30.29 and 35.66°K. At the latter temperature the 
heat capacity decreases abruptly by 15 percent without 
evidence of latent heat. The melting point was found 
to be 139.35°K and the heat of fusion 270.4 cal./mole. 
The heat of vaporization at the boiling point, 185.38°K, 
was determined calorimetrically. The value found was 
3489 cal./mole. The vapor pressure of solid and liquid 
phosphine was measured and the results represented by 
the equations: Solid phosphine 128.67 to 139.35°K; 
logio P(int. cm Hg)=—895.700/T+6.86434. Liquid 
phosphine 139.35 to 185.56°K; logio P(int. cm Hg) 
= — 1027.300/T — 0.0178530T + 0.000029135 T? + 9.73075. 


The entropy of the solid at 49.43°K may be obtained 
either from measurement on the stable or unstable form. 
The values are 8.14 and 8.13 cal./deg. per mole, respec- 
tively. While the extraordinarily close agreement is some- 
what fortuitous, the measurements provide the most 
accurate existing verification of the third law of thermo- 
dynamics by means of two crystalline forms of a substance. 
The molal entropy of the ideal phosphine gas at the boiling 
point was found to be 46.39+0.1 cal./deg. per mole. The 
corresponding value at 298.1°K is 50.35 cal./deg. per mole. 
This may be compared with the value 50.5 calculated by 
Yost with the assistance of spectroscopic data and an 
estimated bond angle. The above entropy values neglect 
the nuclear spin entropy, Rin 16=5.509 cal./deg. per 
mole, which should be added to obtain the absolute 
entropy. The values given are the ones which should be 
used in ordinary thermodynamic calculations. Phosphine 
is known to be a symmetrical top and the two equal 
moments of inertia are known from its band spectrum. 
It has not been found possible to evaluate the other 
moment of inertia in this manner. However, it can be 
calculated from the accurate entropy determination given 
above. The value found is Jc=(7.9+0.8)x10-* gram 
cm?. The P—H distance was found to be 1.45+0.02A and 
the bond angle, 98+5°. 





N this paper we will present the results of a 

low temperature investigation of phosphine. 
During this research a new form of phosphine, 
stable below 49.43°K, was discovered. By means 
of the new modification and the previously 
known form, which has now been shown to be 
unstable at low temperatures, it has been possible 
to obtain what we believe to be the most accurate 
verification of the third law of thermodynamics, 
using two crystalline forms of a substance. 

The entropy of the gas has been evaluated 
from the calorimetric measurements. The en- 
tropy may also be obtained with the assistance 
of spectroscopic data. While the interpretation of 
the band spectrum of phosphine is not as com- 
plete as might be desired, nevertheless it has 
been possible to prove that the third law value 
of the entropy is in agreement with that given 
for the gas by quantum statistics. 


1 Charles A. Coffin Fellow, academic year 1934-35. 


APPARATUS 


The apparatus and methods employed in this 
research are similar to those used in previous 
investigations on condensed gases in this labora- 
tory. The calorimeter is the one which has been 
given the laboratory designation Gold Calo- 
rimeter II.?;* A full description of a similar 
calorimeter has been given by Giauque and 
Egan.‘ The temperatures were measured and 
heat introduced by means of a gold resistance 
thermometer heater. The resistance thermometer 
was calibrated during the measurements against 
copper constantan thermocouple W, which was 
soldered to the bottom of the calorimeter. This 
thermocouple had been compared with a hy- 
drogen gas thermometer,® and oxygen and 


2 Giauque and Wiebe, J. Am. Chem. Soc. 50, 101 (1928). 

3 Blue and Giauque, ibid. 57, 991 (1935). 

4 Giauque and Egan, J. Chem. Phys. 5, 45 (1937). 

5 Giauque, Buffington and Schulze, J. Am. Chem. Soc. 
49, 2343 (1927). 
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hydrogen vapor pressure thermometers,® by 
former investigators in this laboratory. At the 
conclusion of this research, the thermocouple 
was again checked against the triple point of 
hydrogen, hydrogen vapor pressures, the 43.76°K 
transition temperature of oxygen, and the triple 
point and vapor pressures of oxygen. The thermo- 
couple gave temperatures which were low by the 
following amounts: 0.17° at the triple point of 
hydrogen, 13.95°K; 0.15° at the boiling point of 
hydrogen, 20.37°K; 0.13° at the transition tem- 
perature of oxygen, 43.76°K; 0.11° at the triple 
point of oxygen, 54.39°K; 0.01° at the boiling 
point of oxygen, 90.13°K. Suitable corrections 
were applied to the original calibration. 

As is well known, small changes in physical 
state, which may produce a negligible effect on 
the electrical properties of metals at ordinary 
temperatures, commonly produce rather large 
effects at very low temperatures. The small low 
temperature deviations from the original calibra- 
tion which was made nine years before these 
measurements demonstrate that very reliable 
copper constantan thermocouples can be made. 
It may be recalled, however, that Giauque, 
Buffington and Schulze’ tested 58 spools of 
constantan wire at low temperatures, before 
selecting the particular wire from which this 
thermocouple was made. 

It was considered unnecessary to reinvestigate 
the thermocouple above 90°K since not only was 
it in agreement at this temperature but the 
rapidly increasing dE/dT makes small changes 
in e.m.f. relatively less important at higher 
temperatures. 


PREPARATION AND PURITY OF PHOSPHINE 


Phosphine was prepared by the action of 
potassium hydroxide on phosphonium iodide. 
The phosphonium iodide was made from white 
phosphorus, iodine and water by the method of 
Baeyer.’? The reaction is represented by the 
equation 


The yellowish brown product was sublimed into 


6 Giauque, Johnston and Kelley, J. Am. Chem. Soc. 49, 
2367 (1927). 
7 Baeyer, Annalen 155, 266 (1870). 
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the generating flask to give clear colorless 
crystals. 

To the phosphonium iodide, a twenty-five 
percent solution of potassium hydroxide was 
added slowly. The evolved phosphine was washed 
by hydrochloric acid to remove any P2H,, and 
then by a solution of potassium hydroxide to 
remove hydrochloric and hydroiodic acids. 
Finally, the gas was dried by KOH sticks and 
phosphorus pentoxide, and condensed in a liquid- 
air trap protected by a thin metal jacket. 
About 250 cc of liquid phosphine were prepared 
in this manner, and over half of this was dis- 
carded in the fore and end portions of several 
distillations. The solid phosphine was subjected 
to a high vacuum before the product was 
distilled into the calorimeter. 

A series of short heat capacity measurements 
were.made in the temperature region immedi- 
ately below the melting point in order to deter- 
mine the premelting effect due to liquid-soluble 
solid-insoluble impurity. From these measure- 
ments it was calculated that the impurity was 
less than four moles in one hundred thousand 
moles. 

The gas prepared in the above manner was 
always spontaneously inflammable, although the 
moist gas, as it first came from the phosphonium 
iodide, was usually not. When the dry gas was 
moistened by bubbling through a trap containing 
water, it no longer ignited on contact with air. 
Since the dry gas was later shown to be very 
pure, it may be concluded that the pure gas is 
spontaneously inflammable. This confirms the 
similar observation of Ritchie.® 


VAPOR PRESSURE OF PHOSPHINE 


The vapor pressure of phosphine was measured 
by means of a manometer 1.6 cm in diameter. 
The measurements were made by using a 
cathetometer as a comparison instrument in 
connection with a standard meter mounted in 
the manometer case. The pressure observations 
are corrected to international centimeters of 
mercury by means of data given in the /nfer- 
national Critical Tables.* The standard accelera- 


8 Ritchie, Proc. Roy. Soc. London A128, 551 (1930). 
9 International Critical Tables (McGraw-Hill Book Co., 
New York, 1926). 





PROPERTIES OF PHOSPHINE 


TABLE I. Vapor pressure of phosphine. (Boiling Point 
185.38°K, 0°C=273.10°K.) 








T OBs. 
—T CALC. 


P oss. 
—P CALC. 


—0.002 | +0.01 Solid 
.001 0.00 Solid 
.002 | —-.01 Solid 
.000 .000 Triple point 
.002 | — .005 Liquid 
.007 .012 Liquid 
.005 .007 Liquid 
.000 .000 Liquid 
.006 .005 Liquid 
.008 .005 Liquid 
.005 .002 Liquid 
.002 .001 Liquid 
.005 .001 Liquid 


P (int. cm Hg) 
(OBSERVED) 


0.798 
1.312 
1.921 
2.733 
4.477 
6.635 
8.939 
12.380 
17.914 
26.350 
36.674 
50.790 
76.724 


T°K 


128.67 
132.76 
136.10 
139.35 
144.521 
148.960 
152.549 
156.699 
161.724 
167.369 
172.560 
178.043 
185.562 





+++ ++ 


i ae ol 

















tion of gravity was taken as 980.665 cm/sec’. 
The gravitational acceleration for this location, 
979.973 cm/sec?, was taken from the work of 
Sternewarte.'” 

The observations have been represented by 
the following equations: Solid phosphine from 
128.67 to the melting point 139.35°K, 


logio Pcint. em Hg) = —895.700/T + 6.86434. 
Liquid phosphine from 139.35 to 185.56°K, 
logio Print. em Hg) = — 1027.300/T 

—0.01785307 +-0.0000291357?+ 9.73075. 


(1) 


(2) 


The observed and calculated values are given 
in Table I. The boiling point was calculated 
from Eq. (2) to be 185.38°K. 

The vapor pressure of liquid phosphine has 
been measured by Henning and Stock." For 
corresponding pressures, the vapor pressure equa- 
tion of Henning and Stock gives temperatures 
which are higher by the following amounts than 
those given by the above equation: at 150°K, 
0.09°; at 160°K, 0.10°; at 170°K, 0.12°; at 
178°K, 0.18°; at the boiling point, 0.27°. The 
two temperature scales have been shown to be 
consistent to 0.03° from 160°K to 195°K by the 
vapor pressures of HCI.? The higher tempera- 
tures given by Henning and Stock may be due 
to an impurity in their sample. The heat of 
vaporization calculated from the vapor pressure 
equation of this research is in better agreement 
with the calorimetric value than that calculated 

10 Landolt, Bornstein and Roth, Physikalische-chemische 


Tabellen (Verlag Julius Springer, Berlin, 1923). 
" Henning and Stock, Zeits. f. Physik 4, 226 (1921). 


TABLE II. Melting point of phosphine. 








PREs- 
T°K SURE 
THERMO-|(int. cm 
COUPLE Hg) 


T°K 
% RESISTANCE 


/0 
DATE AND TIME MELTED | THERMOMETER 


4/20/35 9:45 p.m. 
int 


10:51 P.M. 9 | 139.321 
4/21/35 12:04 a.m.| Supplied heat 
2:05am.) 12 | 139.333 
2:45 a.m.| Supplied heat 
4:56Aa.M.) 20 | 139.343 
5:26 a.m.| Supplied heat 
8:22am.) 30 | 139.347 
9:08 a.m.| Supplied heat 
11:19am.) 40 | 139.346 
12:01 p.m.| Supplied heat 
2:12pm.) 65 | 139.350 
Accepted values 





Heated into melting 
139.31 | 2.730 
2.730 
2.731 
2.733 


2.735 


139.32 
139.35 
139.35 
139.36 


139.35 


139.35 | 2.733 

















from the equation of Henning and Stock, which 
is evidence that their equation is somewhat in 
error near the boiling point. 


MELTING POINT OF PHOSPHINE 


Table II gives a summary of the observed 
melting point temperatures and pressures. 

Table III contains a summary of measure- 
ments of the melting and boiling point tempera- 
tures for phosphine. 


THE TRANSITION OF PHOSPHINE 


Clusius“ has measured the heat capacity of 
phosphine to the temperatures of liquid hy- 
drogen. He found transitions at 88.5 and 30.3°K. 
He has published his heat capacity results only 
in graphical form. Two somewhat similar regions 
of high heat capacity are shown, one near 35°K 
and the other about 50°K. 

The present investigation confirms the exist- 
ence of the two transitions reported by Clusius 
and the region of high heat capacity at 35°K. 


TABLE III. Melting and boiling point temperatures of 
phosphine. 








MELTING BOILING 
POINT POINT 
°K °K 

140.6 188 
— 186.9 
185.65 


OBSERVER 


Olszewsky"® (1886) 

McIntosh and Steele!* (1904) 
Henning and Stock" (1921) 
Clusius'* (1933) 

This research 





139.7 
139.35 


185.38 








12 Olszewsky, Phil. Mag. (5) 39, 188 (1895). 
13 McIntosh and Steele, Proc. Roy. Soc. 73, 450 (1904). 
14 Clusius, Zeits. f. Elektrochem. 39, 598 (1933). 





However, we have found that another transition 
exists at 49.43°K. The previous observations 
below this temperature refer to an unstable 
modification. The small rise in heat capacity 
shown in Clusius’ graph near 50°K indicates 
that one or two percent of the stable form was 
present during his measurements. All of the 
transitions have been carefully reinvestigated. 


THE TRANSITION TEMPERATURE AT 88.10°K 


From cooling and warming curves it was 
found that phosphine can be superheated or 
supercooled one or two tenths of a degree above 
or below 88.10°K. Table IV contains a summary 
of the measurements. 


THE TRANSITION TEMPERATURE AT 49.43°K AND 
THE RATE OF TRANSFORMATION 


When phosphine is cooled slowly below 
49.43°K it undergoes a transformation to another 
form. If cooled with sufficient rapidity, the high 
temperature form may be completely super- 
cooled and kept at liquid hydrogen temperatures 
with practically no conversion. On warming, the 
supercooled form begins to change spontaneously 
into the stable low temperature form at about 
36°K. The rate of the transformation reaches a 
maximum at about 40°K and becomes too small 
to measure at 47°K. Complete conversion to 
the stable form requires two to three days near 
40°K. Fig. 1 shows the variation of the rate of 


TABLE IV. The transition temperature at 88.10°K. 








T°K 
RESISTANCE 
TIME % TRANSFORMED THERMOMETER THERMOCOUPLE 








0 Cooled approximately half-way through the 
transition 

1:44 l 87.99 | 87.98 
5:00 | Cooled below the transition temperature and 
conversion completed 








5:30] Heat added 
8:12 17 88.20 88.20 
8:23} Heat added 
13:36 33 88.22 88.23 
13:51 | Heat added 
17:20 50 88.25 88.24 
17:28} Heat added 
20:35 67 88.31 88.34 


20:43 | Heat added. Transition complete 

24:34 | Cooled partially back through the transition 
31:30 | 87.99 | 88.03 
Accepted value 88.10+.10° 











Clusius!'4 gives 88.5°K for this temperature. 








152 C. C. STEPHENSON AND W. F. GIAUQUE 





.09 


4 | 7N 
‘— \ 


, = 


30 35 45 50 





DEG.7 MIN. 























40 
TEMPERATURE °K 


Fic. 1, Variation of the rate of transformation of phos- 
phine with temperature ordinates: warming rate of the 
calorimeter in deg. per min. with 50 percent transformed. 


transformation with temperature. This is indi- 
cated approximately by a graph of the warming 
rates at various temperatures. The heat leak 
from the surroundings was negligible during the 
measurements, which were made with between 
50 and 60 percent transformed. 

To determine the transition temperature, 
sufficient energy was added to the low tempera- 
ture form to cause approximately 50 percent 
conversion. The temperatures attained in three 
separate measurements were: 49.44°K, after an 
equilibrium period of 12 hours; 49.40°K, equi- 
librium period, 6 hours; 49.45°K, equilibrium 
period, 7 hours. The average of these values, 
49.43°K, has been taken as the transition 
temperature. 

The transition temperature could not be 
reached by allowing the supercooled phosphine 
to warm spontaneously due to the slow rate of 
conversion; the highest temperature reached in 
this way was 47.3°K. It was possible to heat the 
low temperature form 0.25° above the transition 
temperature before conversion began. 


THE TRANSITION TEMPERATURE AT 30.29°K 


A transition occurred in the unstable form at 
30.29°. No supercooling or superheating could be 
detected at this temperature. Table V contains 
a summary of the temperature measurements. 


THE REGION OF HiGH HEAT CAPACITY 
NEAR 36°K 


Above the transition at 30.29°K the heat 
capacity curve shows an abnormal rise, followed 
by a rather abrupt decrease near 35.66°K. 
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Fic. 2. Cooling rate curve for phosphine near 35.66°K. The 
temperature was 35.66°K at 15.5 minutes on scale. 


Within the limit of accuracy of the observations 
there was no isothermal absorption of heat near 
the region of sudden change and for this reason 
the phenomenon was studied by means of 
cooling and warming rates near 36°K. When the 
rate of cooling or warming is plotted against 
time, a break is found at 35.66°K extending 
over the comparatively small region of +0.03°. 
Fig. 2 shows one of these cooling curves. 

The curve indicates that the heat capacity 
drops by approximately 15 percent within 0.06°, 
and no doubt this region would be smaller if the 
equilibrium were instantaneous. There is no 
indication of an isothermal heat effect, or of 
superheating or supercooling. 


MEASUREMENT OF AMOUNT OF PHOSPHINE 


The amount of phosphine used in the calo- 
rimetric measurements was determined in a 
thermostated bulb of about 5 liters capacity. 
The apparatus and procedure have been de- 
scribed" previously. 

During the experiments, some of the same 
phosphine on which the ealorimetric measure- 
ments were made was used by Long and Gul- 
bransen!® for their data of state measurements 
on this substance. From their data and Ritchie’s® 
density measurements at 0°C we have calcu- 
lated the coefficient, a, of deviation from 
Boyle’s law at one atmosphere and 298.1°K to 
be —0.0074 per atmos where PV = Po V,(1+aP). 
In calculating this value we made use of the 
molecular weight 34.044 rather than the value 





ow and Johnston, J. Am. Chem. Soc. 51, 2300 
a and Gulbransen, J. Am. Chem. Soc. 56, 203 


34.000 obtained by Ritchie and used by Long 
and Gulbransen in presenting their data. Essen- 
tially the calculation assumes that Ritchie’s 
density measurements were correct but that 
the extrapolation of the PV product to zero 
pressure is in error. The molal volume of an 
ideal gas at 0°C=273.10°K and 1 atmos. was 
taken as 22411.5 cc. 2.446 moles of phosphine 
were used for the measurements. 


HEAT CAPACITY OF PHOSPHINE 


The heat capacity results are given in Table 
VI. For the calculation of energy, 1.0004 absolute 
joules were taken equal to 1 international joule 
and 4.185 absolute joules were taken equal to 1 
calorie. The data are shown in Fig. 3. 

Before Series III and IV were taken, the 
calorimeter was kept, in each case, near 40°K 
for 48 hours. Complete conversion was assured 
by the absence of any heat effect due to the 
transition in the unstable form at 30.29°K. 

The heat capacities of the unstable form from 
30.29°K to 35.66°K and the heat of transition 


TABLE V. Transition temperature of phosphine. 








T°K 
RESISTANCE T°K 
% THER- THERMO- 


4/0 
DATE AND TIME TRANSFORMED| MOMETER COUPLE 

















Series 1 
| 
5/7/35 9:04 p.m. Heated into transition 
9:17 P.M. 40 | 30.289 30.25 
9:38 P.M. Heat added 
10:01 P.M. 60 | 30.288 30.25 
10:08 p.m. Heat added 
10:40 p.m. 90 | 30.293 30.27 
Series 2 
| 
5/23/35 9:21 P.M. Heated into transition 
9:48 P.M. 25 | 30.284 30.28 
Series 3 





| 
5/25/35 1:13 a.m. Heated into transition 
1:27 s.M. 65 | 30.290 30.34 





Series 4 





| 
8/7/35 7:00 p.m. Cooled partially 
through transition 


9:30 P.M. | 30.275 | 30.33 
Accepted value | 30.29 
| 














Clusius'* gives 30.3°K for this temperature. 








154 Ce 
at 30.29°K depended to some extent upon the 
rate of cooling between 50 and 35°K. The rate 
of cooling before the measurements on the form 
unstable below 49.43°K is given in Table VII. 

The discussion of the effect of different rates 
of cooling is given later. 

A heat capacity measurement, which has not 
been given in Table VI, was made between 
49.09 and 51.33°K at the end of Series I. This 
measurement was high because it included a 
heat effect due to the transition of a small 
amount of the stable form which was present. 
The data indicate that 1.6 percent of the unstable 
form had gone over to the stable modification. 
Most of this conversion took place during the 
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heat capacity measurements from 36°K to 47°K; 
consequently, the measurements in this region 
for Series I may be in error by one or two 
percent. 

It is estimated from the rate of transformation 
and the time of cooling that below 36°K, the 
high temperature form contained less than 0.2 
percent of the stable modification for Series I, 
less than 0.1 percent for Series V and X, and 
less than 0.4 percent for Series VII. 

The measurements of Series II were made 
after cooling only to 44.7°K. The absence of any 
heat effect at 49.43°K shows that the high 
temperature form underwent no appreciable 
conversion during these measurements. 


TABLE VI. Heat capacity of phosphine. (Molecular weight, 34.044; 0°C=273.10°K.) 








Cp 


CAL./DEG./MOLE SERIES 


C 


P - 
io 4 CAL. /DEG./MOLE SERIES 





Form stable below 49.43°K 





1.076 IV 
1.197 III 
1.562 IV 
1.679 Ill 
2.090 IV 
2.259 Ill 
2.692 IV 
2.847 Ill 
3.300 IV 
3.550 IV 


2.615 
2.694 
2.955 
2.993 
3.175 
3.340 
3.525 
3.404 
3.667 
3.333 


4.169 
4.188 
4.886 
4.893 
5.571 
5.602 
6.314 
6.469 
6.982 


33.31 
33.41 
37.13 
37.29 
41.12 
41.50 
45.10 
45.86 4.291 
47.85 1.517 
49.43 Transition 


3.824 
3.644 
3.835 
4.075 
3.985 
4.432 
3.974 





Form stable above 49.43°K 





2.747 
3.386 
2.930 
3.743 
3.248 
3.898 
2.776 
3.396 
1.803 
2.201 
1.035 


1.987 
2.200 
2.641 
3.054 
3.431 
3.776 
3.992 
4.551 
4.827 
5.040 

5.594 
1.321 


5.809 
Transition 
1.006 8.20 
2.314 8.80 
3.141 8.73 
0.944 8.65 
1.360 9.38 
2.817 10.16 
0.822 10.25 
2.136 10.55 
0.741 11.35 
Abrupt decrease in heat capacity 
0.844 9.61 


9.55 
9.84 
10.18 
10.41 


1.136 
3.018 
3.499 
2.723 





3.255 
2.553 
2.414 
4.217 
4.832 
5.570 
5.203 
5.922 
5.569 
4.841 
Transition 
5.115 
5.389 
5.184 
5.414 
5.193 
4.998 
4.772 
4.573 
4.063 
Melting point 
5.449 


5.293 
5.175 
5.040 
4.919 
4.798 


47.52 
48.66 
51.14 
54.45 
58.19 
63.40 
68.77 
74.33 
80.07 
85.25 
88.10 
92.69 
98.01 
103.50 
109.39 
114.96 
120.28 
125.82 
131.27 
135.72 
139.35 
146.11 
151.69 
157.11 
162.43 
167.64 
172.71 
177.69 4.677 
182.73 4.558 
185.38 Boiling point 
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Fic. 3. Heat capacity in calories per degree per mole of phosphine. 


Points taken from smooth curves through the 
data are summarized in Table VIII. 

Another series of measurements was made 
from 15°K to 50°K after half of the phosphine 
had been allowed to undergo transformation to 
the stable form. These measurements were 
affected by the greater rate of conversion caused 
by the larger amount of the stable form present, 
and it was not possible to obtain results as 
accurate as those given for the other series. The 
results of this series will be discussed below. 

A comparison of the results of Clusius'* with 
those of the present investigation is not possible 
since the preliminary publication of his results 
is in graphical form. 


THE HEAT OF TRANSITION AT 30.29°K 


A summary of the measurements of the heat 
of transition at 30.29°K is given in Table IX. 


The heat of transition was highest for the 
slowest rate of cooling in the above series. When 
the heat of transition is high, the heat capacity 
curve in the region from 30.29 to 35.66°K is 
correspondingly low. No heat capacities were 
taken in Series X, but the rate of cooling was 
the same as in Series V, so the heat capacity of 
that series was assumed in order to make the 
JS C,dT corrections. The small variations men- 
tioned above produce a negligible effect in 
evaluating the entropy. In any case, the effect 
of the variation in the heat of transition is 


TABLE VII. Rate of cooling from 50°K to 35°K. 








SERIES I Vv VII 





Time in minutes 

0 0 
12 6 18 
13 38 
23 62 
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TABLE VIII, Heat capacity of phosphine. Values taken from 
smooth curves through the data. 
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TABLE IX. Heat of transition at 30.29°K. 








TEMPERATURE CORRECTED HEAT AH 
SCpdT  CAL./MOLE 





























Se CG INTERVAL InpUT/ MOLE SERIES 
T= CAL./DEG./MOLE T°K CAL./DEG./MOLE 30.157-30.832 25.06 5.13 19.93 I 
Form stable below 49.43°K 30.094-30.645 23.60 4.25 19.35 xX 
30.036-30.592 23.53 4.17 19.36 xX 
1K ro oo 29.946-30.500 23.19 3.94 19.25 V 
35 2°66 45 631 29.840-30.647 25.93 5.72 20.21 VII 
30 3/56 ; Average value 19.6+0.4 
Form stable above 49.43°K Clusius" gives 19.5 cal./mole. 
15 1.56 90 11.15 
4 — a so TABLE X. Entropy change from 30.29°K to 35.66°K. 
30 6.23 120 11.37 
35 10.90 130 11.53 SERIES I Vv vit 
ro Rpt +4 ; ee Heat of transition/30.29°K 0.658 0.636—(0.667 
50 10.78 160 14:45 30.29°K-35.66°K 1.500 1.559 1.514 
60 11.42 170 14.44 2.158 2.195 2.181 
70 11.79 180 14.46 
80 12.06 185 14.49 











practically compensated by the change in the 
heat capacity. This is‘shown in Table X by 
separate entropy calculations made for three 
different series of measurements over this 
region. 

The general character of the heat capacity 
curves between 25 and 35°K is of the type 
usually associated with molecular rotation in the 
solid state. 

One measurement of the heat of transition 
was made on a mixture of 52 percent of the high 
temperature form and 48 percent of the stable 
form; the relative amounts were determined by 
assuming that the heat capacity of the mixture 
below the transition at 30.29°K is equal to the 
sum of the heat capacities of the two forms. 
The value obtained was 10.34 cal./mole of the 
mixture, or 19.8 cal./mole for 100 percent of the 
high temperature form. 


HEATS OF TRANSITION AND FUSION 


Table XI contains the results of the determina- 
tions of the heat of fusion and heats of transition 
at 49.43°K and 88.10°K. 


THE HEAT OF VAPORIZATION OF PHOSPHINE 


The heat of vaporization was measured by 
adding energy and accepting the gas evolved in 
the five liter measuring bulb referred to above. 
The pressure was kept constant to about a 


TABLE XI. Heat of transition and fusion of phosphine. 
Transition at 49.43°K 























TEMPERATURE CORRECTED HEAT 4H 
INTERVAL Input/MoLE =| fCpdT|_ CAL./MOLE 
49.076-50.209 196.52 11.03 | 185.5 
48.639-51.138 210.04 24.26} 185.8 
Accepted value | 185.7+.2 
Transition at 88.10°K 
87.652-89.625 137.96 22.49} 115.5 
87.555-89.271 135.51 19.74} 115.8 
87.482-89.725 141.46 25.70} 115.8 
87.372-90.009 146.13 30.22 | 115.9 
Accepted value | 115.8+.1 
Fusion at 139.35°K 
137.917-142.576 334.54 64.37 | 270.2 
138.224-141.177 310.75 40.23 | 270.5 
137.891-140.863 310.19 39.55 | 270.6 
Accepted value | 270.4+.2 














Clusius'’ has given 114.0 cal./mole for the heat of transition at 
88.10°K and 267.3 cal./mole for the heat of fusion. 


TABLE XII. Heat of vaporization of phosphine, boiling point, 
185.38°K. 











AMOUNT TIME OF ENERGY AH at 760 MM. 
EVAPORATED, MOLES InpuT, MINUTES CAL./MOLE 
0.20631 45 3486 
0.21263 46 3489 
0.21666 47 ‘4 3486 
0.21840 47 3491 
0.21376 47 3490 


Accepted value 3489+3 
From the vapor pressure equation of 


this research 3472 
From the vapor pressure equation given 
by Henning and Stock" 3408 
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TABLE XIII. The entropy of phosphine, molecular weight 34.044. 








ForM STABLE ABOVE 49.43°K 


ForM STABLE BELOW 49.43°K 











0-15°K, Debye function hcev/k=99 0.495 hev/k=114 0.338 
Graphical 15-30.29°K 2.185 15-49.43°K 4.041 
Transition 19.6/30.29 0.647 185.7/49.43 3.757 
Graphical 30.29-49.43 4.800 — — 
Entropy at 49.43°K 8.13 8.14 
Graphical, 49.43-88.10°K 6.705 

Transition, 115.8/88.10 1.314 

Graphical, 88.10-139.35°K 5.194 

Fusion, 270.4/139.35 1.940 

Graphical, 139.35-185.38°K 4.140 

Vaporization, 3489/185.38 18.82 

Entropy of actual gas at boiling point 46.25 

Correction for gas imperfection +0.14 


Entropy for ideal gas at boiling point 


46.39+0.1 cal. per deg. per mole. 








millimeter of mercury. A summary of the data is 
given in Table XII. 

The heat of vaporization was calculated from 
the vapor pressure equations using the equation 
of state given by Long and Gulbransen.'* This 
equation of state was corrected to the generally 
accepted molecular weight 34.044 rather than 
the value 34.000 used by the above authors. 


THE ENTROPY FROM CALORIMETRIC DATA 


An itemized account of the calculation of the 
entropy of phosphine is given in Table XIII. 
The entropy correction for gas imperfection was 
made by using Berthelot’s equation of state, 
with 7,=324.1°K and P.=64 atmos.” 

The almost exact agreement of the entropies 
at 49.43°K for the two paths is undoubtedly 
fortuitous; the individual values may be in 
error by several hundredths of an entropy unit. 
The close agreement of the two values is an 
excellent verification of the third law of thermo- 
dynamics. 

As mentioned above, one series of heat 
capacity measurements was made from 15 to 
50°K with half each of the stable and unstable 
forms. The entropy value at 49.43°K, obtained 
from these measurements, agrees with the above 
value to 0.1 cal./deg. per mole. As the amount 
converted was measured in terms of the heat 
evolved during the partial transformation, the 
above entropy agreement is of interest in that 


1” International Critical Tables, Vol. 3, p. 248. 


it indicates that the process of conversion of the 
unstable to the stable form is accomplished by 
the formation of distinct phases. Since the 
character of the heat capacity curves of the two 
forms below 49.43°K indicates molecular rotation 
in the unstable form and no molecular rotation 
in the stable form, we had considered the 
possibility that the slow transformation at low 
temperatures might have involved the loss of 
molecular rotation in a gradual and homogeneous 
manner throughout the solid. The random dis- 
tribution accompanying such a process would 
preclude an entropy agreement such as that 
given above. 


THE ENTROPY FROM SPECTROSCOPIC DATA 


Phosphine has a symmetrical pyramidal form, 
in which the principal moments of inertia 
I,4=1ITe¥le. The value I,4=6.221  10-*° gram 
cm? has been made certain by the measurements 
of Wright and Randall.'* The value of Jc has 
not been determined experimentally, but Yost!® 
has calculated a value of 8.2610-*° gram cm? 
by making the plausible assumption that the 
bond angle is 100°. He then calculated Sa9s_.; 
=50.5+1.0 for phosphine. The nuclear spin 
entropy, R In 16=5.509, is not included in this 
value. 

By using the fundamental frequencies, 990, 
1121, and 2327 cm™, observed by Fung and 


18 Wright and Randall, Phys. Rev. 44, 391 (1933). 
19 Yost, J. Chem. Phys. 2, 624 (1934). 








158 Be a 
Barker”® and the frequency 2340 cm calculated 
from valence forces by Howard?! the entropy 
change of phosphine gas from 185.38°K to 
298.1°K is found to be 3.96 cal./deg. per mole. 
This, added to the experimental value Sjs5.33 
=46.39, gives Soo3.1=50.35 cal./deg. per mole, 
which agrees with the value given by Yost. 

It is quite unlikely that any lack of equilibrium 
exists in the solid which would cause a dis- 
crepancy between the calorimetric and spectro- 
scopic values, since the entropies given by the 
two paths below 50°K are in accord. The experi- 


20 Fung and Barker, Phys. Rev. 45, 238 (1934). 
21 Howard, J. Chem. Phys. 3, 207 (1935). 
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mental value of the entropy of phosphine may 
therefore be used with confidence in thermo- 
dynamic calculations. In addition, the constants 
of the phosphine molecule may be obtained from 
the experimental entropy value, the familiar 
Sackur-Tetrode equation, and the value of J,. 
For Ic, the value 7.9 10-°+0.8x10-* gram 
cm? is calculated; for the P—H _ distance, 
1.45+0.02A; and for the bond angle, 98+5°. 
The value found for J¢ definitely eliminates the 
possibility that phosphine is a spherical top, as 
suggested by Fung and Barker. 

We thank Dr. Roy Overstreet for assistance 
with the experimental work. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


New Data on the H.S Spectrum in the Photographic 
Infrared 


The absorption spectrum of hydrogen sulphide gas in the 
region around 10,000A first reported by Cross! has been 
photographed under conditions rather closely approximat- 
ing those of the first investigator. The absorption tube was 
60 feet long. A 500 w projection lamp was used as the 
source. The gas in the tube was under a pressure of about 
1? atmospheres. This pressure is higher than that used by 
Cross as it was found that at lower pressures, the absorp- 
tion was too weak to be measured conveniently. Photo- 
graphs were taken in the first order of a 21 foot concave 
grating giving a dispersion of 2.6A per mm. Eastman 1Q 
and 3Q infrared plates were used with exposures from 2 to 
24 hours. The bands in this region are very weak and only 
the most intense lines appear definitely on the plates. 

The following lines of the band around 10,100 cm~ were 
measured. 


10,168.1 cm™ 10,128.5 
10,160.6 10,123.5 
10,149.3 10,116.9 
10,141.0 10,110.6 
10,135.7 10,103.9 


Reliable measurements of lines beyond 10,168 cm™ were 
not possible due to the presence of the strong water vapor 
band in this region. 

The following lines, also reported by Cross, were the only 
ones observed between 10,013 and 10,216A; this conclusion 
is based on a study of ten spectrograms examined under 
various lighting conditions. 


9983.8 cm™ 9863.4 
9980.9 9851.2 
9977.5 9839.9 
9973.7 9827.3 
9968.8 9813.9 
9964.0 9800.7 
9958.5 9786.2 
9952.4 

9945.7 

9938.7 


While there may be some doubt as to the presence or 
absence of the weak lines measured by Cross, the five strong 
lines that he observed at 9910.8, 9902.9, 9901.8, 9895.5 and 
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9875.5 cm™ with intensities comparable to the other strong 
lines reported, are definitely missing. This result seems to 
disprove the previous! rotational structure analysis. The 
pronounced minimum in the center of the band classes the 
vibration as being along the middle moment of inertia.* 
This has to be parallel to the symmetry plane of the mole- 
cule because otherwise the angle between the two H—S 
bonds would be smaller than 88° 16’ (cf. Cross’ calcula- 
tions) which is highly unlikely from chemical consideration. 

If we assume that the three fundamentals are approxi- 
mately w;3= 1236, w,=2615, and w,=2632 cm™, as seems 
fairly certain now particularly in view of the work of 
Bailey, Thompson and Hale,’ the following parallel type 
bands would be expected in the region in question: 4w,, 
4we, 2we+2w,, also combinations involving fairly high 
multiples of ws. These high multiple ws bands, however, 
have not been observed for any molecule so that we would 
not be justified in expecting them to appear with appre- 
ciable intensities for H2S. The other three bands seem likely 
for the following reason. There seems to be no reason why 
the special type of electric moment observed for H2O 
should hold for H2S, so we would expect w, to appear in the 
infrared (as well as in the Raman effect) together with its 
overtones. The rather complicated structure observed 
around 3.74% makes it likely that the bands there arg a 
superposition of wy; and ws. No resonance between these 
two frequencies is to be expected for symmetry reasons. 
Furthermore, the statement of Bailey, Thompson and Hale 
that no even overtones of w, should occur seems to be a 
misinterpretation of Weizel’s* calculations which apply 
only to perpendicular type bands. In most molecules com- 
binations of the type 2w,-+2w, appear stronger than high 
overtones like 4w, or 4w,. The band around 9900 cm™ might 
thus be due mainly to 2w,+2w, with a possible superposi- 
tion of some other bands. Sufficient evidence is as yet un- 
available to warrant any conclusions as to the classification 
of the band around 10,100 cm. 


FRANK D. CARVIN 


Washington Square College, 
New York University, 
New York, N. Y., 

December 16, 1936. 


1P. C. Cross, Phys. Rev. 47, 7 (1935). 

2? Cf. Figs. 19 and 20 in Dennison's article on polyatomic molecules. 
Rev. Mod. Phys. 3, 280 (1931). 

3 Bailey, Thompson and Hale, J. Chem. Phys. 4, 625 (1936). 

4W. Weizel, Zeits. f. Physik 88, 214 (1934). 














Further Relationships Between Absorption Spectra of 
Rare Earth Solids and Crystal Structure 


The atoms in a solid are so close together that they exert 
powerful electric fields upon one another and thus cause a 
Stark splitting of the energy states of the atoms. The na- 
ture of the splittings will, of course, depend on the states 
involved, the symmetry of the surrounding atoms in the 
solid and on the nature and distances of the surrounding 
atoms. This problem has been treated theorétically by a 
number of investigators notably Bethe and Van Vleck. 
Penney and Schlapp' have treated the specific case of the 
splitting of the basic state of Nd*** in a cubic field. They 
evaluated the constants in their equations so as to give the 
best fit with the experimental work of Gorter and de Haas? 
on the magnetic susceptibility of Nd2(SO4)3s8H2O. We were 
not able to locate the levels where they were predicted by 
that work. We did, however, find levels with the same rela- 
tive separations, but with only about one-third the over-all 
separation. Evaluating the constants in Penney and 
Schlapp’s equations from our levels, we calculated the mag- 
netic susceptibility of the salt at various temperatures and 
obtained results in excellent agreement with the experi- 
mental results of Zernicke and James, S. Meyer, Selwood 
and others but not in agreement with Gorter and de Haas 
or Jackson.’ Further, using the constants as determined for 
Ndt** I calculated the levels for Er+** and found excellent 
agreement with the levels determined spectroscopically by 
Nutting and Meehan.! 

While any rhombic term in the potential equation should 
be easily detected experimentally since it causes splitting 
of a number of the levels, no such splitting was observed and 
we were forced to conclude that no rhombic field existed 
greater than two percent of the cubic field. The Nd*+** 
and Ert+** investigated are in monoclinic crystals. While 
one might expect that the octahedron of oxygens about the 
rare earth ion would give a field predominantly cubic, it 
seems astonishing that a rhombic term should not arise 
from the remaining atoms in the crystal. We are extending 
our investigation to include as many salts as possible. In 
view of the interesting results which we are obtaining it 
seems advisable to report some of them at this time. 

The hydrated sulphates of the rare earths crystallize in 
a homologous series of the form R2(SO4)3-8H2O and are 
completely isomorphic with each other. In these salts, 
which are monoclinic, the rare earth ion is presumably sur- 
rounded by an octahedron of oxygen atoms, four of which 
are part of water molecules and two of which belong to the 
sulphate. Dr. Howe and I have photographed the absorp- 
tion spectrum of Pr2(SO,4)s-8H.2O at various temperatures 
and found that the basic state is split into four levels in 
agreement with theory. Mr. Keller and I have also photo- 
graphed Pr2(SO,4)3-8H2O containing a small amount of 
Nd2(SO,)3-8H20 as an impurity. The basic state of Nd**+t 
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is split exactly the same as in the pure Ndo(SO,)3-8H,0. 
The excited states are also split the same but are shifted to 
the violet by three or four wave numbers. 

The chlorides of the rare earths form a homologous series 
of the form RCIl;-6H2O. The crystals are also monoclinic 
and presumably have the six oxygens of the water mole- 
cules grouped around the rare earth ion in the form of an 
octahedron. One would, therefore, expect that the splitting 
of the states in the chloride and the sulphate would be the 
same. This was found to be true in the case of Gd salts.5 
One might expect it to be true in the case of Nd***, but 
if one examines the homologous series of the chlorides one 
finds that it ends at Nd*** and that Pr+++, Ce+++ and 
Lat*++ form chlorides containing 7H2O. The seventh water 
apparently enters a hole in the lattice brought about by the 
increasing size of the rare earth ion. In Nde2Cl3-6H:20 this 
hole is almost certainly present but is not large enough to 
contain a water molecule. It will, however, permit the 
octahedron to expand into it slightly, and therefore, the 
crystal fields about the rare earth ion will be slightly weaker 
than expected and the over-all splitting of the basic state in 
the chloride will be less than in the sulphate. 

Mr. Hamlin and I have photographed the spectrum of 
NdCl;-6H.O and have found levels at 0, 62 and probably 
at 215 as compared to 0, 77 and 260 cm™, in the sulphate. 
We found no evidence of a rhombic field as the 62 cm™ level 
was not split. Mr. Keller and I have photographed the 
absorption spectra of PrCl;-7H2O containing Nd*** as an 
impurity. Here the water occupies the hole in the lattice 
and we would expect the fields to be stronger. We found that 
the 77 cm™ level was now shifted to 104 cm—.® The excited 
levels of these salts behave in a much more complicated 
manner and cannot be accounted for on the simple theories. 
There are too many lines in the multiplets so that some of 
them have to be of oscillational or vibrational origin. They 
are strongly evident in the case of Pr and are progressively 
less evident as the atomic weight is increased. Beyond Gd 
only a few, if any, can be of that type. We are extending 
these investigations to salts of Ho, Tm and Er and hope to 
have more results available shortly. 

F. H. SPEDDING 


Baker Research Fellow, 
Cornell University, 
Ithaca, New York, 
January 14, 1937. 


1 Phys. Rev. 41, 194 (1932). 

2 Leiden Commu. 278b. ' ali ™ 

3 The experimental results on the magnetic susceptibility of Nd2(SOs)s 
-8H2O are not satisfactory in that they do not agree among themselves. 
The low values obtained by Gorter and de Haas, and Jackson, which we 
do not check, also give an excellent fit with theory so that their relative 
values must be nearly correct. Perhaps the absolute values are in error. 
In any case more experimental work on this salt would be welcome. 

4 Phys. Rev. 50, 574 (1936). 

5 J. Chem. Phys. 5, 33 (1937). : : 

6 Drs. Meehan and Nutting have informed me privately that they 
also found that the lines of the other rare earths are markedly shifted 
when crystallized with Ce2(SO«)s salts which belong to a different crystal 
type. 
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